Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonlinear processes”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A sensitivity analysis of twinning crystal plasticity finite element model using single crystal and poly crystal Zircaloy

The popularity of crystal plasticity finite element method (CPFEM) models is increasing due to their ability to predict the mechanical response of crystalline materials such as metals and metal alloys more accurately than traditional continuum mechanics models. This is since the crystal plasticity models consider the effect of atomic structure, microstructural morphology, and properties of individual grains. These CPFEM models use a large number of material parameters in order to capture the mesoscale physics which comes with the downside of the tedious calibration process. In this paper, a CPFEM code was developed to include the twinning induced grain reorientation and subsequent crystallographic slip for HPC material. The developed code is incorporated in a large-scale, parallelized nonlinear solver WARP3D. Further, a sensitivity analysis with respect to 22 material parameters was then conducted using single crystal and polycrystal representative volume element (RVE) of Zircaloy material. Loading was applied along five different crystallographic orientations for single crystal RVE and along three directions namely, rolling (RD), transverse (TD), and normal (ND) direction for polycrystal RVE. Results obtained from the sensitivity analysis were used for the calibration of material parameters for Zircaloy. Finally, developed code along with calibrated material parameters was used to investigate the effect of the hydride phase formation in Zircaloy which is a typical case observed for nuclear applications. It was found that the volume fraction of the hydride phase has a significant impact on the mechanical properties of Zircaloy.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Marginally stable current sheets in collisionless magnetic reconnection

Noncollisional current sheets that form during the nonlinear development of spontaneous magnetic reconnection are characterized by a small thickness, of the order of the electron skin depth. They can become unstable to the formation of plasmoids, which allows the magnetic reconnection process to reach high reconnection rates. In this work, we investigate the marginal stability conditions for the development of plasmoids when the forming current sheet is purely collisionless and in the presence of a strong guide field. We analyze the geometry that characterizes the reconnecting current sheet, and what promotes its elongation. Once the reconnecting current sheet is formed, we identify the regimes for which it is plasmoid unstable. Finally, our study shows that plasmoids can be obtained, in this context, from current sheets with an aspect ratio much smaller than in the collisional regime, and that the plasma flow channel of the marginally stable current layers maintains an inverse aspect ratio of 0.1.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A set of test problems for nonlinear optimization

The software is a set of test problems for nonlinear optimization algorithms, including subroutines such as linear algebra routines and automatic differentiation algorithms. The test problems come from chemical engineering open literature, and describe optimization tasks related to the design and operation of processes such as carbon capture, Hydrogen production, heat exchange, and distillation.

Parker, Robert↗
DOE OSTI2023

Phonon-informed Neural Thermal Scattering (NeTS) Optimization for Crystalline Graphite and Beryllium Metal

Fast neutrons born from fission lose energy through scattering interactions in the process of slowing-down. As neutrons thermalize to the order of $k$ $b$ $T$ (where $k$ $b$ is the Boltzmann constant, and $T$ is the temperature of the medium), their de Broglie wavelength and energy approaches the order of inter-atomic spacing and quantized lattice vibrations, i.e., phonons. At thermal energies, the thermal scattering law (TSL), i.e., $S$($α, β$), captures crystal binding contributions to the total reaction rate, or cross section. This dimensionless material property describes the energy ($β$) and momentum ($α$) exchanges available in a medium. Currently, $S$($α, β$) is evaluated in the Full Law Analysis Scattering System Hub (FLASSH) code for discrete inputs and stored as ENDF/B File 7 for 0-phonon elastic (MT 2) and n-phonon inelastic (MT 4) processes. Further processing recasts $S$($α, β$) into cumulative distribution functions for sampling post-collision scattering kinematics. In practice, interpolation schemes are employed to access data between tabulated values. An improvement to this juncture of the nuclear data pipeline is supplying cross sections on-the-fly (OTF), as has been developed for the un-resolved resonance region to minimize non-physical interpolation errors. This capability may improve simulation accuracy for accident and transient analyses, where rapidly varying changes in temperature and pressure are difficult to predict beforehand. To do so, deep artificial neural networks (ANNs) can be employed which collapse non-linear, complex data into a lightweight dictionary of neural weights and biases. This has been successfully demonstrated for the hydrogen in light water $S$($α, β$) dataset in the form of a Neural Thermal Scattering (NeTS) module. In this work, the NeTS framework is extended to consider the impact of material-dependent dynamical features on optimal neural pre-processing and architecture design decisions, such as number of neurons per hidden layer, residual skip connections and neural depth. New NeTS modules for crystalline graphite and beryllium metal illuminate a novel correlation between dynamical nonlinearity and optimal neural parametrization when deploying $S$($α, β$) on-the-fly.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Hybrid Bismuth Halide with Rich Polymorphism and Second Harmonic Generation Response

