Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Regularization”

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 235 records · Page 13

A regularized deconvolution method for turbulent closure modeling in implicitly filtered large-eddy simulation

Turbulence closure models that are based on a regularized deconvolution method (RDM) are proposed for application to implicitly filtered Large-Eddy Simulation (LES). This method reconstructs sub-grid scale contributions by approximately recovering filtered quantities through an optimization approach. Physical principles on the boundedness and conservation of reconstructed quantities are preserved by constraining the optimization problem. Further, to account for effects of scales that are not resolved in implicitly filtered LES, a reconstruction method is proposed that combines sub-grid projection and energy-similarity to enrich unresolved scales. With this method, closure models are developed for turbulent scalar fluxes and turbulence–chemistry interaction that are intrinsically consistent with the LES formulation. These closure models are evaluated in application to a piloted turbulent jet flame with inhomogeneous inlet composition, in which combustion is represented by a two-scalar manifold model. The RDM-formulation is benchmarked against simulations obtained using a presumed PDF method and eddy diffusivity closure model. In partially-premixed flame regions, it is observed that RDM improves the prediction of turbulent mixing, and faster grid convergence is observed for results with RDM compared to simulations employing a presumed PDF method.

42 ENGINEERING↗

PINN surrogate of Li-ion battery models for parameter inference, Part II: Regularization and application of the pseudo-2D model

Bayesian parameter inference is useful to improve Li-ion battery diagnostics and can help formulate battery aging models. However, it is computationally intensive and cannot be easily repeated for multiple cycles, multiple operating conditions, or multiple replicate cells. To reduce the computational cost of Bayesian calibration, numerical solvers for physics-based models can be replaced with faster surrogates. A physics-informed neural network (PINN) is developed as a surrogate for the pseudo-2D (P2D) battery model calibration. For the P2D surrogate, additional training regularization was needed as compared to the PINN single-particle model (SPM) developed in Part I. Both the PINN SPM and P2D surrogate models are exercised for parameter inference and compared to data obtained from a direct numerical solution of the governing equations. A parameter inference study highlights the ability to use these PINNs to calibrate scaling parameters for the cathode Li diffusion and the anode exchange current density. By realizing computational speed-ups of ~2250x for the P2D model, as compared to using standard integrating methods, the PINN surrogates enable rapid state-of-health diagnostics. Finally, in the low-data availability scenario, the testing error was estimated to ~2 mV for the SPM surrogate and ~10 mV for the P2D surrogate which could be mitigated with additional data.

25 ENERGY STORAGE↗

Simulation of gas flow in micro-porous media with the regularized lattice Boltzmann method

One primary challenge for prediction of gas flow in the unconventional gas reservoir at the pore-scale such as shale and tight gas reservoirs is the geometric complexity of the micro-porous media. In this paper, a regularized multiple-relaxation-time (MRT) lattice Boltzmann (LB) model is applied to analyze gas flow in 2-dimensional micro-porous medium reconstructed by quartet structure generation set (QSGS) on pore-scale. In this paper, the velocity distribution inside the porous structure is presented and analyzed, and the effects of the porosity and specific surface area on the rarefied gas flow and apparent permeability are examined and investigated. The simulation results indicate that the gas exhibits different flow behaviours at various pressure conditions and the gas permeability is strongly related to the pressure. Finally, the increased porosity or the decreased specific surface area leads to the increase of the gas apparent permeability, and the gas flow is more sensitive to the pore morphological properties at low-pressure conditions.

02 PETROLEUM↗

A regularized-interface method as a unified formulation for simulations of high-pressure multiphase flows

The injection of multi-species fluids into high-pressure and high-temperature environments beyond the species' critical points is commonly found in engineering applications. At these conditions, for immiscible species, both subcritical interfacial dynamics and supercritical mixing can coexist due to variations in temperature around the mixture critical point. The modeling of these complex transcritical phenomena for large-scale configurations is so far not possible. To address this issue, we propose the Regularized-Interface Method (RIM) as a unified formulation that can describe both sub- and supercritical processes as well as the transition between them. The proposed method is derived via filtering of the nanoscale interface-resolving formulation based on van der Waals' linear gradient theory. Thus, this approach allows for the consistent modeling of interfacial dynamics that vanishes at supercritical conditions, while significantly reducing the temporal and spatial resolution constraints of the original nanoscale formulation. The resulting RIM formulation is examined in interface-capturing simulations of sub-, trans-, and supercritical fuel injection processes, involving droplets and jets. Furthermore, these results highlight the importance of resolving spatio-temporal transitions from subcritical interfacial dynamics to supercritical mixing in high-pressure multiphase simulations, in contrast to commonly employed diffused-interface methods, where interfacial dynamics are often neglected.

