SEARCH · Engineering Papers
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.
We consider one-dimensional Carleman and Godunov–Sultangazin equations and obtain local equilibrium conditions for solutions of the Cauchy problem with finite energy and periodic initial data. Moreover, we prove the exponential stabilization to the equilibrium state.
In this paper, we describe properties of the Dedekind η-function, constructions arising from it, and their applications in various topics of the number theory and algebra. We discuss connections with the group representation theory and the study of the structure of spaces of modular forms. Special attention is paid to the special class of modular forms, namely, η-functions with multiplicative coefficients.
We formulate a complete system of relations of the local gradient electromagnetothermomechanics of electrically conductive nonferromagnetic polarizable solid media. The nonlocal character of the constitutive relations of the proposed mathematical model is explained by the presence of electric quadrupole moments in the polarization current. As a result of taking into account these moments, the space of parameters of the thermodynamic state of the body is expanded by including a pair of additional conjugate parameters, namely, the quadrupole moment and the gradient of the vector of electric-field intensity. It is shown that the developed model takes into account the electromechanical interaction for materials with high level of symmetry (isotropic materials) and describes the flexoelectric and thermopolarization effects. We also present the key system of equations for a physically and geometrically linear medium.
We present sufficient conditions for the existence of a weak solution of the Dirichlet problem for the heat-conduction equation with random action described by an integral over the general stochastic measure.
In order to express a polyhedron as the Minkowski sum of a polytope and a polyhedral cone, Motzkin (Beiträge zur Theorie der linearen Ungleichungen. Dissertation, University of Basel, 1936) devised a homogenization technique that translates the polyhedron to a polyhedral cone in one higher dimension. Refining his technique, we present a conical representation of a set in the Euclidean space. Then, we use this representation to reach four main results: First, we establish a convex programming based framework for determining a maximal element—an element with the maximum number of positive components—of a non-negative convex set—a convex set in the non-negative Euclidean orthant. Second, we develop a linear programming problem for finding a relative interior point of a polyhedron. Third, we propose two procedures for identifying a strictly complementary solution in linear programming. Finally, we generalize Motzkin’s (Beiträge zur Theorie der linearen Ungleichungen. Dissertation, University of Basel, 1936) representation theorem for a class of closed convex sets in the Euclidean space.
We consider the discretization of the stationary Navier–Stokes/Darcy system in a two-dimensional domain by the non-conforming finite volume element method. We use the standard formulation of the Navier–Stokes/Darcy system in the primitive variables and take as approximation space the non-conforming P{sub 1} elements for velocity and piezometric head and piecewise constant elements for the hydrostatic pressure. We prove that the unique solution of the non-conforming finite volume element method converges to the true solution with optimal order for velocity and piezometric head in discrete H{sup 1} norm and for pressure in discrete L{sup 2} norm, respectively. Finally, some numerical experiments are presented to validate our theoretical results.
The paper is motivated by a model for the injection of air into a porous medium that contains a solid fuel. In previous work, a system of three evolutionary partial differential equations that models combustion of light porous foam under air injection was considered. The existence and uniqueness of traveling waves were studied and the wave sequences appearing in Riemann solutions were identified. This analysis was done under assumption that the combustion wave velocity is positive. In the present work, this hypothesis was neglected and the results generalized including the case of negative combustion wave speed. The existence of such waves was proved and the uniqueness investigated for some particular cases using Melnikov’s integral. The wave sequences appearing in Riemann solutions were identified and the numerical examples using finite difference scheme were presented.
We consider the use of surfaces as topological spaces in order to provide an answer to the quantization problem of a discrete memoryless channel (DMC). The aim is to identify the geometric and algebraic properties of the surface associated with the embedding of such a DMC. The procedure is based on the following steps: knowing the complete bipartite graph associated with a given DMC: (1) to determine the minimum and the maximum genus of the corresponding surfaces in which the given graph is embedded; (2) knowing the genus, to specify the roots as the elements of a Farey sequence of a planar algebraic curve; (3) to determine the solutions of a second-order Fuchsian differential equation as the generators of the corresponding Fuchsian group. We consider the cases of the hypergeometric and Heun equations, with three and four regular singular points, respectively. By means of this procedure, the fundamental region associated with the Fuchsian group is identified and this is the region where the planar algebraic curve is uniformized. A generalization of the quantization problem related to an m-ary input n-ary output symmetric channel viewed as a complete bipartite graph K{sub m,n}, by use of the embedding of this channel in compact surfaces is straightforward.
The goal of this paper is to develop a highly accurate and efficient numerical method for the solution of a time-dependent partial differential equation with a piecewise constant coefficient, on a finite interval with periodic boundary conditions. The resulting algorithm can be used, for example, to model the diffusion of heat energy in one space dimension, in the case where the spatial domain represents a medium consisting of two homogeneous materials. The resulting model has, to our knowledge, not yet been solved in closed form through analytical methods, and is difficult to solve using existing numerical methods, thus suggesting an alternative approach. The approach presented in this paper is to represent the solution as a linear combination of wave functions that change frequencies at the interfaces between different materials. It is demonstrated through numerical experiments that using the Uncertainty Principle to construct a basis of such functions, in conjunction with a spectral method, a mathematical model for heat diffusion through different materials can be solved much more efficiently than with conventional time-stepping methods.
In this paper, we investigate a mathematical model for age-structured forest–beetle interactions that includes harvesting of trees. The aim is to broaden the understanding of the synergistic effects of harvesting and insect infestation on the age structure of the forests and the harvesting benefit. In the first part of this study, we consider different scenarios of the forest infestation by beetles and observe that the quantitative age profile of the forest depends significantly on whether the beetle population is in its endemic or epidemic states. In the second part, we also include harvesting of the forest trees and analyze two different harvesting strategies: cutting all trees older than a certain age, and cutting a fixed proportion of trees older than a certain age. Numerical simulations are implemented to determine the optimal cutting age for both harvesting strategies. The numerical simulations reveal that, independent of the steady state of the beetle population (that is, no beetles, endemic or epidemic state) clear cutting all trees older than a given age provides a higher harvesting benefit. Our numerical simulations further indicate that to obtain a fixed harvesting yield, a forest under a beetle epidemic state has to be cut at a younger age than if the forest were at an endemic beetle state or a no-beetle state.
We consider a nonlinear one-dimensional Stefan problem for a semi-infinite material x>0, with phase change temperature T{sub f}. We assume that the heat capacity and the thermal conductivity satisfy a Storm’s condition. A convective boundary condition and a heat flux over-specified condition on the fixed face x=0 are considered. Unknown thermal coefficients are determined for the free boundary problem and for the associate moving boundary problem and we give sufficient conditions to obtain a parametric representation of a similarity type solution. Moreover, we give formulae for the thermal coefficients in both cases.
In this study, a Legendre spectral tau method is revisited to handle the multi-term time-fractional diffusion equations (MTT-FDEs). An error estimate and rigorous convergence analysis are carried out using some critical theorems. The applicability and accuracy of the solution method are demonstrated by a numerical example to verify the theoretical analysis.
A vector-borne disease model with general incidence rates is proposed and investigated in this paper, where both vector and host are stratied by infection ages in the form of a hyperbolic system of partial differential equations coupled with ordinary differential equations.The existence, uniqueness, nonnegativeness, and boundedness of solution of the model are studied for biologically reasonable purpose. Furthermore, a global threshold dynamics of the system is established by constructing suitable Lyapunov functionals, which is determined by the basic reproduction number R{sub 0}: the infection-free equilibrium is globally asymptotically stable when R{sub 0} < 1 while the endemic equilibrium is globally asymptotically stable when R{sub 0} > 1.
In this paper, the new mass-preserving time second-order explicit–implicit domain decomposition (DD) schemes for solving parabolic equations with various coefficients are proposed. In the schemes, first, the interface values on the interfaces of sub-domains are explicitly calculated, in which two kinds of explicit difference schemes are used. Second, the interior values in sub-domains are determined by the implicit schemes. Third, the interface values are recalculated explicitly to preserve mass. The important feature is that the developed schemes are more than first order in time step and the correction operators at the third step ensure the schemes’ mass conservative. The stability and convergence are analyzed for the two schemes. The optimal error estimates are obtained in the H{sup 1} semi-norm under a stability condition. Numerical experiments confirm the theoretical results.
In this paper, we consider the problem of finding solutions of nonlinear inclusion problems in Banach space. Using convex optimization techniques introduced by Robinson (Numer Math 19:341–347, 1972), a convergence theorem for Kantorovich-like methods is given, which improves the results of Yamamoto (Jpn J Appl Math 3(2):295–313, 1986; Numer Math 51(5):545–557, 1987) and Robinson (Numer Math 19:341–347, 1972). The result is compared with previously known results. Numerical examples further justify the theoretical results.
In this paper, we propose a numerical method to approach the solution of a Fredholm integral equation with a weakly singular kernel by applying the convolution product as a regularization operator and the Fourier series as a projection. Preliminary numerical results show that the order of convergence of the method is better than the one of conventional projection methods.
In this paper, we introduce a new set of wavelets, namely Müntz wavelets. These wavelets are defined based on Müntz Legendre polynomials which have the property of orthogonality. We stabilize the wavelet coefficients by utilizing Jacobi polynomials to construct the modified Müntz wavelets. Next, we present a collocation scheme based on Müntz wavelets for solving fractional differential equations. The fractional derivative is described in the Caputo sense. Some error estimates are given and numerical examples are included to demonstrate the applicability and accuracy of the proposed method. This includes the elastic-dissipative mathematical model of linear mechanical devices known as Boltzman operator.
In this paper, using generalized resolvents and retractions onto shrinking closed subsets, we present new extragradient algorithms for finding the solution of a J-variational-like inequality, the zero points of a family of maximal monotone operators, and the solution set of mixed equilibrium problems for a J-α-inverse-strongly monotone-like operator in a 2-uniformly convex and 2-uniformly smooth Banach space. By introducing the new definitions, we prove strong convergence of generated iterates in the extragradient methods. Using FMINCON optimization toolbox in MATLAB, we give some numerical examples and compare them with several existence results in the literature to illustrate the usability of our results.