Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical differentiation”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

Numerical studies of motion of vortex filaments - Implementing the asymptotic analysis

A computational code is developed for the integro-differential equations governing the motion of the centerlines of vortex filaments submerged in a background potential flow. These equations, which are derived from the method of matched asymptotic analysis, include the effect of the decaying large-magnitude circumferential and axial velocity components in the vortical cores. Numerical examples are presented to assess the effect of a large axial velocity and that of nonsimilar initial profiles in the vortical cores. The initial configurations of the filaments are chosen so as to fulfill the basic assumption of the asymptotic analysis, which is that the effective vortical core size is much smaller than all the other length scales in the flowfield, e.g., the radius of curvature and the interfilament distance. The computations are continued until the basic assumption is no longer valid, that is when the merging or intersection of filaments has begun. A classification of the various types of local or global merging or intersection of filaments is made and demonstrated by numerical examples. It is then shown that the asymptotic solution not only provides the initial data but also can be used to formulate the appropriate boundary conditions for the numerical solution of a merged region.

Liu, C. H.↗
NASA NTRS1973

Numerical integration of shell equations using the field method.

The 'field method' for the numerical solution of even-order linear boundary-value problems in ordinary differential equations is formulated. This method converts the boundary-value problem into two successive initial-value problems, which may be solved by standard forward integration techniques. The method has been implemented in a computer program to calculate the static response of ring-stiffened branched shells of revolution to asymmetric loads. For such problems the field method eliminates the well-known numerical problem of 'long subintervals' and also executes operations significantly faster than other numerical integration methods.

Cohen, G. A.↗
DOE OSTI2021

Long noncoding RNA MEG3 expressed in human dental pulp regulates LPS-Induced inflammation and odontogenic differentiation in pulpitis

Highlights: • LncRNA MEG3, upregulated in inflamed dental pulp, is involved in the inflammation and regeneration in pulpitis. • The expression of lncRNA MEG3 was semiquantitatively detected in the pulp tissue by the innovative RNAscope® technique. • MEG3 downregulation can inhibit inflammatory cytokines secretion and also promote the odontogenic differentiation of hDPCs. • LncRNA MEG3 can be a potential therapeutic target in the inflammation and regeneration of dentin-pulp complex in pulpitis. Pulpitis refers to inflammation of the inner pulp by invading microbes, and tissue repair occurs due to odontogenic differentiation of human dental pulp cells (hDPCs) with multidifferentiation potential. Long noncoding RNAs (lncRNAs) can modulate numerous pathological and biological processes; however, the role of lncRNAs in the inflammation and regeneration of the dentin-pulp complex in pulpitis is unclear. Here, we performed high-throughput sequencing to identify differentially expressed lncRNAs between human normal and inflamed pulp and concluded that lncMEG3 (lncRNA maternally expressed gene 3, MEG3) was signicantly upregulated in both inflamed pulp and LPS-treated hDPCs. MEG3 expression in the pulp tissue was detected using the RNAscope® technique. RNA pulldown assays identified the MEG3-interacting proteins and the potential mechanisms. With MEG3 knockdown, we investigated the role of MEG3 in the secretion of inflammatory cytokines in LPS-treated hDPCs and odontogenic differentiation of hDPCs. MEG3 downregulation inhibited the secretion of TNF-α, IL-1β and IL-6 in LPS-treated hDPCs, and the p38/MAPK signaling pathway may be related to this effect. MEG3 knockdown promoted odontogenic differentiation of hDPCs by regulating the Wnt/β-catenin signaling pathway. Our study suggested that MEG3 has a negative effect on inflammation and regeneration of the dentin-pulp complex in pulpitis.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Modeling betatron radiation using particle-in-cell codes for plasma wakefield accelerator diagnostics

The analysis of plasma wakefield acceleration experimental measurements, particularly in the characterization of photons emitted through the betatron radiation mechanism, requires the development of accurate numerical models. These computational models are crucial for supporting modern instrumentation designed to measure the single-shot, double-differential angular-energy radiation spectra resulting from interactions between beams and plasmas. Motivated by the needs of such applications, this paper presents detailed numerical models of betatron radiation generated in beam-plasma acceleration experiments. These models are based on the integration of the Liénard-Wiechert (LW) potentials, applied to computed particle trajectories. The particle trajectories are generated using three distinct methods: first, by tracking particles through idealized fields in the blowout regime of PWFA; second, by obtaining trajectories using the fast quasistatic particle-in-cell (PIC) code quickpic; and third, obtaining trajectories from the fully self-consistent PIC code osiris. To ensure the accuracy and reliability of these models, the paper includes various benchmark tests using analytical expressions, as well as employing the PIC code epoch, which takes an alternative approach by using a Monte Carlo quantum electrodynamics (QED)-based radiation model. Additionally, the paper presents simulations of the expected experimental betatron radiation spectra, taking into account parameters relevant to PWFA and plasma photocathode experiments at the SLAC FACET-II facility.

