Search NASA⌕ Search

SEARCH · Search NASA

Results for “Closed form solution”

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 685 records · Page 38

An Engineering Solution for Using Coarse Meshes in the Simulation of Delamination with Cohesive Zone Models

This paper presents a methodology to determine the parameters used in the simulation of delamination in composite materials using decohesion finite elements. A closed-form expression is developed to define the stiffness of the cohesive layer. A novel procedure that allows the use of coarser meshes of decohesion elements in large-scale computations is proposed. The procedure ensures that the energy dissipated by the fracture process is correctly computed. It is shown that coarse-meshed models defined using the approach proposed here yield the same results as the models with finer meshes normally used in the simulation of fracture processes.

Turon, Albert↗

Chemical Design of Pb-Free Relaxors for Giant Capacitive Energy Storage

Dielectric capacitors have captured substantial attention for advanced electrical and electronic systems. Developing dielectrics with high energy density and high storage efficiency is challenging owing to the high compositional diversity and the lack of general guidelines. Herein, we propose a map that captures the structural distortion ($\delta$) and tolerance factor (t) of perovskites to design Pb-free relaxors with extremely high capacitive energy storage. Our map shows how to select ferroelectric with large $\delta$ and paraelectric components to form relaxors with a t value close to 1 and thus obtaining eliminated hysteresis and large polarization under a high electric breakdown. Taking the Bi 0.5 Na 0.5 TiO 3 -based solid solution as an example, we demonstrate that composition-driven predominant order-disorder characteristic of local atomic polar displacements endows the relaxor with a slushlike structure and strong local polar fluctuations at several nanoscale. This leads to a giant recoverable energy density of 13.6 J cm -3 , along with an ultrahigh efficiency of 94%, which is far beyond the current performance boundary reported in Pb-free bulk ceramics. In conclusion, our work provides a solution through rational chemical design for obtaining Pb-free relaxors with outstanding energy-storage properties.

25 ENERGY STORAGE↗

AI‐Accelerated Optimization of Self‐Assembled Organic Mixed Ionic‐Electronic Conductors (OMIEC) (Final Report)

This document describes research activities, products and outcomes of a DOE-funded program to help accelerate development of Organic Mixed Ionic Electronic Conductors (OMIECs) using neutron scattering and high-throughput experimentation. OMIECs are organic materials that conduct both ionic and electronic charge carriers. For these new OMIEC materials, self-assembling ion-conducting block copolymers (BCPs) are used as a structural template for electronic conducting polymers. This forms OMIECs with long-range structural order that can help facilitate long-range electronic transport. The conductive properties of OMIECs are closely associated with their structure, which is affected by solution conditions and polymer macromolecular designs. Thus high-throughput experimentation has been implemented to explore this large design space effectively. The BCP-CP OMIEC systems explored are composed of ion-conducting diblock or triblock copolymers containing polyethylene oxide (PEO), di(ethylene glycol) ethyl ether acrylate (DEGEEA) or poly(ethylene glycol) methyl ether acrylate) (PEGMEA) hydrophilic blocks. These blocks will be coupled with either polypropylene oxide (PPO) or polyheptafluorobutyl acrylate (PHFBA) hydrophobic blocks to drive assembly. The electronic conducting component consists of several different types of conjugated polymers. Small angle scattering of neutrons and X-rays (SANS/SAXS) as well as electrochemical analysis have been coupled with modern algorithms for autonomous research using artificial intelligence (AI).

36 MATERIALS SCIENCE↗

A viable suboptimal controller for robotic manipulators

A simple algebraic solution is obtained to an optimal control problem in joint space for a general robotic manipulator. Open loop and closed loop control laws are obtained. The resulting control algorithm is similar in form to a commonly used heuristic controller using cubic polynomial trajectories.

Turner, T. L.↗

Unsteady supersonic aerodynamic theory for interfering surfaces by the method of potential gradient

A generalized solution of the hyperbolic wave equation was further developed to relate the velocity components at a field point to the potential gradient distribution in the dependence domain. Singular integrals were evaluated in closed form, with numerical integration methods for more complex but analytic functions. Idealization of the lifting surfaces by trapezoidal elements with two sides parallel to the streamlines is computationally efficient. Streamwise integrals were performed analytically, and spanwise integrals were neccessary only on element leading and trailing sides. All integrands vanish on the Mach cone. Pressure distribution on a double delta wing and generalized aerodynamic coefficients for three AGARD planforms were calculated and compared with available results.

