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 901 records · Page 50

Parameter estimation problems for distributed systems using a multigrid method

The problem of estimating spatially varying coefficients of partial differential equations is considered from observation of the solution and of the right hand side of the equation. It is assumed that the observations are distributed in the domain and that enough observations are given. A method of discretization and an efficient multigrid method for solving the resulting discrete systems are described. Numerical results are presented for estimation of coefficients in an elliptic and a parabolic partial differential equation.

Ta'asan, Shlomo↗

Flat optics for image differentiation

Image processing has become a critical technology in a variety of science and engineering disciplines. Although most image processing is performed digitally, optical analog processing has the advantages of being low-power and high-speed, but it requires a large volume. Here, we demonstrate flat optics for direct image differentiation, allowing us to significantly shrink the required optical system size. We first demonstrate how the differentiator can be combined with traditional imaging systems such as a commercial optical microscope and camera sensor for edge detection with a numerical aperture up to 0.32. We next demonstrate how the entire processing system can be realized as a monolithic compound flat optic by integrating the differentiator with a metalens. The compound nanophotonic system manifests the advantage of thin form factor as well as the ability to implement complex transfer functions, and could open new opportunities in applications such as biological imaging and computer vision.

36 MATERIALS SCIENCE↗

FASOR - A second generation shell of revolution code

An integrated computer program entitled Field Analysis of Shells of Revolution (FASOR) currently under development for NASA is described. When completed, this code will treat prebuckling, buckling, initial postbuckling and vibrations under axisymmetric static loads as well as linear response and bifurcation under asymmetric static loads. Although these modes of response are treated by existing programs, FASOR extends the class of problems treated to include general anisotropy and transverse shear deformations of stiffened laminated shells. At the same time, a primary goal is to develop a program which is free of the usual problems of modeling, numerical convergence and ill-conditioning, laborious problem setup, limitations on problem size and interpretation of output. The field method is briefly described, the shell differential equations are cast in a suitable form for solution by this method and essential aspects of the input format are presented. Numerical results are given for both unstiffened and stiffened anisotropic cylindrical shells and compared with previously published analytical solutions.

Cohen, G. A.↗

Stochastic journeys of cell progenies through compartments and the role of self-renewal, symmetric and asymmetric division

Abstract Division and differentiation events by which cell populations with specific functions are generated often take place as part of a developmental programme, which can be represented by a sequence of compartments. A compartment is the set of cells with common characteristics; sharing, for instance, a spatial location or a phenotype. Differentiation events are transitions from one compartment to the next. Cells may also die or divide. We consider three different types of division events: (i) where both daughter cells inherit the mother’s phenotype (self-renewal), (ii) where only one of the daughters changes phenotype (asymmetric division), and (iii) where both daughters change phenotype (symmetric division). The self-renewal probability in each compartment determines whether the progeny of a single cell, moving through the sequence of compartments, is finite or grows without bound. We analyse the progeny stochastic dynamics with probability generating functions. In the case of self-renewal, by following one of the daughters after any division event, we may construct lifelines containing only one cell at any time. We analyse the number of divisions along such lines, and the compartment where lines terminate with a death event. Analysis and numerical simulations are applied to a five-compartment model of the gradual differentiation of hematopoietic stem cells and to a model of thymocyte development: from pre-double positive to single positive (SP) cells with a bifurcation to either SP4 or SP8 in the last compartment of the sequence.

Science & Technology - Other Topics↗

Spurious Solutions Of Nonlinear Differential Equations

Report utilizes nonlinear-dynamics approach to investigate possible sources of errors and slow convergence and non-convergence of steady-state numerical solutions when using time-dependent approach for problems containing nonlinear source terms. Emphasizes implications for development of algorithms in CFD and computational sciences in general. Main fundamental conclusion of study is that qualitative features of nonlinear differential equations cannot be adequately represented by finite-difference method and vice versa.

Yee, H. C.↗

Accurate data-driven surrogates of dynamical systems for forward propagation of uncertainty

Stochastic collocation (SC) is a well-known non-intrusive method of constructing surrogate models for uncertainty quantification. In dynamical systems, SC is especially suited for full-field uncertainty propagation that characterizes the distributions of the high-dimensional solution fields of a model with stochastic input parameters. However, due to the highly nonlinear nature of the parameter-to-solution map in even the simplest dynamical systems, the constructed SC surrogates are often inaccurate. Here, this work presents an alternative approach, where we apply the SC approximation over the dynamics of the model, rather than the solution. By combining the data-driven sparse identification of nonlinear dynamics framework with SC, we construct dynamics surrogates and integrate them through time to construct the surrogate solutions. We demonstrate that the SC-over-dynamics framework leads to smaller errors, both in terms of the approximated system trajectories as well as the model state distributions, when compared against full-field SC applied to the solutions directly. We present numerical evidence of this improvement using three test problems: a chaotic ordinary differential equation, and two partial differential equations from solid mechanics.

