SEARCH · Engineering Papers

Results for “numerical approximations”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 523 records · Page 29

All sources NASA NTRS DOE OSTI

NASA NTRS1982

On the statistical uncertainties associated with line profile fitting

The paper considers the statistical uncertainties in the best-fit parameters associated with the fitting of observational data by model line profiles of specific form. The parameters involved are amplitude, center, and width. A general procedure is described for any profile type of parametric form. The case of Gaussian profiles is treated in detail, and numerical solutions are presented for equally spaced data of equal weight. Approximate analytical solutions are derived that give good general agreement with the numerical solutions over the whole range of interest. Approximate solutions are also obtained for the case of Lorentzian profiles. It is shown that the Lorentzian results can be used in conjunction with the Gaussian results to estimate the corresponding parameter uncertainties in the case of the general Voigt profile.

Landman, D. A.↗

NASA NTRS1985

The numerical viscosity of entropy stable schemes for systems of conservation laws

Discrete approximations to hyperbolic systems of conservation laws are studied. The amount of numerical viscosity present in such schemes, is quantified and related to their entropy stability by means of comparison. To this end, conservative schemes which are also entropy conservative are constructed. These entropy conservative schemes enjoy second-order accuracy; moreover, they admit a particular interpretation within the finite-element frameworks, and hence can be formulated on various mesh configurations. It is then shown that conservative schemes are entropy stable if and only if they contain more viscosity than the mentioned above entropy conservative ones.

Tadmor, E.↗

NASA NTRS1987

The numerical viscosity of entropy stable schemes for systems of conservation laws. I

Discrete approximations to hyperbolic systems of conservation laws are studied. The amount of numerical viscosity present in such schemes is quantified and related to their entropy stability by means of comparison.To this end, conservative schemes which are also entropy-conservative are constructed. These entropy-conservative schemes enjoy second-order accuracy; moreover, they can be interpreted as piecewise-linear finite-element methods, and hence can be formulated on various mesh configurations. It is then shown that conservative schemes are entropy stable, if and (for three-point schemes) only they contain more viscosity than that present in the above-mentioned entropy-conservative ones.

Tadmor, Eitan↗

NASA NTRS1987

Approximate Simulation of Turbulence

Computed spectra resemble von Karman spectra of frequencies of interest. Numerical technique yields simulated atmospheric-turbulence spectra that closely approximates von Karman spectra within frequency ranges of interest for aircraft response. Suitable for computation in that they represent stable systems and roll off as f-2 at high frequencies outside range of interest.

Campbell, C. W.↗

NASA NTRS2001

B-spline Method in Fluid Dynamics

B-spline functions are bases for piecewise polynomials that possess attractive properties for complex flow simulations : they have compact support, provide a straightforward handling of boundary conditions and grid nonuniformities, and yield numerical schemes with high resolving power, where the order of accuracy is a mere input parameter. This paper reviews the progress made on the development and application of B-spline numerical methods to computational fluid dynamics problems. Basic B-spline approximation properties is investigated, and their relationship with conventional numerical methods is reviewed. Some fundamental developments towards efficient complex geometry spline methods are covered, such as local interpolation methods, fast solution algorithms on cartesian grid, non-conformal block-structured discretization, formulation of spline bases of higher continuity over triangulation, and treatment of pressure oscillations in Navier-Stokes equations. Application of some of these techniques to the computation of viscous incompressible flows is presented.

Botella, Olivier↗

NASA NTRS1990

Morphometric analysis of epidermal differentiation in primary roots of Zea mays

Epidermal differentiation in primary roots of Zea mays was divided into six cell types based on cellular shape and cytoplasmic appearance. These six cell types are: 1) apical protoderm, located at the tip of the root pole and characterized by periclinally flattened cells; 2) cuboidal protoderm, located approximately 230 microns from the root pole and characterized by cuboidal cells; 3) tabular epidermis, located approximately 450 microns from the root pole and characterized by anticlinally flattened cells; 4) cuboidal epidermis, located approximately 900 microns from the root pole and characterized by cuboidal cells having numerous small vacuoles; 5) vacuolate cuboidal epidermis, located approximately 1,500 microns from the root pole and characterized by cuboidal cells containing several large vacuoles; and 6) columnar epidermis, located approximately 2,200 microns from the root pole (i.e., at the beginning of the zone of elongation) and characterized by elongated cells. We also used stereology to quantify the cellular changes associated with epidermal differentiation. The quiescent center and the apical protoderm have significantly different ultrastructures. The relative volume of dictyosomes increases dramatically during the early stages of epidermal differentiation. This increase correlates inversely with the amount of coverage provided by the root cap and mucilage.

NASA Discipline Plant Biology↗

NASA NTRS1995