Hybrid structures have emerged as a promising class of optical materials, due to their ability to couple the robustness of inorganic and the tunability of organic compounds. However, their application in nonlinear optics (NLO) remains limited, largely due to underexplored factors that drive the formation of noncentrosymmetric structures and NLO property characterization. In this work, we explore the formation, structural and temperature polymorphism, and optical properties of the (Et 3 NH) 3 Bi 2 Br 9 composition, which crystallizes as either noncentrosymmetric or centrosymmetric polymorph. Structural analysis showed that the alignment of [Bi 2 Br 9 ] 3– units dictates the symmetry of phases, as well as the nonlinear optical properties, with the triclinic polymorph exhibiting a second harmonic generation (SHG) response both in visible and IR regions (1.29 × KH 2 PO 4 and 0.08 × AgGaS 2 ). Thermal analysis reveals polymorphic phase transitions and low melting points, making them melt-processable and ionic liquid candidates.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Ultrafast nonequilibrium dynamics and high-harmonic generation in two-dimensional quantum spin Hall materials

For this work, we develop the theoretical framework of nonequilibrium ultrafast photonics in monolayer quantum spin Hall insulators supporting a multitude of topological states. In these materials, ubiquitous strong light-matter interactions in the femtosecond scale lead to nonadiabatic quantum dynamics, resulting in topology-dependent nonlinear optoelectronic transport phenomena. We investigate the mechanism driving topological Dirac fermions interacting with strong ultrashort light pulses and uncover various experimentally accessible physical quantities that encode fingerprints of the quantum material's topological electronic state from the high-harmonic generated spectrum. Our work sets the theoretical cornerstones to realize the full potential of time-resolved harmonic spectroscopy for understanding nonequilibrium processes in quantum topological systems and identifying topological invariants in two-dimensional quantum spin Hall solid state systems.

2-dimensional systems↗
DOE OSTI2021

Electrocatalytic Reduction of CO 2 to CO over Ag(110) and Cu(211) Modeled by Grand-Canonical Density Functional Theory

We report the results of modeling CO 2 reduction (CO 2 R) to CO over Ag(110) and Cu(211) surfaces at different applied potentials using grand-canonical density functional theory (GC-DFT), a method specifically designed to accurately model electrochemical systems. In addition to demonstrating GC-DFT’s ability to accurately model electrochemical processes, we also compare it with the computational hydrogen electrode (CHE) approach. GC-DFT predicts that the geometries of these reacting systems strongly depend on the applied potential and the Helmholtz free energies vary nonlinearly with the applied potential, which contradicts a central assumption of the CHE approach. The CHE approach neglects the change in the number of electrons on the electrode surface at different applied potentials, which reduces its accuracy as the potential changes from the potential of zero charge. Our results further demonstrate that the grand free energies of the reaction intermediates not only depend on the value of the applied potential but also on the metal surface type, adsorption site, and adsorbate. GC-DFT’s ability to predict the effect of the applied potential on adsorbate geometry enables it to evaluate different possible reaction mechanisms at different applied potentials. For instance, GC-DFT predicts that the first step of CO2R likely switches from proton-coupled electron transfer to sequential electron transfer and then proton transfer at more reducing potentials, a result that cannot be determined by the CHE because it assumes that all electron transfers are coupled to proton transfers and neglects the effect of the applied potential on the adsorbate geometry.

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

Design, tuning, and blackbox optimization of laser systems

Chirped pulse amplification (CPA) and subsequent nonlinear optical (NLO) systems constitute the backbone of myriad advancements in semiconductor manufacturing, communications, biology, defense, and beyond. Accurately and efficiently modeling CPA+NLO-based laser systems is challenging because of the complex coupled processes and diverse simulation frameworks. Our modular start-to-end model unlocks the potential for exciting new optimization and inverse design approaches reliant on data-driven machine learning methods, providing a means to create tailored CPA+NLO systems unattainable with current models. To demonstrate this new, to our knowledge, technical capability, we present a study on the LCLS-II photo-injector laser, representative of a high-power and spectro-temporally non-trivial CPA+NLO system.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Reverse Polarizability of Rare Earth Ions (La 3+ , Gd 3+ , Lu 3+ , Y 3+ ) in Tellurite Glasses and Glass Ceramics for Optical Limiting