42 ENGINEERING↗

A stochastic solver based on the residence time algorithm for crystal plasticity models

Abstract The deformation of crystalline materials by dislocation motion takes place in discrete amounts determined by the Burgers vector. Dislocations may move individually or in bundles, potentially giving rise to intermittent slip. This confers plastic deformation with a certain degree of variability that can be interpreted as being caused by stochastic fluctuations in dislocation behavior. However, crystal plasticity (CP) models are almost always formulated in a continuum sense, assuming that fluctuations average out over large material volumes and/or cancel out due to multi-slip contributions. Nevertheless, plastic fluctuations are known to be important in confined volumes at or below the micron scale, at high temperatures, and under low strain rate/stress deformation conditions. Here, we develop a stochastic solver for CP models based on the residence-time algorithm that naturally captures plastic fluctuations by sampling among the set of active slip systems in the crystal. The method solves the evolution equations of explicit CP formulations, which are recast as stochastic ordinary differential equations and integrated discretely in time. The stochastic CP model is numerically stable by design and naturally breaks the symmetry of plastic slip by sampling among the active plastic shear rates with the correct probability. This can lead to phenomena such as intermittent slip or plastic localization without adding external symmetry-breaking operations to the model. The method is applied to body-centered cubic tungsten single crystals under a variety of temperatures, loading orientations, and imposed strain rates.

42 ENGINEERING↗

Symmetry-breaking-induced rare fluctuations in a time-delay dynamic system

Inspired by the experimental and numerical findings, we study the dynamic instabilities of two coupled nonlinear delay differential equations that are used to describe the coherent oscillations between the top and bottom boundary layers in turbulent Rayleigh–Bénard convection. By introducing two sensitivity parameters for the instabilities of the top and bottom boundary layers, we find three different types of solutions, namely in-phase single-period oscillations, multi-period oscillations and chaos. The chaos solution contains rare but large amplitude fluctuations. The statistical properties of these fluctuations are consistent with those observed in the experiment for the massive eruption of thermal plumes, which causes random reversals of the large-scale circulation in turbulent Rayleigh–Bénard convection. Finally, our study thus provides new insights into the origin of rare massive eruptions and sudden changes of large-scale flow pattern that are often observed in convection systems of geophysical and astrophysical scales.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Deep learning closure models for large-eddy simulation of flows around bluff bodies

Near-wall flow simulation remains a central challenge in aerodynamics modelling: Reynolds-averaged Navier–Stokes predictions of separated flows are often inaccurate, and large-eddy simulation (LES) can require prohibitively small near-wall mesh sizes. A deep learning (DL) closure model for LES is developed by introducing untrained neural networks into the governing equations and training in situ for incompressible flows around rectangular prisms at moderate Reynolds numbers. The DL-LES models are trained using adjoint partial differential equation (PDE) optimization methods to match, as closely as possible, direct numerical simulation (DNS) data. They are then evaluated out-of-sample – for aspect ratios, Reynolds numbers and bluff-body geometries not included in the training data – and compared with standard LES models. The DL-LES models outperform these models and are able to achieve accurate LES predictions on a relatively coarse mesh (downsampled from the DNS mesh by factors of four or eight in each Cartesian direction). We study the accuracy of the DL-LES model for predicting the drag coefficient, near-wall and far-field mean flow, and resolved Reynolds stress. A crucial challenge is that the LES quantities of interest are the steady-state flow statistics; for example, a time-averaged velocity component $\langle {u}_i\rangle (x) = \lim _{t \rightarrow \infty } ({1}/{t}) \int _0^t u_i(s,x)\, {\rm d}s$ . Calculating the steady-state flow statistics therefore requires simulating the DL-LES equations over a large number of flow times through the domain. It is a non-trivial question whether an unsteady PDE model with a functional form defined by a deep neural network can remain stable and accurate on $t \in [0, \infty )$ , especially when trained over comparatively short time intervals. Our results demonstrate that the DL-LES models are accurate and stable over long time horizons, which enables the estimation of the steady-state mean velocity, fluctuations and drag coefficient of turbulent flows around bluff bodies relevant to aerodynamics applications.