Aerodynamic sound of flow past an airfoil

The long term objective of this project is to develop a computational method for predicting the noise of turbulence-airfoil interactions, particularly at the trailing edge. We seek to obtain the energy-containing features of the turbulent boundary layers and the near-wake using Navier-Stokes Simulation (LES or DNS), and then to calculate the far-field acoustic characteristics by means of acoustic analogy theories, using the simulation data as acoustic source functions. Two distinct types of noise can be emitted from airfoil trailing edges. The first, a tonal or narrowband sound caused by vortex shedding, is normally associated with blunt trailing edges, high angles of attack, or laminar flow airfoils. The second source is of broadband nature arising from the aeroacoustic scattering of turbulent eddies by the trailing edge. Due to its importance to airframe noise, rotor and propeller noise, etc., trailing edge noise has been the subject of extensive theoretical (e.g. Crighton & Leppington 1971; Howe 1978) as well as experimental investigations (e.g. Brooks & Hodgson 1981; Blake & Gershfeld 1988). A number of challenges exist concerning acoustic analogy based noise computations. These include the elimination of spurious sound caused by vortices crossing permeable computational boundaries in the wake, the treatment of noncompact source regions, and the accurate description of wave reflection by the solid surface and scattering near the edge. In addition, accurate turbulence statistics in the flow field are required for the evaluation of acoustic source functions. Major efforts to date have been focused on the first two challenges. To this end, a paradigm problem of laminar vortex shedding, generated by a two dimensional, uniform stream past a NACA0012 airfoil, is used to address the relevant numerical issues. Under the low Mach number approximation, the near-field flow quantities are obtained by solving the incompressible Navier-Stokes equations numerically at chord Reynolds number of 104. The far-field noise is computed using Curle's extension to the Lighthill analogy (Curle 1955). An effective method for separating the physical noise source from spurious boundary contributions is developed. This allows an accurate evaluation of the Reynolds stress volume quadrupoles, in addition to the more readily computable surface dipoles due to the unsteady lift and drag. The effect of noncompact source distribution on the far-field sound is assessed using an efficient integration scheme for the Curle integral, with full account of retarded-time variations. The numerical results confirm in quantitative terms that the far-field sound is dominated by the surface pressure dipoles at low Mach number. The techniques developed are applicable to a wide range of flows, including jets and mixing layers, where the Reynolds stress quadrupoles play a prominent or even dominant role in the overall sound generation.

Wang, Meng↗

NASA NTRS1990

Robustness of controllers designed using Galerkin type approximations

One of the difficulties in designing controllers for infinite-dimensional systems arises from attempting to calculate a state for the system. It is shown that Galerkin type approximations can be used to design controllers which will perform as designed when implemented on the original infinite-dimensional system. No assumptions, other than those typically employed in numerical analysis, are made on the approximating scheme.

Morris, K. A.↗

DOE OSTI2024

Sparse Cholesky factorization for solving nonlinear PDEs via Gaussian processes

In recent years, there has been widespread adoption of machine learning-based approaches to automate the solving of partial differential equations (PDEs). Among these approaches, Gaussian processes (GPs) and kernel methods have garnered considerable interest due to their flexibility, robust theoretical guarantees, and close ties to traditional methods. They can transform the solving of general nonlinear PDEs into solving quadratic optimization problems with nonlinear, PDE-induced constraints. However, the complexity bottleneck lies in computing with dense kernel matrices obtained from pointwise evaluations of the covariance kernel, and its partial derivatives, a result of the PDE constraint and for which fast algorithms are scarce. The primary goal of this paper is to provide a near-linear complexity algorithm for working with such kernel matrices. We present a sparse Cholesky factorization algorithm for these matrices based on the near-sparsity of the Cholesky factor under a novel ordering of pointwise and derivative measurements. The near-sparsity is rigorously justified by directly connecting the factor to GP regression and exponential decay of basis functions in numerical homogenization. We then employ the Vecchia approximation of GPs, which is optimal in the Kullback-Leibler divergence, to compute the approximate factor. This enables us to compute ϵ-approximate inverse Cholesky factors of the kernel matrices with complexity O(N log d (N/ϵ)) in space and O(N log 2d (N/ϵ)) in time. We integrate sparse Cholesky factorizations into optimization algorithms to obtain fast solvers of the nonlinear PDE. We numerically illustrate our algorithm’s near-linear space/time complexity for a broad class of nonlinear PDEs such as the nonlinear elliptic, Burgers, and Monge-Ampère equations. In summary, we provide a fast, scalable, and accurate method for solving general PDEs with GPs and kernel methods.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Conservative numerical schemes with optimal dispersive wave relations: Part I. Derivation and analysis