Interface capturing↗

Experimental and numerical investigations of a two-body floating-point absorber wave energy converter in regular waves

This paper presents experimental and numerical studies on the hydrodynamics of a two-body floating-point absorber (FPA) wave energy converter (WEC) under both extreme and operational wave conditions. In this study, the responses of the WEC in heave, surge, and pitch were evaluated for various regular wave conditions. For extreme condition analysis, we assume the FPA system has a survival mode that locks the power-take-off (PTO) mechanism in extreme waves, and the WEC moves as a single body in this scenario. A series of Reynolds-averaged Navier-Stokes (RANS) simulations was performed for the survival condition analysis, and the results were validated with the measurements from experimental wave tank tests. For the FPA system in operational conditions, both a boundary element method (BEM) and the RANS simulations were used to analyze the motion response and power absorption performance. Additional viscous damping, primarily induced by flow separation and vortex shedding, is included in the BEM simulation as a quadratic drag force proportional to the square of body velocity. The inclusion of viscous drag improves the accuracy of BEM's prediction of the heave response and the power absorption performance of the FPA system, which agrees well with experimental data and the RANS simulation results over a broad range of incident wave periods, except near resonance in large wave height scenarios. Overall, the experimental and numerical results suggest that nonlinear effects, caused by viscous damping and interaction between waves and the FPA, significantly influence the system response and power absorption performance. Nonlinear effects were found to be particularly significant when the wave height was large and the period was near or shorter than the resonant period of the FPA system. The study is expected to be helpful for understanding the nonlinear interaction between waves and the FPA system so that the structure of the FPA can be adequately designed.

16 TIDAL AND WAVE POWER↗

Tuning the microstructures of decavanadate-based supramolecular hybrids via regularly changing the spacers of bis(triazole) ligands

With tuning the ligands from bte, btp, btb to bth, four new decavanadate-based metal–organic hybrid compounds, [Zn(bte)(H{sub 2}O){sub 4}][Zn{sub 2}(bte)(H{sub 2}O){sub 10}](V{sub 10}O{sub 28})·8H{sub 2}O, [Zn{sub 2}(btp){sub 4}(H{sub 2}O){sub 6}](H{sub 2}V{sub 10}O{sub 28})·4H{sub 2}O, [Zn(H{sub 2}O){sub 6}][Zn{sub 2}(btb){sub 2}V{sub 10}O{sub 28}(H{sub 2}O){sub 6}]·4H{sub 2}O, and [Zn{sub 2}(bth)(H{sub 2}O){sub 10}](H{sub 2}V{sub 10}O{sub 28})·6H{sub 2}O (bte=1,2-bis(1,2,4-triazol-1-yl)ethane, btp=1,3-bis(1,2,4-triazol-1-y1)propane, btb=1,4-bis(1,2,4-triazol-1-y1)butane, bth=1,6-bis(1,2,4-triazol-1-y1)hexane), have been synthesized under conventional conditions. The four compounds represent the first examples of decavanadate-based metal–organic hybrids constructed by Zn–bis(triazole) complexes. Their structural analyses show that the four compounds possess different Zn–bis(triazole) structural motifs and various finally structures, which verifies that regular changing the spacers of ligands is an effective strategy to tuning the structures of polyoxometalate-based hybrids. Also, the electrochemical studies show that the compounds have good electrocatalytic activities towards oxidation of nitrite molecules ascribed to V-centers. - Graphical abstract: Four compounds representing the first examples of V{sub 10}O{sub 28}-based hybrids constructed by Zn–bis(triazole) complexes have been synthesized by changing the spacers of the ligands and their electrocatalytic properties have been investigated. - Highlights: • The first examples of V{sub 10}O{sub 28}-based hybrids constructed by Zn-bis(triazole) complexes. • Verifying that changing the spacers of ligands is a strategy to tuning structures. • Showing good electrocatalytic activities toward oxidation of nitrite molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Experimental and numerical comparisons of a dual-flap floating oscillating surge wave energy converter in regular waves