All-optical modulation using inherent third-order optical nonlinearity of a medium has garnered considerable interest in photonics and optoelectronics. Herein, nonlinear optical (NLO) properties of tellurite glasses and glass ceramics (GCs) containing four different rare earths (RE = La, Gd, Lu, and Y) have been deliberated in near-infrared regions under an ultrafast regime. The La-based glass exhibits ~10 times higher nonlinear refraction (n 2 ) and absorption (α 2 ) than reported NLO materials. The NLO susceptibility [χ (3) ] trend in the studied glasses is La > Gd > Lu > Y, matching with RE 3+ polarizability. Furthermore, Ln 2 Te 6 O 15 nanocrystallite-embedded transparent GCs exhibit a larger NLO coefficient due to the enhanced local field from oxygen vacancies in crystallites. Interestingly, the trend of χ (3) in GCs follows the sequence of Y > Lu > Gd > La, precisely opposite to the glasses. This observation challenges the general polarizability approach of RE 3+ ions, emphasizing that quadratic hyperpolarizability of RE 3+ is pivotal for NLO properties of GCs. Among the studied matrices, Y-containing GCs showed the lowest optical limiting (OL) threshold (5.4 mJ/cm 2 at 800 nm), much lower than those of the reported NLO materials, suggesting its potential as a femtosecond NIR-laser safety material. A combination of large α 2 and n 2 from the studied matrices indicates their advantage for harmonic generation, potentially aiding in the design of ultrafast signal processing devices.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Jacobian-free Newton–Krylov method for the simulation of non-thermal plasma discharges with high-order time integration and physics-based preconditioning

A preconditioning framework for the numerical simulation of non-thermal streamer discharges is developed using the Jacobian-free Newton-Krylov (JFNK) method. A reduced plasma fluid model is considered, consisting of electrons, one positive ion, one negative ion, and the electrostatic potential. Here, the plasma kinetics model includes ionization, electron-ion recombination, electron attachment, electron detachment, and ion-ion recombination. The governing equations are made dimensionless, discretized in space with finite differences, and integrated in time with a fully implicit method based on high-order backward differentiation formulas. The preconditioning framework is based on a linearized form of the governing equations and physics-based operator splitting. The efficiency of the preconditioning strategy is assessed through two test cases: streamer propagation between parallel plates and an axisymmetric pin-to-pin discharge. The fully implicit approach overcomes traditional restrictions in the time step size due to processes such as electron drift, electron diffusion, and dielectric relaxation. Excellent performance is observed through relevant statistics of the JFNK solver, although the number of linear iterations increases for the pin-to-pin discharge when nonlinear numerical boundary conditions are imposed at the electrodes. Performance studies show scalability with O(100-1000) processors for O(10M) unknowns with ample room for optimization.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Bayesian High-Rank Hankel Matrix Completion for Nonlinear Synchrophasor Data Recovery

Phasor measurement units (PMUs) provide high temporal-resolution synchrophasor measurements for power system monitoring and control. The frequent data quality issues, such as missing and bad data, prevent the incorporation of synchrophasor data in real-time operations. Most existing data-driven data recovery methods assume the power system dynamics can be approximated by a linear dynamical system, and the recovery performance degrades significantly when the power system is experiencing nonlinear dynamics during significant events. Here, this paper proposes a data-driven Bayesian nonlinear synchrophasor data recovery method (Ba-NSDR) that can recover a consecutive time period of simultaneous data losses or errors across all channels, even when the underlying system is highly nonlinear. The idea is to lift the Hankel matrix of the spatial-temporal synchrophasor data to a higher dimension such that the lifted Hankel matrix is low-rank in that space and can be processed with the kernel trick. Our proposed Bayesian method then infers the probabilistic distributions of synchrophasor from the partial observations. Some distinctive features of Ba-NSDR include an uncertainty index to measure the accuracy of the recovery result and the robustness to parameter selections. Our method is verified on both synthetic and recorded event datasets.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Latent Space Dynamics Identification

LaSDI is a data-driven physical simulation software that forms a latent space for a given high-fidelity model and discovers a set of ordinary differential equations for the latent space dynamics. It allows a fast and accurate solution process, which is useful for multi-query decision making applications, such as design optimization and uncertainty quantification. The performance of the LaSDI framework is demonstrated on four different problems, i.e., 1D and 2D Burgers equations, nonlinear heat conduction, and radial advection problems. Both linear and nonlinear compression techniques, such as neural network and proper orthogonal decomposition, are used to form a latent space. A concept of local dynamics identification procedure is introduced to enable a parametric model, which enhances the accuracy level over a given parameter space.

Fries, William↗
DOE OSTI2022

304L Can Crush Validation Studies