An energy-conserving and an energy-and-enstrophy conserving numerical schemes are derived by approximating the Hamiltonian formulation of the inviscid shallow water flows based on the vorticity-divergence variables. These schemes also conserve the first-order moments such as mass and vorticity, as usual. The conservative properties of the schemes stem from the skew-symmetry and singularities of the Poisson brackets, which are carefully retained in the discrete approximations. Here, the schemes operate on unstructured orthogonal dual meshes, over bounded or unbounded domains, and they are also shown to possess the same optimal dispersive wave relations as those of the Z-grid scheme, which is a consequence of the use of the vorticity and divergence variables.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

A unified hp-HDG framework for Friedrichs' PDE systems

Here, this work proposes a unified hp-adaptivity framework for hybridized discontinuous Galerkin (HDG) method for a large class of partial differential equations (PDEs) of Friedrichs' type. In particular, we present unified hp-HDG formulations for abstract one-field and two-field structures and prove their well-posedness. In order to handle non-conforming interfaces we simply take advantage of HDG built-in mortar structures. With split-type mortars and the approximation space of trace, a numerical flux can be derived via Godunov approach and be naturally employed without any additional treatment. As a consequence, the proposed formulations are parameter-free. We perform several numerical experiments for time-independent and linear PDEs including elliptic, hyperbolic, and mixed-type to verify the proposed unified hp-formulations and demonstrate the effectiveness of hp-adaptation. Two adaptivity criteria are considered: one is based on a simple and fast error indicator, while the other is rigorous but more expensive using an adjoint-based error estimate. The numerical results show that these two approaches are comparable in terms of convergence rate even for problems with strong gradients, discontinuities, and singularities.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Performance of multifunctional piezoelectric energy harvesting microgyroscopes with material degradation

Most microelectromechanical systems (MEMS) are made from brittle materials that are prone to cracking. Damage from the formation of cracks can occur anytime during the manufacturing process or due to mechanical loading of MEMS devices. Cracks can increase the structures flexibility and lead to changes in performance or total failure. Understanding the potential influences cracks in MEMS is necessary for the longevity of these systems. The energy harvesting microgyroscope studied in this work is susceptible to crack propagation in its silicon substrate and lead zirconate titanate piezoelectric patch. To this end, an electromechanical model is developed to capture the impacts of cracks on the performance of a damaged energy harvesting microgyroscope. Numerical methods are used to approximate the damage associated with distributed crack networks. Then, the static and dynamic responses of the damaged microgyroscope are analyzed for two possible multifunctional setups. Here, the current study revealed that cracks can significantly affect the static and dynamic behaviors of the microgyroscope energy harvester. Combinations of cracking and piezoelectric property degradation cause a decrease in the static pull-in voltage and shift the resonance frequencies of the system. Additionally, it is shown that the systems sensitivity to changes in the base rotation and energy harvesting capabilities vary widely depending on the severity of material degradation.

36 MATERIALS SCIENCE↗

DOE OSTI2021

A classical, non-singular, bouncing universe

We present a model for a classical, non-singular bouncing cosmology without violation of the null energy condition (NEC). The field content is General Relativity plus a real scalar field with a canonical kinetic term and only renormalizable, polynomial-type self-interactions for the scalar field in the Jordan frame. The universe begins vacuum-energy dominated and is contracting at t=-∞. We consider a closed universe with a positive spatial curvature, which is responsible for the universe bouncing without any NEC violation. An Rφ 2 coupling between the Ricci scalar and the scalar field drives the scalar field from the initial false vacuum to the true vacuum during the bounce. The model is sub-Planckian throughout its evolution and every dimensionful parameter is below the effective-field-theory scale M P , so we expect no ghost-type or tachyonic instabilities. This model solves the horizon problem and extends co-moving particle geodesics to past infinity, resulting in a geodesically complete universe without singularities. We solve the Friedman equations and the scalar-field equation of motion numerically, and analytically under certain approximations.

Astronomy & Astrophysics↗

DOE OSTI2020

Parametrically driven nonlinear Dirac equation with arbitrary nonlinearity

Here, the damped and parametrically driven nonlinear Dirac equation with arbitrary nonlinearity parameter $\kappa$ is analyzed, when the external force is periodic in space and given by $\mathcal{f(x) = r}$ cos(Kx), both numerically and in a variational approximation using five collective coordinates (time dependent shape parameters of the wave function). Our variational approximation satisfies exactly the low-order moment equations. Because of competition between the spatial period of the external force λ = 2π/K , and the soliton width $l_s$, which is a function of the nonlinearity $\kappa$ as well as the initial frequency ω 0 of the solitary wave, there is a transition (at fixed ω 0 ) from trapped to unbound behavior of the soliton, which depends on the parameters $\mathcal{r}$ and K of the external force and the nonlinearity parameter $\kappa$. We previously studied this phenomena when $\kappa$ = 1 where we showed that for λ $\gg l_s$ the soliton oscillates in an effective potential, while forλ $\ll l_s$ it moves uniformly as a free particle. In this paper we focus on the $\kappa$ dependence of the transition from oscillatory to particle behavior and explicitly compare the curves of the transition regime found in the collective coordinate approximation as a function of $\mathcal{r}$ and K when $\kappa$ = 1/2,1,2 at fixed value of the frequency ω 0 . Since the solitary wave gets narrower for fixed ω 0 as a function of $\kappa$, we expect and indeed find that the regime where the solitary wave is trapped is extended as we increase $\kappa$.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