The floating oscillating surge wave energy converter (FOSWEC) is a wave energy converter that was designed, built, and tested to develop an open-access data set for the purpose of numerical model validation. Here, this paper details the experimental testing of the 1:33-scale FOSWEC in a directional wave basin, and compares experimental data to numerical simulations using the wave energy converter simulator (WEC-Sim) open-source code. The FOSWEC consists of a floating platform moving in heave, pitch, and surge, and two pitching flaps. Power is extracted through relative motion between each of the flaps and the platform. The device was designed to constrain different degrees of freedom so that it could be configured into a variety of operating conditions with varying dynamics. The FOSWEC was tested in a range of different conditions including: static offset, free decay, forced oscillation, wave excitation, and dynamic response to regular waves. In this paper, results from the range of experimental tests are presented and compared to numerical simulations using the WEC-Sim code.

16 TIDAL AND WAVE POWER↗

How Well Does Kohn–Sham Regularizer Work for Weakly Correlated Systems?

Kohn–Sham regularizer (KSR) is a differentiable machine learning approach to finding the exchange-correlation functional in Kohn–Sham density functional theory that works for strongly correlated systems. Here we test KSR for a weak correlation. We propose spin-adapted KSR (sKSR) with trainable local, semilocal, and nonlocal approximations found by minimizing density and total energy loss. We assess the atoms-to-molecules generalizability by training on one-dimensional (1D) H, He, Li, Be, and Be 2+ and testing on 1D hydrogen chains, LiH, BeH 2 , and helium hydride complexes. The generalization error from our semilocal approximation is comparable to other differentiable approaches, but our nonlocal functional outperforms any existing machine learning functionals, predicting ground-state energies of test systems with a mean absolute error of 2.7 mH.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Large Methane Emissions From the Pantanal During Rising Water‐Levels Revealed by Regularly Measured Lower Troposphere CH 4 Profiles

Abstract The Pantanal region of Brazil is the largest seasonally flooded tropical grassland and, according to local chamber measurements, a substantial CH 4 source. CH 4 emissions from wetlands have recently become of heightened interest because global atmospheric 13 CH 4 data indicate they may contribute to the resumption of atmospheric CH 4 growth since 2007. We have regularly measured vertical atmospheric profiles for 2 years in the center of the Pantanal with the objectives to obtain an estimate of CH 4 emissions using an atmospheric approach, and provide information about flux seasonality and its relation to controlling factors. Boundary layer‐free troposphere differences observed in the Pantanal are large compared to other wetlands. Total emissions based on a planetary boundary layer budgeting technique are 2.0–2.8 TgCH 4 yr −1 (maximum flux ∼0.4 gCH 4 m −2 d −1 ) while those based on a Bayesian inversion using an atmospheric transport model are ∼3.3 TgCH 4 yr −1 . Compared to recent estimates for Amazonia (∼41 ± 3 TgCH 4 yr −1 , maximum flux ∼0.3 gCH 4 m −2 d −1 ) these emissions are not that large. Our Pantanal data suggest a clear flux seasonality with CH 4 being released in large amounts just after water levels begin to rise again after minimum levels have been reached. CH 4 emissions decline substantially once the maximum water level has been reached. While predictions with prognostic wetland CH 4 emission models agree well with the magnitude of the fluxes, they disagree with the phasing. Our approach shows promise for detecting and understanding longer‐term trends in CH 4 emissions and the potential for future wetlands CH 4 emissions climate feedbacks.

Gloor, M.↗

Confined dense circumstellar material surrounding a regular type II supernova