Jones, W. P.↗

Analytic solutions for single and multiple cylinders of gravitating polytropes in magnetostatic equilibrium

Exact analytic solutions for the static equilibrium of a gravitating plasma polytrope in the presence of magnetic fields are presented. The means of generating various equilibrium configurations to illustrate directly the complex physical relationships between pressure, magnetic fields, and gravity in self-gravitating systems is demonstrated. One of the solutions is used to model interstellar clouds suspended by magnetic fields against the galactic gravity such as may be formed by the Parker (1966) instability. It is concluded that the pinching effect of closed loops of magnetic fields in the clouds may be a dominant agent in further collapsing the clouds following their formation.

Lerche, I.↗

An efficient computational method for characterizing the effects of random surface errors on the average power pattern of reflectors

Based on the works of Ruze (1966) and Vu (1969), a novel mathematical model has been developed to determine efficiently the average power pattern degradations caused by random surface errors. In this model, both nonuniform root mean square (rms) surface errors and nonuniform illumination functions are employed. In addition, the model incorporates the dependence on F/D in the construction of the solution. The mathematical foundation of the model rests on the assumption that in each prescribed annular region of the antenna, the geometrical rms surface value is known. It is shown that closed-form expressions can then be derived, which result in a very efficient computational method for the average power pattern. Detailed parametric studies are performed with these expressions to determine the effects of different random errors and illumination tapers on parameters such as gain loss and sidelobe levels. The results clearly demonstrate that as sidelobe levels decrease, their dependence on the surface rms/wavelength becomes much stronger and, for a specified tolerance level, a considerably smaller rms/wavelength is required to maintain the low sidelobes within the required bounds.

Rahmat-Samii, Y.↗

Transonic viscous flow calculations for a turbine cascade with a two equation turbulence model

A numerical method for the study of steady, transonic, turbulent viscous flow through plane turbine cascades is presented. The governing equations are written in Favre-averaged form and closed with a first order model. The turbulent quantities are expressed according to a two-equation kappa-epsilon model where low Reynolds number and compressibility effects are included. The solution is obtained by using a pseudo-unsteady method with improved perturbation propagation properties. The equations are discretized in space by using a finite volume formulation. An explicit multistage dissipative Runge-Kutta algorithm is then used to advance the flow equations in the pseudo-time. First results of calculations compare fairly well with experimental data.

Boretti, A. A.↗

Vortex/surface interaction

This paper considers the interaction of a vortex generated upstream in a flow field with a downstream aerodynamic surface that possesses a large chord. The flow is assumed to be steady, incompressible, inviscid and irrotational, and the surface to be semiinfinite. The vortex is considered to be a straight vortex filament. To lowest order the problem is modeled using potential theory, where the 3D Laplace's equation for the velocity potential on the surface is solved exactly. The closed-form equation for pressure distribution obtained from this theory is found to have a square root singularity at the leading-edge. It also converges, as x goes to infinity, to the solution of the 2D point-vortex/infinite plane problem. The pressure coefficient presents an anti-symmetric behavior, near the leading-edge and a symmetric behavior as x goes to infinity.

Bodstein, G. C. R.↗

Ultrascale Visualization of Climate Data

Fueled by exponential increases in the computational and storage capabilities of high-performance computing platforms, climate simulations are evolving toward higher numerical fidelity, complexity, volume, and dimensionality. These technological breakthroughs are coming at a time of exponential growth in climate data, with estimates of hundreds of exabytes by 2020. To meet the challenges and exploit the opportunities that such explosive growth affords, a consortium of four national laboratories, two universities, a government agency, and two private companies formed to explore the next wave in climate science. Working in close collaboration with domain experts, the Ultrascale Visualization Climate Data Analysis Tools (UV-CDAT) project aims to provide high-level solutions to a variety of climate data analysis and visualization problems.

analysis tools↗

Combustion response modeling for composite solid propellants

A computerized mathematical model of the combustion response function of composite solid propellants was developed with particular attention to the contributions of the solid phase heterogeneity. The one-dimensional model treats the solid phase as alternating layers of ammonium perchlorate and binder, with an exothermic melt layer at the surface. Solution of the Fourier heat equation in the solid provides temperature and heat flux distributions with space and time. The problem is solved by conserving the heat flux at the surface from that produced by a suitable model of the gas phase. An approximation of the BDP flame model is utilized to represent the gas phase. By the use of several reasonable assumptions, it is found that a significant portion of the problem can be solved in closed form. A method is presented by which the model can be applied to tetramodal particle size distributions. A computerized steady-state version of the model was completed, which served to validate the various approximations and lay a foundation for the combustion response modeling. The combustion response modeling was completed in a form which does not require an iterative solution, and some preliminary results were acquired.