Mechanics↗

Structural basis of Qng1-mediated salvage of the micronutrient queuine from queuosine-5′-monophosphate as the biological substrate

Abstract Eukaryotic life benefits from—and ofttimes critically relies upon—the de novo biosynthesis and supply of vitamins and micronutrients from bacteria. The micronutrient queuosine (Q), derived from diet and/or the gut microbiome, is used as a source of the nucleobase queuine, which once incorporated into the anticodon of tRNA contributes to translational efficiency and accuracy. Here, we report high-resolution, substrate-bound crystal structures of the Sphaerobacter thermophilus queuine salvage protein Qng1 (formerly DUF2419) and of its human ortholog QNG1 (C9orf64), which together with biochemical and genetic evidence demonstrate its function as the hydrolase releasing queuine from queuosine-5′-monophosphate as the biological substrate. We also show that QNG1 is highly expressed in the liver, with implications for Q salvage and recycling. The essential role of this family of hydrolases in supplying queuine in eukaryotes places it at the nexus of numerous (patho)physiological processes associated with queuine deficiency, including altered metabolism, proliferation, differentiation and cancer progression.

59 BASIC BIOLOGICAL SCIENCES↗

Thermoelastic damping and its effect on flutter of stressed panels situated in a supersonic airflow

The effects of material damping on flutter of stressed rectangular panels are studied within the context of linear thermoelasticity theory. The closed-form expression for the thermoelastic (material) damping coefficient is obtained as a function of frequency, panel temperature and dimensions, and material properties. The solution of the stability boundary value problem is obtained by use of a generalized Galerkin method in the cross stream direction which reduces the governing partial differential equations to a system of ordinary differential equations in the streamwise direction. These equations are then solved exactly. Numerical results are given for the thermoelastic damping coefficients and for the flutter speeds of partially and fully clamped panels subjected to midplane stress.

Shieh, R. C.↗

Numerical solutions of atmospheric flow over semielliptical simulated hills

Atmospheric motion over obstacles on plane surfaces to compute simulated wind fields over terrain features was studied. Semielliptical, two dimensional geometry and numerical simulation of flow over rectangular geometries is also discussed. The partial differential equations for the vorticity, stream function, turbulence kinetic energy, and turbulence length scale were solved by a finite difference technique. The mechanism of flow separation induced by a semiellipse is the same as flow over a gradually sloping surface for which the flow separation is caused by the interaction between the viscous force, the pressure force, and the turbulence level. For flow over bluff bodies, a downstream recirculation bubble is created which increases the aspect ratio and/or the turbulence level results in flow reattachment close behind the obstacle.

Shieh, C. F.↗

Estimation of atmospheric species concentrations from remote sensing data

The atmospheric species concentration distribution-estimation problem that arises from the measurement configuration of remote sensing, which involves a field of view that moves across the region, is solved by means of the partial differential equation for atmospheric diffusion and the Wiener-Hopf theory. A numerical study is also undertaken for the estimation of concentration distribution downwind of a hypothetical, continuous ground-level source of pollutants.

Omatu, S.↗

Parameter identification for an abstract Cauchy problem by quasilinearization

A parameter identification problem is considered in the context of a linear abstract Cauchy problem with a parameter-dependent evolution operator. Conditions are investigated under which the gradient of the state with respect to a parameter possesses smoothness properties which lead to local convergence of an estimation algorithm based on quasi-linearization. Numerical results are presented concerning estimation of unknown parameters in delay-differential equations.

Brewer, Dennis W.↗

Origin and thermal evolution of Mars

The thermal evolution of Mars is governed by subsolidus mantle convection beneath a thick lithosphere. Models of the interior evolution are developed by parameterizing mantle convective heat transport in terms of mantle viscosity, the superadiabatic temperature rise across the mantle, and mantle heat production. Geological, geophysical, and geochemical observations of the compositon and structure of the interior and of the timing of major events in Martian evolution are used to constrain the model computations. Such evolutionary events include global differentiation, atmospheric outgassing, and the formation of the hemispherical dichotomy and Tharsis. Numerical calculations of fully three-dimensional, spherical convection in a shell the size of the Martian mantle are performed to explore plausible patterns of Martian mantel convection and to relate convective features, such as plumes, to surface features, such as Tharsis. The results from the model calculations are presented.

Schubert, Gerald↗