With the advent of new wide-field, high-cadence optical transient surveys, our understanding of the diversity of core-collapse supernovae has grown tremendously in the last decade. However, the pre-supernova evolution of massive stars, that sets the physical backdrop to these violent events, is theoretically not well understood and difficult to probe observationally. Here we report the discovery of the supernova iPTF 13dqy = SN 2013fs a mere ~3 hr after explosion. Our rapid follow-up observations, which include multiwavelength photometry and extremely early (beginning at ~6 hr post-explosion) spectra, map the distribution of material in the immediate environment (≲ 10 15 cm) of the exploding star and establish that it was surrounded by circumstellar material (CSM) that was ejected during the final ~1 yr prior to explosion at a high rate, around 10 -3 solar masses per year. The complete disappearance of flash-ionised emission lines within the first several days requires that the dense CSM be confined to within ≲10 15 cm, consistent with radio non-detections at 70–100 days. The observations indicate that iPTF 13dqy was a regular Type II SN; thus, the finding that the probable red supergiant (RSG) progenitor of this common explosion ejected material at a highly elevated rate just prior to its demise suggests that pre-supernova instabilities may be common among exploding massive stars.

79 ASTRONOMY AND ASTROPHYSICS↗

Regularized machine learning on molecular graph model explains systematic error in DFT enthalpies

Abstract A major goal of materials research is the discovery of novel and efficient heterogeneous catalysts for various chemical processes. In such studies, the candidate catalyst material is modeled using tens to thousands of chemical species and elementary reactions. Density Functional Theory (DFT) is widely used to calculate the thermochemistry of these species which might be surface species or gas-phase molecules. The use of an approximate exchange correlation functional in the DFT framework introduces an important source of error in such models. This is especially true in the calculation of gas phase molecules whose thermochemistry is calculated using the same planewave basis set as the rest of the surface mechanism. Unfortunately, the nature and magnitude of these errors is unknown for most practical molecules. Here, we investigate the error in the enthalpy of formation for 1676 gaseous species using two different DFT levels of theory and the ‘ground truth values’ obtained from the NIST database. We featurize molecules using graph theory. We use a regularized algorithm to discover a sparse model of the error and identify important molecular fragments that drive this error. The model is robust to rigorous statistical tests and is used to correct DFT thermochemistry, achieving more than an order of magnitude improvement.

36 MATERIALS SCIENCE↗

Model-based iterative reconstruction with adaptive regularization for artifact reduction in electron tomography

Obtaining high-quality 3D reconstructions from electron tomography of crystalline particles embedded in lighter support elements is crucial for various material systems such as catalysts for fuel cell applications. However, significant challenges arise due to the limited tilt range, sparse and low signal-to-noise ratio of the measurements. In addition, small metal particles can cause strong streaking and shading artifacts in the 3D reconstructions when using conventional reconstruction algorithms due to the presence of Bragg diffraction and the large scattering cross-section difference between the materials of the particles and the background support regions. These artifacts lead to errors in the downstream characterization affecting extraction of critical features such as the size of the metal particles, their distribution and the volume of the lighter support regions. In this paper, we present a two-stage algorithm based on metal artifact reduction, utilizing model-based iterative reconstruction methods with adaptive adjustment of regularization parameters. Our approach yields high-quality 3D reconstructions compared to traditional algorithms, accurately capturing both the metal particles as well as the background support. We demonstrate the effectiveness of our algorithm through simulated and experimental bright-field electron tomography data, showing significant improvements in reconstruction quality compared to traditional methods.

97 MATHEMATICS AND COMPUTING↗

Sequential linear integer programming for integer optimal control with total variation regularization

We propose a trust-region method that solves a sequence of linear integer programs to tackle integer optimal control problems regularized with a total variation penalty. The total variation penalty implies that the considered integer control problems admit minimizers. We introduce a local optimality concept for the problem, which arises from the infinite-dimensional perspective. In the case of a one-dimensional domain of the control function, we prove convergence of the iterates produced by our algorithm to points that satisfy first-order stationarity conditions for local optimality. We demonstrate the theoretical findings on a computational example.

97 MATHEMATICS AND COMPUTING↗

A mixed $\ell_1$ regularization approach for sparse simultaneous approximation of parameterized PDEs