Accurate prediction of ductile behavior of structural alloys up to and including failure is essential in component or system failure assessment, which is necessary for nuclear weapons alteration and life extensions programs of Sandia National Laboratories. Modeling such behavior requires computational capabilities to robustly capture strong nonlinearities (geometric and material), rate- dependent and temperature-dependent properties, and ductile failure mechanisms. This study's objective is to validate numerical simulations of a high-deformation crush of a stainless steel can. The process consists of identifying a suitable can geometry and loading conditions, conducting the laboratory testing, developing a high-quality Sierra/SM simulation, and then drawing comparisons between model and measurement to assess the fitness of the simulation in regards to material model (plasticity), finite element model construction, and failure model. Following previous material model calibration, a J 2 plasticity model with a microstructural BCJ failure model is employed to model the test specimen made of 304L stainless steel. Simulated results are verified and validated through mesh and mass-scaling convergence studies, parameter sensitivity studies, and a comparison to experimental data. The converged mesh and degree of mass-scaling are the mesh discretization with 140,372 elements, and a mass scaling with a target time increment of 1.0e-6 seconds and time step scale factor of 0.5, respectively. Results from the coupled thermal-mechanical explicit dynamic analysis are comparable to the experimental data. Simulated global force vs displacement (F/D) response predicts key points such as yield, ultimate, and kinks of the experimental F/D response. Furthermore, the final deformed shape of the can and field data predicted from the analysis are similar to that of the deformed can, as measured by 3D optical CMM scans and DIC data from the experiment.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Evaluation of Reduced Order Model for HT-9 Creep and Modifications to Current HT-9 Creep Model in BISON

In order to leverage existing reduced order models (ROM) for the modeling of LANL-developed HT-9 cladding constitutive behavior, mechanistic-based models have been implemented into BISON. It is posited that the increased fidelity of the ROM will show improvement over engineering-scale models in terms of comparison against experimental measurements. INL will assist LANL in implementing an HT9 ROM in BISON. Once the HT-9 mechanistic constitutive model can be leveraged in nuclear performance simulations, the numerical results of fast reactor models will be analyzed and compared against the Fuels Irradiation & Physics Database (FIPD) and separate effects tests to attempt a validation process. INL and LANL will work together to perform any necessary improvement that is identified during the implementation and use processes. Successful completion of the milestone will enable the general use of validated mechanistic ROMs for HT-9 cladding in BISON. This report provides a brief introduction to the Los Alamos Reduced Order Model Applied to Nonlinear Constitutive Equations (LAROMANCE) code, mechanical testing used to illustrate its creep and plastic deformation predictions for HT-9 cladding, and a comparison with existing models implemented into BISON. Current models have been modified to more appropriately account for primary thermal creep and results for corresponding testing are included in this report.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Instrument Control Software for Spectrophotometers

SRNL has developed instrument control software for diode array spectrophotometers that can be used for real time monitoring of dynamic systems. This software has evolved through several iterations as instruments have been developed for use in processing facilities at Savannah River Site. There are several noteworthy features of the software. The wavelength calibration is checked and updated in real time against reference emission lines that are designed into the instrument light source. Nonlinearities and offsets associated with the intensity response of the instruments are characterized, and absorbance spectra are appropriately corrected. These actions are necessary for accurate calculation of an absorbance spectrum from the intensity responses of two individual diode array detectors. The corrections also result in the spectra obtained for a particular sample by two (or more) different diode array detector pairs being nearly identical. This property assures consistent measurements at multiple locations with different instruments. It also permits the use of a single analyte calibration across many instruments without modification and without loss of accuracy

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Improving snow albedo modeling in the E3SM land model (version 2.0) and assessing its impacts on snow and surface fluxes over the Tibetan Plateau