A Quantum Approach for Implementing Fixed-Point Arithmetic in Solving Ordinary Differential Equations

Differential equations (DEs) serve as fundamental tools in mathematical modeling across scientific disciplines, yet classical numerical solvers face limitations with large-scale or computationally intensive problems. This study explores a quantum-inspired approach to solving DEs, combining quantum-inspired techniques with classical methods. It focuses on fixed-point arithmetic on quantum circuits, utilizing basic quantum gates to manipulate DE solutions. We expand upon the techniques introduced by Zanger et al. [Quantum, 5, 502 (2021)] by offering a precise computation for a fixed-point signed multiplication scheme, while also presenting a quantum circuit capable of executing the fixed-point division algorithm. We demonstrate the feasibility of our approach through the simulation of a linear Ordinary Differential Equation (ODE), where initial conditions and parameters are encoded into quantum circuits using fixed-point representation. By executing sequences of quantum gates mimicking numerical integration steps, we obtain approximate solutions to the ODE with specified fixed-point precision.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2020

Development and Validation of Passive Yaw in the Open-Source WEC-Sim Code

A passive yaw implementation is developed, validated, and explored for the WEC-Sim, an open-source wave energy converter modeling tool that works within MATLAB/Simulink. The Reference Model 5 (RM5) is selected for this investigation, and a WEC-Sim model of the device is modified to allow yaw motion. A boundary element method (BEM) code was used to calculate the excitation force coefficients for a range of wave headings. An algorithm was implemented in WEC-Sim to determine the equivalent wave heading from a body’s instantaneous yaw angle and interpolate the appropriate excitation coefficients to ensure the correct time-domain excitation force. This approach is able to determine excitation force for a body undergoing large yaw displacement. For the mathematically simple case of regular wave excitation, the dynamic equation was integrated numerically and found to closely approximate the results from this implementation in WEC-Sim. A case study is presented for the same device in irregular waves. In this case, computation time is increased by 32x when this interpolation is performed at every time step. To reduce this expense, a threshold yaw displacement can be set to reduce the number of interpolations performed. A threshold of 0.01° was found to increase computation time by only 22x without significantly affecting time domain results. Similar amplitude spectra for yaw force and displacements are observed for all threshold values less than 1°, for which computation time is only increased by 2.2x.

50 EE - Wind and Water Power Program - Water (EE-4↗

DOE OSTI2025

Extended Galerkin Neural Network Approximation of Singular Variational Problems with Error Control

We present extended Galerkin neural networks, a variational framework for approximating general boundary value problems (BVPs) with error control. The main contributions of this work are (1) a rigorous theory guiding the construction of new weighted least squares variational formulations suitable for use in neural network approximation of general BVPs, and (2) an “extended” feedforward network architecture which incorporates and is even capable of learning singular solution structures, thus greatly improving approximability of singular solutions. Furthermore, numerical results are presented for several problems, including steady Stokes flow around reentrant corners and in convex corners with Moffatt eddies in order to demonstrate efficacy of the method.

a posteriori error estimate↗

DOE OSTI2025

A Quantum Approach for Implementing Fixed-Point Arithmetic in Solving Ordinary Differential Equations

Differential equations (DEs) serve as fundamental tools in mathematical modeling across scientific disciplines, yet classical numerical solvers face limitations with large-scale or computationally intensive problems. This study explores a quantum-inspired approach to solving DEs, combining quantum- inspired techniques with classical methods. It focuses on fixed- point arithmetic on quantum circuits, utilizing basic quantum gates to manipulate DE solutions. We expand upon the techniques introduced by Zanger et al. [Quantum, 5, 502 (2021)] by offering a precise computation for a fixed-point signed multiplication scheme, while also presenting a quantum circuit capable of executing the fixed-point division algorithm. We demonstrate the feasibility of our approach through the simulation of a linear Ordinary Differential Equation (ODE), where initial conditions and parameters are encoded into quantum circuits using fixed- point representation. By executing sequences of quantum gates mimicking numerical integration steps, we obtain approximate solutions to the ODE with specified fixed-point precision.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