Source record↗

A disturbance based control/structure design algorithm

Some authors take a classical approach to the simultaneous structure/control optimization by attempting to simultaneously minimize the weighted sum of the total mass and a quadratic form, subject to all of the structural and control constraints. Here, the optimization will be based on the dynamic response of a structure to an external unknown stochastic disturbance environment. Such a response to excitation approach is common to both the structural and control design phases, and hence represents a more natural control/structure optimization strategy than relying on artificial and vague control penalties. The design objective is to find the structure and controller of minimum mass such that all the prescribed constraints are satisfied. Two alternative solution algorithms are presented which have been applied to this problem. Each algorithm handles the optimization strategy and the imposition of the nonlinear constraints in a different manner. Two controller methodologies, and their effect on the solution algorithm, will be considered. These are full state feedback and direct output feedback, although the problem formulation is not restricted solely to these forms of controller. In fact, although full state feedback is a popular choice among researchers in this field (for reasons that will become apparent), its practical application is severely limited. The controller/structure interaction is inserted by the imposition of appropriate closed-loop constraints, such as closed-loop output response and control effort constraints. Numerical results will be obtained for a representative flexible structure model to illustrate the effectiveness of the solution algorithms.

Mclaren, Mark D.↗

Carrier tracking loop performance in the presence of strong CW interference

The analyses of Bruno (1973) and Blanchard (1974, 1975) of the effects of strong radio-frequency interference on the behavior of a phase-locked loop are examined in detail. The lock constraints derived directly from a steady-state trial solution by Bruno and those derived with simplifying assumptions from the same trial solution by Blanchard are compared, and circumstances which may be accommodated by the results of Bruno but not those of Blanchard are pointed out. An explicit approximation is then derived which extends the closed-form approximation of Bruno for the loop phase error to the entire lock region. Results are illustrated for a loop filter representative of a phase-locked loop in a deep space network receiver.

Levitt, B. K.↗

Remote sensing of earth terrain

In remote sensing, the encountered geophysical media such as agricultural canopy, forest, snow, or ice are inhomogeneous and contain scatters in a random manner. Furthermore, weather conditions such as fog, mist, or snow cover can intervene the electromagnetic observation of the remotely sensed media. In the modelling of such media accounting for the weather effects, a multi-layer random medium model has been developed. The scattering effects of the random media are described by three-dimensional correlation functions with variances and correlation lengths corresponding to the fluctuation strengths and the physical geometry of the inhomogeneities, respectively. With proper consideration of the dyadic Green's function and its singularities, the strong fluctuation theory is used to calculate the effective permittivities which account for the modification of the wave speed and attenuation in the presence of the scatters. The distorted Born approximation is then applied to obtain the correlations of the scattered fields. From the correlation of the scattered field, calculated is the complete set of scattering coefficients for polarimetric radar observation or brightness temperature in passive radiometer applications. In the remote sensing of terrestrial ecosystems, the development of microwave remote sensing technology and the potential of SAR to measure vegetation structure and biomass have increased effort to conduct experimental and theoretical researches on the interactions between microwave and vegetation canopies. The overall objective is to develop inversion algorithms to retrieve biophysical parameters from radar data. In this perspective, theoretical models and experimental data are methodically interconnected in the following manner: Due to the complexity of the interactions involved, all theoretical models have limited domains of validity; the proposed solution is to use theoretical models, which is validated by experiments, to establish the region in which the radar response is most sensitive to the parameters of interest; theoretically simulated data will be used to generate simple invertible models over the region. For applications to the remote sensing of sea ice, the developed theoretical models need to be tested with experimental measurements. With measured ground truth such as ice thickness, temperature, salinity, and structure, input parameters to the theoretical models can be obtained to calculate the polarimetric scattering coefficients for radars or brightness temperature for radiometers and then compare theoretical results with experimental data. Validated models will play an important role in the interpretation and classification of ice in monitoring global ice cover from space borne remote sensors in the future. We present an inversion algorithm based on a recently developed inversion method referred to as the Renormalized Source-Type Integral Equation approach. The objective of this method is to overcome some of the limitations and difficulties of the iterative Born technique. It recasts the inversion, which is nonlinear in nature, in terms of the solution of a set of linear equations; however, the final inversion equation is still nonlinear. The derived inversion equation is an exact equation which sums up the iterative Neuman (or Born) series in a closed form and, thus, is a valid representation even in the case when the Born series diverges; hence, the name Renormalized Source-Type Integral Equation Approach.