Yadav, M. [University of California, Los Angeles, ↗
DOE OSTI2026

Higher-order LaSDI: Reduced order modeling with multiple time derivatives

Solving complex partial differential equations (PDEs) is essential across scientific disciplines but often requires numerical models that can be prohibitively expensive in time-sensitive applications. Reduced-order models (ROMs) address this challenge by exploiting low-dimensional structure to create fast approximations. The Latent Space Dynamics Identification (LaSDI) framework has demonstrated success in learning ROMs for parameterized PDE families, but remains limited to first-order systems. Here, in this paper, we propose Higher-Order LaSDI (HLaSDI), which extends the LaSDI framework to PDEs with arbitrary order of time derivatives. This generalization significantly expands the applicability of LaSDI-based methods to systems previously outside their scope, including hyperbolic PDEs. We demonstrate HLaSDI’s accuracy and efficiency on several linear and nonlinear benchmark problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Effect of time and spatial domains on monolayer 2D material interface thermal conductance measurement using ns ET-Raman

Raman-based approaches to characterizing thermal properties of nanoscale and microscale materials have proven to be powerful experimental methods. Here, this work explores the potential of energy-transport state-resolved Raman (ET-Raman) for the precise determination of interfacial thermal conductance (G) between monolayer 2D materials and their substrates. ET-Raman employs a continuous wave laser and nanosecond (ns) pulsed laser to induce both steady state and transient heating while simultaneously collecting 2D material temperature information from the Raman signal. A 3D numerical model built using the finite volume method simulates the two differentiated thermal transport processes. The simulated temperature rises are compared with experimental Raman temperature information to extract G. In this study, we investigate the temporal effects of the pulse frequency and pulse width of the ns laser on precise characterization of G. By running multiple simulations for a range of pulsed laser frequencies and pulse widths, we determine that ET-Raman is most sensitive to precise interfacial thermal conductance measurements with short pulse widths (~below 20 ns) and with pulsing frequencies that ensure sufficient cooling time between consecutive pulses (~below 500 kHz). We also conclude that substrate thickness and laser spot size play critical roles in heat distribution and subsequent sensitivity to meaningful G measurements. As both laser spot size and SiO 2 substrate thickness decrease, detection of interfacial thermal transport becomes more difficult. These results offer essential guidance in the design of future Raman-based approaches to thermal measurements using pulsed lasers.

2D materials↗
DOE OSTI2025

Quantum entropy as a harbinger of factorizability

Deeply inelastic scattering (DIS) is a powerful probe for investigating the QCD structure of hadronic matter and testing the standard model (SM). DIS can be described through QCD factorization theorems which separate contributions to the scattering interaction arising from disparate scales — e.g ., with nonperturbative matrix elements associated with long distances and a perturbative hard scattering kernel applying to short-distance parton-level interactions. The fundamental underpinnings of factorization may be recast in the quantum-theoretic terms of entanglement, (de)coherence, and system localization in a fashion which sheds complementary light on the dynamics at work in DIS from QCD bound states. In this Letter, we propose and quantitatively test such a quantum-information theoretic approach for dissecting factorization in DIS and its domain of validity; we employ metrics associated with quantum entanglement such as a differential quantum entropy and associated Kullback-Leibler (KL) divergences in numerical tests. We deploy these methods on an archetypal quark-spectator model of the proton, for which we monitor quantum decoherence in DIS as underlying model parameters are varied. On this basis, we demonstrate quantitatively how factorization-breaking effects may be imprinted on quantum entropies in a kinematic regime where leading-twist factorization increasingly receives large corrections from finite- Q 2 effects; our findings suggest potential applications of quantum simulation to QCD systems and their interactions.

Deep inelastic scattering↗
DOE OSTI2024

Learning nonlinear operators in latent spaces for real-time predictions of complex dynamics in physical systems

Abstract Predicting complex dynamics in physical applications governed by partial differential equations in real-time is nearly impossible with traditional numerical simulations due to high computational cost. Neural operators offer a solution by approximating mappings between infinite-dimensional Banach spaces, yet their performance degrades with system size and complexity. We propose an approach for learning neural operators in latent spaces, facilitating real-time predictions for highly nonlinear and multiscale systems on high-dimensional domains. Our method utilizes the deep operator network architecture on a low-dimensional latent space to efficiently approximate underlying operators. Demonstrations on material fracture, fluid flow prediction, and climate modeling highlight superior prediction accuracy and computational efficiency compared to existing methods. Notably, our approach enables approximating large-scale atmospheric flows with millions of degrees, enhancing weather and climate forecasts. Here we show that the proposed approach enables real-time predictions that can facilitate decision-making for a wide range of applications in science and engineering.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Absolute stimulated Raman side scatter in direct-drive laser-produced plasmas

Production of suprathermal electrons by stimulated Raman scattering (SRS) is a principal concern for contemporary direct-drive inertial confinement fusion experiments at the National Ignition Facility and similar systems since such electrons penetrate and preheat the target core, preventing efficient implosion. The higher temperatures and longer scale lengths in these experiments favor SRS over two-plasmon decay, which predominated in earlier experiments. In particular, current experiments are expected to exceed the threshold for absolute Raman side scatter, which would then dominate the interaction since it grows temporally until saturated by nonlinear mechanisms such as hot-electron production. Until recently, analyses of SRS side scatter have treated the case of a single laser beam incident on a plasma, but the direct-drive approach to laser fusion employs a multitude of beams to drive the implosion. Here in this paper, we present an analysis that can be applied to an arbitrary number of beams with varied angles of incidence and polarizations. In the case of a single beam, it allows a physically motivated derivation and verification of an analytic threshold formula. In the general case of multiple beams and arbitrary orientation and polarizations, the threshold is found by numerical integration of a set of first-order linear partial differential equations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Autodifferentiable Ensemble Kalman Filters

Data assimilation is concerned with sequentially estimating a temporally evolving state. This task, which arises in a wide range of scientific and engineering applications, is particularly challenging when the state is high-dimensional and the state-space dynamics are unknown. This paper introduces a machine learning framework for learning dynamical systems in data assimilation. Here, our auto-differentiable ensemble Kalman filters (AD-EnKFs) blend ensemble Kalman filters for state recovery with machine learning tools for learning the dynamics. In doing so, AD-EnKFs leverage the ability of ensemble Kalman filters to scale to high-dimensional states and the power of automatic differentiation to train high-dimensional surrogate models for the dynamics. Numerical results using the Lorenz-96 model show that AD-EnKFs outperform existing methods that use expectation-maximization or particle filters to merge data assimilation and machine learning. In addition, AD-EnKFs are easy to implement and require minimal tuning.

autodifferentiation↗
DOE OSTI2023

Association of Serum Bile Acid and Unsaturated Fatty Acid Profiles with the Risk of Diabetic Retinopathy in Type 2 Diabetic Patients

Aim: We aimed to identify the ability of serum bile acids (BAs) and unsaturated fatty acids (UFAs) profiles to predict the development of diabetic retinopathy (DR) in type 2 diabetes mellitus (T2DM) patients. Methods: We first used univariate and multivariate analysis to compare 15 serum BA and 11 UFA levels in healthy control (HC) group (n = 82), T2DM patients with DR (n = 58) and T2DM patients without DR (n = 60). Forty T2DM patients were considered for validation. Then, the receiver operating characteristic curve (ROC) and decision curve analysis were used to assess the diagnostic value and clinical benefit of serum biomarkers alone, clinical variables alone or in combination, and the area under the curve (AUC), integrated discrimination improvement (IDI), and net reclassification improvement (NRI) were used to further assess whether the addition of biomarkers significantly improved the predictive ability of the model. Results: Orthogonal partial least squares-discriminant analysis (OPLS-DA) of serum BAs and UFAs separated the three cohorts including HC, T2DM patients with or without DR. The difference in serum BA and UFA profiles of T2DM patients with or without DR was mainly manifested in the three metabolites of taurolithocholic acid (TLCA), tauroursodeoxycholic acid (TUDCA) and arachidonic acid (AA). Together, they had an AUC of 0.785 (0.918 for validation cohort) for predicting DR in T2DM patients. After adjusting for numerous confounding factors, TLCA, TUDCA, and AA were independent predictors that differentiated T2DM with or without DR. The results of AUC, IDI, and NRI demonstrated that adding these three biomarkers to a model with clinical variables statistically increased their predictive value and were replicated in our independent validation cohort. Conclusion: These findings highlight the association of three metabolites, TLCA, TUDCA and AA, with DR and may indicate their potential value in the pathogenesis of DR.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1974

Contribution to the theory of tidal oscillations of an elastic earth. External tidal potential

The differential equations of the tidal oscillations of the earth were established under the assumption that the interior of the earth is laterally inhomogeneous. The theory was developed using vectorial and dyadic symbolism to shorten the exposition and to reduce the differential equations to a symmetric form convenient for programming and for numerical integration. The formation of tidal buldges on the surfaces of discontinuity and the changes in the internal density produce small periodic variations in the exterior geopotential which are reflected in the motion of artificial satellites. The analoques of Love elastic parameters in the expansion of exterior tidal potential reflect the asymmetric and inhomogeneous structure of the interior of the earth.

Musen, P.↗
NASA NTRS1974

Geothermal reservoir simulation

The prediction of long-term geothermal reservoir performance and the environmental impact of exploiting this resource are two important problems associated with the utilization of geothermal energy for power production. Our research effort addresses these problems through numerical simulation. Computer codes based on the solution of partial-differential equations using finite-element techniques are being prepared to simulate multiphase energy transport, energy transport in fractured porous reservoirs, well bore phenomena, and subsidence.

Mercer, J. W., Jr.↗
NASA NTRS1978

A three-dimensional model of corotating streams in the solar wind. 1: Theoretical foundations

The theoretical and mathematical background pertinent to the study of steady, corotating solar wind structure in all three spatial dimensions (3-D) is discussed. The dynamical evolution of the plasma in interplanetary space (defined as the region beyond roughly 35 solar radii where the flow is supersonic) is approximately described by the nonlinear, single fluid, polytropic (magneto-) hydrodynamic equations. Efficient numerical techniques for solving this complex system of coupled, hyperbolic partial differential equations are outlined. The formulation is inviscid and nonmagnetic, but methods allow for the potential inclusion of both features with only modest modifications. One simple, highly idealized, hydrodynamic model stream is examined to illustrate the fundamental processes involved in the 3-D dynamics of stream evolution. Spatial variations in the rotational stream interaction mechanism were found to produce small nonradial flows on a global scale that lead to the transport of mass, energy, and momentum away from regions of relative compression and into regions of relative rarefaction.

Pizzo, V. J.↗
NASA NTRS1978

Calculation of electrostatic fields in periodic structures of complex shape

A universal algorithm is presented for calculating electrostatic fields in an infinite periodic structure consisting of electrodes of arbitrary shape which are located in mirror-symmetrical manner along the axis of electron-beam propagation. The method is based on the theory of R-functions, and the differential operators which are derived on the basis of the functions. Numerical results are presented and the accuracy of the results is examined.

Kravchenko, V. F.↗
NASA NTRS1979

Analytical modeling of flash-back phenomena

To understand the flame flash-back phenomena more extensively, an analytical model was formed and a numerical program was written and tested to solve the set of differential equations describing the model. Results show that under a given set of conditions flame propagates in the boundary layer on a flat plate when the free stream is at or below 1.8 m/s.

Feng, C. C.↗
NASA NTRS1981

Multigrid simulation of asymptotic curved-duct flows using a semi-implicit numerical technique

Asymptotic flows inside curved ducts of rectangular as well as polar cross section are analyzed using the Navier-Stokes equations in terms of the axial velocity and vorticity and the cross-flow stream function. Numerical solutions of the three second-order coupled elliptic partial differential equations governing this flow are obtained efficiently using the coupled alternating-direction implicit (ADI) method as well as the multigrid strongly-implicit (SI) scheme. For the flow configuration studied, the ADI method is found to be more sensitive to the time steps used than is the SI scheme. Use of the multigrid-coupled-strongly-implicit (MG-SI) scheme makes it possible to efficiently obtain fine-grid solutions for configurations having strong secondary flow. It is shown that, for this asymptotic curved-duct flow, the similarity parameter of significance is the Dean's number K rather than the Reynolds number Re. Results are obtained for curved ducts with square cross sections for K up to 900, which here corresponds to Re = 9,000 for this internal flow configuration.

Ghia, K. N.↗