With the highest albedo of the land surface, snow plays a vital role in Earth's surface energy budget and water cycle. Snow albedo is primarily controlled by snow grain properties (e.g., size and shape) and light-absorbing particles (LAPs) such as black carbon (BC) and dust. The mixing state of LAPs in snow also has impacts on LAP-induced snow albedo reduction and surface radiative forcing (RF). However, most land surface models assume that snow grain shape is spherical and LAPs are externally mixed with the snow grains. This study improves the snow radiative transfer model in the Energy Exascale Earth System Model version 2.0 (E3SM v2.0) Land Model (ELM v2.0) by considering non-spherical snow grain shapes (i.e., spheroid, hexagonal plate, and Koch snowflake) and internal mixing of dust–snow, and it systematically evaluates the impacts on the surface energy budget and water cycle over the Tibetan Plateau (TP). A series of ELM simulations with different treatments of snow grain shape, mixing state of BC–snow and dust–snow, and sub-grid topographic effects (TOP) on solar radiation are performed. Compared with two remote sensing snow products derived from the Moderate Resolution Imaging Spectroradiometer, the control ELM simulation (ELM_Control) with the default configurations of spherical snow grain shape, internal mixing of BC–snow, external mixing of dust–snow, and without TOP as well as the ELM simulation with new model features (ELM_New) can both capture the overall snow distribution reasonably. Additionally, ELM_New overall shows smaller biases in snow cover fraction than ELM_Control in spring when snowmelt is important for water management. The estimated LAP-induced RF in ELM_New ranges from 0 to 19.3 W m -2 with the area-weighted average value of 1.5 W m -2 that is comparable to the reported values in existing studies. The Koch snowflake shape, among other non-spherical shapes, shows the largest difference from the spherical shape in spring when snow processes related to the surface energy budget and water cycle have high importance. The impacts of the mixing state of LAP in snow are smaller than the shape effects and depend on snow grain shape. Compared to external mixing, internal mixing of LAP–snow can lead to larger snow albedo reduction and snowmelt, which further affect the surface energy budget and water cycle. The individual contributions of non-spherical snow shape, mixing state of LAP–snow, and local topography impacts on the snow and surface fluxes have different signs and magnitudes, and their combined effects may be negative or positive due to complex and nonlinear interactions among the factors. Overall, the changes in net solar radiation in spring due to individual and combined effects range from -28.6 to 16.9 W m -2 and -29.7 to 12.2 W m -2 , respectively. This study advances understanding of the role of snow grain shape and mixing state of LAP–snow in land surface processes and offers guidance for improving snow simulations and RF estimates in Earth system models under climate change.

58 GEOSCIENCES↗
DOE OSTI2025

Single-shot electrical detection of short-wavelength magnon pulse transmission in a magnonic thin-film waveguide

Abstract The advance of magnon spintronics requires understanding of time-domain magnon pulse transmission in order to develop high-speed information processing protocols. In this work, we demonstrate single-shot electrical detection of narrow-band magnon pulse transmission in a yttrium iron garnet thin-film delay line. The high signal-to-background ratio of magnon transmission band allows us to directly probe the magnon transmission electrically using a fast oscilloscope and to study its spectral evolution using Fast Fourier Transform (FFT) of the time-domain transmitted signal. At elevated input power, we show a magnon transmission reduction and a spectral distortion, which can be understood by the nonlinear magnon excitation in the transmission band defined by the antenna geometry. In addition, we also find that the higher- (lower-) frequency magnon spectral component exhibits a lower (higher) magnon group velocity, showing a dispersion agreeing with the Damon-Eshbach dependence. Our results provide important guidance of magnon pulse engineering for their applications in spin wave computing and coherent magnon information processing.

Song, Moojune↗
DOE OSTI2022

A 1D–0D–3D coupled model for simulating blood flow and transport processes in breast tissue

In this work, we present mixed dimensional models for simulating blood flow and transport processes in breast tissue and the vascular tree supplying it. These processes are considered, to start from the aortic inlet to the capillaries and tissue of the breast. Large variations in biophysical properties and flow conditions exist in this system necessitating the use of different flow models for different geometries and flow regimes. In total, we consider four different model types. First, a system of 1D nonlinear hyperbolic partial differential equations (PDEs) is considered to simulate blood flow in larger arteries with highly elastic vessel walls. Second, we assign 1D linearized hyperbolic PDEs to model the smaller arteries with stiffer vessel walls. The third model type consists of ODE systems (0D models). It is used to model the arterioles and peripheral circulation. Finally, homogenized 3D porous media models are considered to simulate flow and transport in capillaries and tissue within the breast volume. Sink terms are used to account for the influence of the venous and lymphatic systems. Combining the four model types, we obtain two different 1D–0D–3D coupled models for simulating blood flow and transport processes: The first model results in a fully coupled 1D–0D–3D model covering the complete path from the aorta to the breast combining a generic arterial network with a patient specific breast network and geometry. The second model is a reduced one based on the separation of the generic and patient specific parts. The information from a calibrated fully coupled model is used as inflow condition for the patient specific sub-model allowing a significant computational cost reduction. Several numerical experiments are conducted to calibrate the generic model parameters and to demonstrate realistic flow simulations compared to existing data on blood flow in the human breast and vascular system. Moreover, we use two different breast vasculature and tissue data sets to illustrate the robustness of our reduced sub-model approach.

multiscale blood flow models↗