Yueh, Herng-Aung↗

Formulating Precursors for Coating Metals and Ceramics

A protocol has been devised for formulating low-vapor-pressure precursors for protective and conversion coatings on metallic and ceramic substrates. The ingredients of a precursor to which the protocol applies include additives with phosphate esters, or aryl phosphate esters in solution. Additives can include iron, chromium, and/or other transition metals. Alternative or additional additives can include magnesium compounds to facilitate growth of films on substrates that do not contain magnesium. Formulation of a precursor begins with mixing of the ingredients into a high-vapor-pressure solvent to form a homogeneous solution. Then the solvent is extracted from the solution by evaporation - aided, if necessary, by vacuum and/or slight heating. The solvent is deemed to be completely extracted when the viscosity of the remaining solution closely resembles the viscosity of the phosphate ester or aryl phosphate ester. In addition, satisfactory removal of the solvent can be verified by means of a differential scanning calorimetry essay: the absence of endothermic processes for temperatures below 150 C would indicate that the residual solvent has been eliminated from the solution beyond a detectable dilution level.

Morales, Wilfredo↗

Electron Transfer in Microemulsion-Based Electrolytes

The use of flowing electrochemical reactors, for example, in redox flow batteries and in various electrosynthesis processes, is increasing. This technology has the potential to be of central significance in the increased deployment of renewable electricity for carbon-neutral processes. A key element of optimizing efficiency of electrochemical reactors is the combination of high solution conductivity and reagent solubility. In this work, we show a substantial rate of charge transfer for an electrochemical reaction occurring in a microemulsion containing electroactive material is loaded inside the nonpolar (toluene) subphase of the microemulsion. The measured rate constant translates to an exchange current density comparable to that in redox flow batteries. The rate could be controlled by the surfactant, which maintains partitioning of reactants and products by forming an interfacial region with ions in the aqueous phase in close proximity. The hypothesized echanism is evocative of membrane-bound enzymatic reactions. Achieving sufficient rates of electrochemical reaction is the product of an effort designed to establish a reaction condition that meets the requirements of electrochemical reactors using microemulsions to realize a separation of conducting and reactive elements of the solution, opening a door to the broad use of microemulsions to effect controlled electrochemical reactions as steps in more complex processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Lattice Green’s Functions for High-Order Finite Difference Stencils

Lattice Green's Functions (LGFs) are fundamental solutions to discretized linear operators, and as such they are a useful tool for solving discretized elliptic PDEs on domains that are unbounded in one or more directions. The majority of existing numerical solvers that make use of LGFs rely on a second-order discretization and operate on domains with free-space boundary conditions in all directions. Under these conditions, fast expansion methods are available that enable precomputation of 2D or 3D LGFs in linear time, avoiding the need for brute-force multi-dimensional quadrature of numerically unstable integrals. Here we focus on higher-order discretizations of the Laplace operator on domains with more general boundary conditions, by (1) providing an algorithm for fast and accurate evaluation of the LGFs associated with high-order dimension-split centered finite differences on unbounded domains, and (2) deriving closed-form expressions for the LGFs associated with both dimension-split and Mehrstellen discretizations on domains with one unbounded dimension. Through numerical experiments we demonstrate that these techniques provide LGF evaluations with near machine-precision accuracy, and that the resulting LGFs allow for numerically consistent solutions to high-order discretizations of the Poisson's equation on fully or partially unbounded 3D domains.

97 MATHEMATICS AND COMPUTING↗

Exact solution of the relativistic Dunkl oscillator in (2+1) dimensions

In this paper we study the (2+1)-dimensional Dirac–Dunkl oscillator coupled to an external magnetic field. Our Hamiltonian is obtained from the standard Dirac oscillator coupled to an external magnetic field by changing the partial derivatives by the Dunkl derivatives. We solve the Dunkl–Klein–Gordon-type equations in polar coordinates in a closed form. The angular part eigenfunctions are given in terms of the Jacobi–Dunkl polynomials and the radial functions in terms of the Laguerre functions. Also, we compute the energy spectrum of this problem and show that, in the non-relativistic limit, it properly reduces to the Hamiltonian of the two dimensional harmonic oscillator.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