We introduce and assess a novel sparse polynomial technique for the simultaneous approximation of parameterized partial differential equations (PDEs) with deterministic and stochastic inputs. Our approach treats the numerical solution as a jointly sparse reconstruction problem through the reformulation of the standard basis pursuit denoising, where the set of jointly sparse vectors is infinite. To achieve global reconstruction of sparse solutions to parameterized elliptic PDEs over both physical and parametric domains, we combine the standard measurement scheme developed for compressed sensing in the context of bounded orthonormal systems with a novel mixed-norm based $\ell_1$ regularization method that exploits both energy and sparsity. Moreover, we are able to prove that, with minimal sample complexity, error estimates comparable to the best $s$-term and quasi-optimal approximations are achievable, while requiring only {\em a priori} bounds on polynomial truncation error with respect to the energy norm.Finally, we perform extensive numerical experiments on several high-dimensional parameterized elliptic PDE models to demonstrate the superior recovery properties of the proposed approach.

97 MATHEMATICS AND COMPUTING↗

Pauli-Villars regularization in nonperturbative Hamiltonian approach on the light front

The advantage of Pauli-Villars regularization in quantum field theory quantized on the light front is explained. Simple examples of scalar λφ{sup 4} field theory and Yukawa-type model are used. We give also an example of nonperturbative calculation in the theory with Pauli-Villars fields, using for that a model of anharmonic oscillator modified by inclusion of ghost variables playing the role similar to Pauli-Villars fields.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Defect reduction in overgrown semi-polar (11-22) GaN on a regularly arrayed micro-rod array template

We demonstrate a great improvement in the crystal quality of our semi-polar (11-22) GaN overgrown on regularly arrayed micro-rod templates fabricated using a combination of industry-matched photolithography and dry-etching techniques. As a result of our micro-rod configuration specially designed, an intrinsic issue on the anisotropic growth rate which is a great challenge in conventional overgrowth technique for semi-polar GaN has been resolved. Transmission electron microscopy measurements show a different mechanism of defect reduction from conventional overgrowth techniques and also demonstrate major advantages of our approach. The dislocations existing in the GaN micro-rods are effectively blocked by both a SiO{sub 2} mask on the top of each GaN micro-rod and lateral growth along the c-direction, where the growth rate along the c-direction is faster than that along any other direction. Basal stacking faults (BSFs) are also effectively impeded, leading to a distribution of BSF-free regions periodically spaced by BSF regions along the [-1-123] direction, in which high and low BSF density areas further show a periodic distribution along the [1-100] direction. Furthermore, a defect reduction model is proposed for further improvement in the crystalline quality of overgrown (11-22) GaN on sapphire.

36 MATERIALS SCIENCE↗

Toroidal regularization of the guiding center Lagrangian

In the Lagrangian theory of guiding center motion, an effective magnetic field B* = B+ (m/e)v ∥ ∇ x b appears prominently in the equations of motion. Because the parallel component of this field can vanish, there is a range of parallel velocities where the Lagrangian guiding center equations of motion are either ill-defined or very badly behaved. Moreover, the velocity dependence of B* greatly complicates the identification of canonical variables and therefore the formulation of symplectic integrators for guiding center dynamics. Here, this letter introduces a simple coordinate transformation that alleviates both these problems simultaneously. In the new coordinates, the Liouville volume element is equal to the toroidal contravariant component of the magnetic field. Consequently, the large-velocity singularity is completely eliminated. Moreover, passing from the new coordinate system to canonical coordinates is extremely simple, even if the magnetic field is devoid of flux surfaces. We demonstrate the utility of this approach in regularizing the guiding center Lagrangian by presenting a new and stable one-step variational integrator for guiding centers moving in arbitrary time-dependent electromagnetic fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Using Poisson-regularized inversion of Bremsstrahlung emission to extract full electron energy distribution functions from x-ray pulse-height detector data

We measured Electron Energy Distribution Functions (EEDFs) from below 200 eV to over 8 keV and spanning five orders-of-magnitude in intensity, produced in a low-power, RF-heated, tandem mirror discharge in the PFRC-II apparatus. The EEDF was obtained from the x-ray energy distribution function (XEDF) using a novel Poisson-regularized spectrum inversion algorithm applied to pulse-height spectra that included both Bremsstrahlung and line emissions. The XEDF was measured using a specially calibrated Amptek Silicon Drift Detector (SDD) pulse-height system with 125 eV FWHM at 5.9 keV. Finally, the algorithm is found to out-perform current leading x-ray inversion algorithms when the error due to counting statistics is high.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