Quasi-generalized variables

NASA Technical Reports Server (NTRS)

Baumgarten, J.; Ostermeyer, G. P.

1986-01-01

The numerical solution of a system of differential and algebraic equations is difficult, due to the appearance of numerical instabilities. A method is presented here which permits numerical solutions of such a system to be obtained which satisfy the algebraic constraint equations exactly without reducing the order of the differential equations. The method is demonstrated using examples from mechanics.

Sedimentary Geothermal Feasibility Study: October 2016

DOE Office of Scientific and Technical Information (OSTI.GOV)

Augustine, Chad; Zerpa, Luis

The objective of this project is to analyze the feasibility of commercial geothermal projects using numerical reservoir simulation, considering a sedimentary reservoir with low permeability that requires productivity enhancement. A commercial thermal reservoir simulator (STARS, from Computer Modeling Group, CMG) is used in this work for numerical modeling. In the first stage of this project (FY14), a hypothetical numerical reservoir model was developed, and validated against an analytical solution. The following model parameters were considered to obtain an acceptable match between the numerical and analytical solutions: grid block size, time step and reservoir areal dimensions; the latter related to boundarymore » effects on the numerical solution. Systematic model runs showed that insufficient grid sizing generates numerical dispersion that causes the numerical model to underestimate the thermal breakthrough time compared to the analytic model. As grid sizing is decreased, the model results converge on a solution. Likewise, insufficient reservoir model area introduces boundary effects in the numerical solution that cause the model results to differ from the analytical solution.« less

An approach of traffic signal control based on NLRSQP algorithm

NASA Astrophysics Data System (ADS)

Zou, Yuan-Yang; Hu, Yu

2017-11-01

This paper presents a linear program model with linear complementarity constraints (LPLCC) to solve traffic signal optimization problem. The objective function of the model is to obtain the minimization of total queue length with weight factors at the end of each cycle. Then, a combination algorithm based on the nonlinear least regression and sequence quadratic program (NLRSQP) is proposed, by which the local optimal solution can be obtained. Furthermore, four numerical experiments are proposed to study how to set the initial solution of the algorithm that can get a better local optimal solution more quickly. In particular, the results of numerical experiments show that: The model is effective for different arrival rates and weight factors; and the lower bound of the initial solution is, the better optimal solution can be obtained.

A Numerical Simulation of Scattering from One-Dimensional Inhomogeneous Dielectric Random Surfaces

NASA Technical Reports Server (NTRS)

Sarabandi, Kamal; Oh, Yisok; Ulaby, Fawwaz T.

1996-01-01

In this paper, an efficient numerical solution for the scattering problem of inhomogeneous dielectric rough surfaces is presented. The inhomogeneous dielectric random surface represents a bare soil surface and is considered to be comprised of a large number of randomly positioned dielectric humps of different sizes, shapes, and dielectric constants above an impedance surface. Clods with nonuniform moisture content and rocks are modeled by inhomogeneous dielectric humps and the underlying smooth wet soil surface is modeled by an impedance surface. In this technique, an efficient numerical solution for the constituent dielectric humps over an impedance surface is obtained using Green's function derived by the exact image theory in conjunction with the method of moments. The scattered field from a sample of the rough surface is obtained by summing the scattered fields from all the individual humps of the surface coherently ignoring the effect of multiple scattering between the humps. The statistical behavior of the scattering coefficient sigma(sup 0) is obtained from the calculation of scattered fields of many different realizations of the surface. Numerical results are presented for several different roughnesses and dielectric constants of the random surfaces. The numerical technique is verified by comparing the numerical solution with the solution based on the small perturbation method and the physical optics model for homogeneous rough surfaces. This technique can be used to study the behavior of scattering coefficient and phase difference statistics of rough soil surfaces for which no analytical solution exists.

Numerical simulation of KdV equation by finite difference method

NASA Astrophysics Data System (ADS)

Yokus, A.; Bulut, H.

2018-05-01

In this study, the numerical solutions to the KdV equation with dual power nonlinearity by using the finite difference method are obtained. Discretize equation is presented in the form of finite difference operators. The numerical solutions are secured via the analytical solution to the KdV equation with dual power nonlinearity which is present in the literature. Through the Fourier-Von Neumann technique and linear stable, we have seen that the FDM is stable. Accuracy of the method is analyzed via the L2 and L_{∞} norm errors. The numerical, exact approximations and absolute error are presented in tables. We compare the numerical solutions with the exact solutions and this comparison is supported with the graphic plots. Under the choice of suitable values of parameters, the 2D and 3D surfaces for the used analytical solution are plotted.

GLOBAL PROPERTIES OF FULLY CONVECTIVE ACCRETION DISKS FROM LOCAL SIMULATIONS

DOE Office of Scientific and Technical Information (OSTI.GOV)

Bodo, G.; Ponzo, F.; Rossi, P.

2015-08-01

We present an approach to deriving global properties of accretion disks from the knowledge of local solutions derived from numerical simulations based on the shearing box approximation. The approach consists of a two-step procedure. First, a local solution valid for all values of the disk height is constructed by piecing together an interior solution obtained numerically with an analytical exterior radiative solution. The matching is obtained by assuming hydrostatic balance and radiative equilibrium. Although in principle the procedure can be carried out in general, it simplifies considerably when the interior solution is fully convective. In these cases, the construction ismore » analogous to the derivation of the Hayashi tracks for protostars. The second step consists of piecing together the local solutions at different radii to obtain a global solution. Here we use the symmetry of the solutions with respect to the defining dimensionless numbers—in a way similar to the use of homology relations in stellar structure theory—to obtain the scaling properties of the various disk quantities with radius.« less

Multiresolution representation and numerical algorithms: A brief review

NASA Technical Reports Server (NTRS)

Harten, Amiram

1994-01-01

In this paper we review recent developments in techniques to represent data in terms of its local scale components. These techniques enable us to obtain data compression by eliminating scale-coefficients which are sufficiently small. This capability for data compression can be used to reduce the cost of many numerical solution algorithms by either applying it to the numerical solution operator in order to get an approximate sparse representation, or by applying it to the numerical solution itself in order to reduce the number of quantities that need to be computed.

An Efficient Numerical Approach for Nonlinear Fokker-Planck equations

NASA Astrophysics Data System (ADS)

Otten, Dustin; Vedula, Prakash

2009-03-01

Fokker-Planck equations which are nonlinear with respect to their probability densities that occur in many nonequilibrium systems relevant to mean field interaction models, plasmas, classical fermions and bosons can be challenging to solve numerically. To address some underlying challenges in obtaining numerical solutions, we propose a quadrature based moment method for efficient and accurate determination of transient (and stationary) solutions of nonlinear Fokker-Planck equations. In this approach the distribution function is represented as a collection of Dirac delta functions with corresponding quadrature weights and locations, that are in turn determined from constraints based on evolution of generalized moments. Properties of the distribution function can be obtained by solution of transport equations for quadrature weights and locations. We will apply this computational approach to study a wide range of problems, including the Desai-Zwanzig Model (for nonlinear muscular contraction) and multivariate nonlinear Fokker-Planck equations describing classical fermions and bosons, and will also demonstrate good agreement with results obtained from Monte Carlo and other standard numerical methods.

A meshless method using radial basis functions for numerical solution of the two-dimensional KdV-Burgers equation

NASA Astrophysics Data System (ADS)

Zabihi, F.; Saffarian, M.

2016-07-01

The aim of this article is to obtain the numerical solution of the two-dimensional KdV-Burgers equation. We construct the solution by using a different approach, that is based on using collocation points. The solution is based on using the thin plate splines radial basis function, which builds an approximated solution with discretizing the time and the space to small steps. We use a predictor-corrector scheme to avoid solving the nonlinear system. The results of numerical experiments are compared with analytical solutions to confirm the accuracy and efficiency of the presented scheme.

A Semi-Analytical Solution to Time Dependent Groundwater Flow Equation Incorporating Stream-Wetland-Aquifer Interactions

NASA Astrophysics Data System (ADS)

Boyraz, Uǧur; Melek Kazezyılmaz-Alhan, Cevza

2017-04-01

Groundwater is a vital element of hydrologic cycle and the analytical & numerical solutions of different forms of groundwater flow equations play an important role in understanding the hydrological behavior of subsurface water. The interaction between groundwater and surface water bodies can be determined using these solutions. In this study, new hypothetical approaches are implemented to groundwater flow system in order to contribute to the studies on surface water/groundwater interactions. A time dependent problem is considered in a 2-dimensional stream-wetland-aquifer system. The sloped stream boundary is used to represent the interaction between stream and aquifer. The rest of the aquifer boundaries are assumed as no-flux boundary. In addition, a wetland is considered as a surface water body which lies over the whole aquifer. The effect of the interaction between the wetland and the aquifer is taken into account with a source/sink term in the groundwater flow equation and the interaction flow is calculated by using Darcy's approach. A semi-analytical solution is developed for the 2-dimensional groundwater flow equation in 5 steps. First, Laplace and Fourier cosine transforms are employed to obtain the general solution in Fourier and Laplace domain. Then, the initial and boundary conditions are applied to obtain the particular solution. Finally, inverse Fourier transform is carried out analytically and inverse Laplace transform is carried out numerically to obtain the final solution in space and time domain, respectively. In order to verify the semi-analytical solution, an explicit finite difference algorithm is developed and analytical and numerical solutions are compared for synthetic examples. The comparison of the analytical and numerical solutions shows that the analytical solution gives accurate results.

Exact Closed-form Solutions for Lamb's Problem

NASA Astrophysics Data System (ADS)

Feng, Xi; Zhang, Haiming

2018-04-01

In this article, we report on an exact closed-form solution for the displacement at the surface of an elastic half-space elicited by a buried point source that acts at some point underneath that surface. This is commonly referred to as the 3-D Lamb's problem, for which previous solutions were restricted to sources and receivers placed at the free surface. By means of the reciprocity theorem, our solution should also be valid as a means to obtain the displacements at interior points when the source is placed at the free surface. We manage to obtain explicit results by expressing the solution in terms of elementary algebraic expression as well as elliptic integrals. We anchor our developments on Poisson's ratio 0.25 starting from Johnson's (1974) integral solutions which must be computed numerically. In the end, our closed-form results agree perfectly with the numerical results of Johnson (1974), which strongly confirms the correctness of our explicit formulas. It is hoped that in due time, these formulas may constitute a valuable canonical solution that will serve as a yardstick against which other numerical solutions can be compared and measured.

Exact closed-form solutions for Lamb's problem

NASA Astrophysics Data System (ADS)

Feng, Xi; Zhang, Haiming

2018-07-01

In this paper, we report on an exact closed-form solution for the displacement at the surface of an elastic half-space elicited by a buried point source that acts at some point underneath that surface. This is commonly referred to as the 3-D Lamb's problem for which previous solutions were restricted to sources and receivers placed at the free surface. By means of the reciprocity theorem, our solution should also be valid as a means to obtain the displacements at interior points when the source is placed at the free surface. We manage to obtain explicit results by expressing the solution in terms of elementary algebraic expression as well as elliptic integrals. We anchor our developments on Poisson's ratio 0.25 starting from Johnson's integral solutions which must be computed numerically. In the end, our closed-form results agree perfectly with the numerical results of Johnson, which strongly confirms the correctness of our explicit formulae. It is hoped that in due time, these formulae may constitute a valuable canonical solution that will serve as a yardstick against which other numerical solutions can be compared and measured.

An iterative transformation procedure for numerical solution of flutter and similar characteristics-value problems

NASA Technical Reports Server (NTRS)

Gossard, Myron L

1952-01-01

An iterative transformation procedure suggested by H. Wielandt for numerical solution of flutter and similar characteristic-value problems is presented. Application of this procedure to ordinary natural-vibration problems and to flutter problems is shown by numerical examples. Comparisons of computed results with experimental values and with results obtained by other methods of analysis are made.

Projection scheme for a reflected stochastic heat equation with additive noise

NASA Astrophysics Data System (ADS)

Higa, Arturo Kohatsu; Pettersson, Roger

2005-02-01

We consider a projection scheme as a numerical solution of a reflected stochastic heat equation driven by a space-time white noise. Convergence is obtained via a discrete contraction principle and known convergence results for numerical solutions of parabolic variational inequalities.

Unsteady magnetohydrodynamics mixed convection flow in a rotating medium with double diffusion

DOE Office of Scientific and Technical Information (OSTI.GOV)

Jiann, Lim Yeou; Ismail, Zulkhibri; Khan, Ilyas

2015-05-15

Exact solutions of an unsteady Magnetohydrodynamics (MHD) flow over an impulsively started vertical plate in a rotating medium are presented. The effects of thermal radiative and thermal diffusion on the fluid flow are also considered. The governing equations are modelled and solved for velocity, temperature and concentration using Laplace transforms technique. Expressions of velocity, temperature and concentration profiles are obtained and their numerical results are presented graphically. Skin friction, Sherwood number and Nusselt number are also computed and presented in tabular forms. The determined solutions can generate a large class of solutions as special cases corresponding to different motions withmore » technical relevance. The results obtained herein may be used to verify the validation of obtained numerical solutions for more complicated fluid flow problems.« less

Numerical modelling of the Black Sea eigen-oscillations on a curvilinear boundary fitted coordinate system

NASA Astrophysics Data System (ADS)

Koychev Demirov, Encho

1994-12-01

The paper presents a numerical solution of barotropic and two-layer eigen-oscillation problems for the Black Sea on a boundary fitted coordinate system. This solution is compared with model and empirical data obtained by other workers. Frequencies of the eigen-oscillations found by the numerical solution of spectral problem are compared with the data obtained by spectral analysis of the sea-level oscillations measured near the town of Achtopol and Cape Irakli in stormy sea on 17-21 February 1979. Extreme oscillations of the sea-level result from resonant amplifications of three eigenmodes of the Black Sea of 68.3 -1, 36.6 -1 and 27.3 -1 cycles h -1 frequency.

Numerical solution of a coupled pair of elliptic equations from solid state electronics

NASA Technical Reports Server (NTRS)

Phillips, T. N.

1983-01-01

Iterative methods are considered for the solution of a coupled pair of second order elliptic partial differential equations which arise in the field of solid state electronics. A finite difference scheme is used which retains the conservative form of the differential equations. Numerical solutions are obtained in two ways, by multigrid and dynamic alternating direction implicit methods. Numerical results are presented which show the multigrid method to be an efficient way of solving this problem.

Second kind Chebyshev operational matrix algorithm for solving differential equations of Lane-Emden type

NASA Astrophysics Data System (ADS)

Doha, E. H.; Abd-Elhameed, W. M.; Youssri, Y. H.

2013-10-01

In this paper, we present a new second kind Chebyshev (S2KC) operational matrix of derivatives. With the aid of S2KC, an algorithm is described to obtain numerical solutions of a class of linear and nonlinear Lane-Emden type singular initial value problems (IVPs). The idea of obtaining such solutions is essentially based on reducing the differential equation with its initial conditions to a system of algebraic equations. Two illustrative examples concern relevant physical problems (the Lane-Emden equations of the first and second kind) are discussed to demonstrate the validity and applicability of the suggested algorithm. Numerical results obtained are comparing favorably with the analytical known solutions.

A mathematical solution for the parameters of three interfering resonances

NASA Astrophysics Data System (ADS)

Han, X.; Shen, C. P.

2018-04-01

The multiple-solution problem in determining the parameters of three interfering resonances from a fit to an experimentally measured distribution is considered from a mathematical viewpoint. It is shown that there are four numerical solutions for a fit with three coherent Breit-Wigner functions. Although explicit analytical formulae cannot be derived in this case, we provide some constraint equations between the four solutions. For the cases of nonrelativistic and relativistic Breit-Wigner forms of amplitude functions, a numerical method is provided to derive the other solutions from that already obtained, based on the obtained constraint equations. In real experimental measurements with more complicated amplitude forms similar to Breit-Wigner functions, the same method can be deduced and performed to get numerical solutions. The good agreement between the solutions found using this mathematical method and those directly from the fit verifies the correctness of the constraint equations and mathematical methodology used. Supported by National Natural Science Foundation of China (NSFC) (11575017, 11761141009), the Ministry of Science and Technology of China (2015CB856701) and the CAS Center for Excellence in Particle Physics (CCEPP)

Numerical Solution of the Electron Transport Equation in the Upper Atmosphere

DOE Office of Scientific and Technical Information (OSTI.GOV)

Woods, Mark Christopher; Holmes, Mark; Sailor, William C

A new approach for solving the electron transport equation in the upper atmosphere is derived. The problem is a very stiff boundary value problem, and to obtain an accurate numerical solution, matrix factorizations are used to decouple the fast and slow modes. A stable finite difference method is applied to each mode. This solver is applied to a simplifieed problem for which an exact solution exists using various versions of the boundary conditions that might arise in a natural auroral display. The numerical and exact solutions are found to agree with each other to at least two significant digits.

An extension of the Derrida-Lebowitz-Speer-Spohn equation

NASA Astrophysics Data System (ADS)

Bordenave, Charles; Germain, Pierre; Trogdon, Thomas

2015-12-01

We show how the derivation of the Derrida-Lebowitz-Speer-Spohn equation can be prolonged to obtain a new equation, generalizing the models obtained in the paper by these authors. We then investigate its properties from both an analytical and numerical perspective. Specifically, a numerical method is presented to approximate solutions of the prolonged equation. Using this method, we investigate the relationship between the solutions of the prolonged equation and the Tracy-Widom GOE distribution.

Analytical and Numerical solutions of a nonlinear alcoholism model via variable-order fractional differential equations

NASA Astrophysics Data System (ADS)

Gómez-Aguilar, J. F.

2018-03-01

In this paper, we analyze an alcoholism model which involves the impact of Twitter via Liouville-Caputo and Atangana-Baleanu-Caputo fractional derivatives with constant- and variable-order. Two fractional mathematical models are considered, with and without delay. Special solutions using an iterative scheme via Laplace and Sumudu transform were obtained. We studied the uniqueness and existence of the solutions employing the fixed point postulate. The generalized model with variable-order was solved numerically via the Adams method and the Adams-Bashforth-Moulton scheme. Stability and convergence of the numerical solutions were presented in details. Numerical examples of the approximate solutions are provided to show that the numerical methods are computationally efficient. Therefore, by including both the fractional derivatives and finite time delays in the alcoholism model studied, we believe that we have established a more complete and more realistic indicator of alcoholism model and affect the spread of the drinking.

Stochastic coalescence in finite systems: an algorithm for the numerical solution of the multivariate master equation.

NASA Astrophysics Data System (ADS)

Alfonso, Lester; Zamora, Jose; Cruz, Pedro

2015-04-01

The stochastic approach to coagulation considers the coalescence process going in a system of a finite number of particles enclosed in a finite volume. Within this approach, the full description of the system can be obtained from the solution of the multivariate master equation, which models the evolution of the probability distribution of the state vector for the number of particles of a given mass. Unfortunately, due to its complexity, only limited results were obtained for certain type of kernels and monodisperse initial conditions. In this work, a novel numerical algorithm for the solution of the multivariate master equation for stochastic coalescence that works for any type of kernels and initial conditions is introduced. The performance of the method was checked by comparing the numerically calculated particle mass spectrum with analytical solutions obtained for the constant and sum kernels, with an excellent correspondence between the analytical and numerical solutions. In order to increase the speedup of the algorithm, software parallelization techniques with OpenMP standard were used, along with an implementation in order to take advantage of new accelerator technologies. Simulations results show an important speedup of the parallelized algorithms. This study was funded by a grant from Consejo Nacional de Ciencia y Tecnologia de Mexico SEP-CONACYT CB-131879. The authors also thanks LUFAC® Computacion SA de CV for CPU time and all the support provided.

Sensitivity of inelastic response to numerical integration of strain energy. [for cantilever beam

NASA Technical Reports Server (NTRS)

Kamat, M. P.

1976-01-01

The exact solution to the quasi-static, inelastic response of a cantilever beam of rectangular cross section subjected to a bending moment at the tip is obtained. The material of the beam is assumed to be linearly elastic-linearly strain-hardening. This solution is then compared with three different numerical solutions of the same problem obtained by minimizing the total potential energy using Gaussian quadratures of two different orders and a Newton-Cotes scheme for integrating the strain energy of deformation. Significant differences between the exact dissipative strain energy and its numerical counterpart are emphasized. The consequence of this on the nonlinear transient responses of a beam with solid cross section and that of a thin-walled beam on elastic supports under impulsive loads are examined.

Analytical guidance law development for aerocapture at Mars

NASA Technical Reports Server (NTRS)

Calise, A. J.

1992-01-01

During the first part of this reporting period research has concentrated on performing a detailed evaluation, to zero order, of the guidance algorithm developed in the first period taking the numerical approach developed in the third period. A zero order matched asymptotic expansion (MAE) solution that closely satisfies a set of 6 implicit equations in 6 unknowns to an accuracy of 10(exp -10), was evaluated. Guidance law implementation entails treating the current state as a new initial state and repetitively solving the MAE problem to obtain the feedback controls. A zero order guided solution was evaluated and compared with optimal solution that was obtained by numerical methods. Numerical experience shows that the zero order guided solution is close to optimal solution, and that the zero order MAE outer solution plays a critical role in accounting for the variations in Loh's term near the exit phase of the maneuver. However, the deficiency that remains in several of the critical variables indicates the need for a first order correction. During the second part of this period, methods for computing a first order correction were explored.

Numerical solution of fluid flow and heat tranfer problems with surface radiation

NASA Technical Reports Server (NTRS)

Ahuja, S.; Bhatia, K.

1995-01-01

This paper presents a numerical scheme, based on the finite element method, to solve strongly coupled fluid flow and heat transfer problems. The surface radiation effect for gray, diffuse and isothermal surfaces is considered. A procedure for obtaining the view factors between the radiating surfaces is discussed. The overall solution strategy is verified by comparing the available results with those obtained using this approach. An analysis of a thermosyphon is undertaken and the effect of considering the surface radiation is clearly explained.

Numerical analysis of the transient response of an axisymmetric ablative char layer considering internal flow effects

NASA Technical Reports Server (NTRS)

Pittman, C. M.; Howser, L. M.

1972-01-01

The differential equations governing the transient response of the char layer of an ablating axisymmetric body, internal pyrolysis gas flow effects being considered, have been derived. These equations have been expanded into finite difference form and programed for numerical solution on a digital computer. Numerical results compare favorably with simplified exact solutions. The complete numerical analysis was used to obtain solutions for two representative body shapes subjected to a typical entry heating environment. Pronounced effects of the lateral flow of pyrolysis gases on the mass flow field within the char layer and the associated surface and pyrolysis interface recession rates are shown.

A numerical solution of a singular boundary value problem arising in boundary layer theory.

PubMed

Hu, Jiancheng

2016-01-01

In this paper, a second-order nonlinear singular boundary value problem is presented, which is equivalent to the well-known Falkner-Skan equation. And the one-dimensional third-order boundary value problem on interval [Formula: see text] is equivalently transformed into a second-order boundary value problem on finite interval [Formula: see text]. The finite difference method is utilized to solve the singular boundary value problem, in which the amount of computational effort is significantly less than the other numerical methods. The numerical solutions obtained by the finite difference method are in agreement with those obtained by previous authors.

Advanced numerical methods for three dimensional two-phase flow calculations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Toumi, I.; Caruge, D.

1997-07-01

This paper is devoted to new numerical methods developed for both one and three dimensional two-phase flow calculations. These methods are finite volume numerical methods and are based on the use of Approximate Riemann Solvers concepts to define convective fluxes versus mean cell quantities. The first part of the paper presents the numerical method for a one dimensional hyperbolic two-fluid model including differential terms as added mass and interface pressure. This numerical solution scheme makes use of the Riemann problem solution to define backward and forward differencing to approximate spatial derivatives. The construction of this approximate Riemann solver uses anmore » extension of Roe`s method that has been successfully used to solve gas dynamic equations. As far as the two-fluid model is hyperbolic, this numerical method seems very efficient for the numerical solution of two-phase flow problems. The scheme was applied both to shock tube problems and to standard tests for two-fluid computer codes. The second part describes the numerical method in the three dimensional case. The authors discuss also some improvements performed to obtain a fully implicit solution method that provides fast running steady state calculations. Such a scheme is not implemented in a thermal-hydraulic computer code devoted to 3-D steady-state and transient computations. Some results obtained for Pressurised Water Reactors concerning upper plenum calculations and a steady state flow in the core with rod bow effect evaluation are presented. In practice these new numerical methods have proved to be stable on non staggered grids and capable of generating accurate non oscillating solutions for two-phase flow calculations.« less

Forebody and base region real gas flow in severe planetary entry by a factored implicit numerical method. II - Equilibrium reactive gas

NASA Technical Reports Server (NTRS)

Davy, W. C.; Green, M. J.; Lombard, C. K.

1981-01-01

The factored-implicit, gas-dynamic algorithm has been adapted to the numerical simulation of equilibrium reactive flows. Changes required in the perfect gas version of the algorithm are developed, and the method of coupling gas-dynamic and chemistry variables is discussed. A flow-field solution that approximates a Jovian entry case was obtained by this method and compared with the same solution obtained by HYVIS, a computer program much used for the study of planetary entry. Comparison of surface pressure distribution and stagnation line shock-layer profiles indicates that the two solutions agree well.

Traveling-Wave Solutions of the Kolmogorov-Petrovskii-Piskunov Equation

NASA Astrophysics Data System (ADS)

Pikulin, S. V.

2018-02-01

We consider quasi-stationary solutions of a problem without initial conditions for the Kolmogorov-Petrovskii-Piskunov (KPP) equation, which is a quasilinear parabolic one arising in the modeling of certain reaction-diffusion processes in the theory of combustion, mathematical biology, and other areas of natural sciences. A new efficiently numerically implementable analytical representation is constructed for self-similar plane traveling-wave solutions of the KPP equation with a special right-hand side. Sufficient conditions for an auxiliary function involved in this representation to be analytical for all values of its argument, including the endpoints, are obtained. Numerical results are obtained for model examples.

An analysis of a charring ablator with thermal nonequilibrium, chemical kinetics, and mass transfer

NASA Technical Reports Server (NTRS)

Clark, R. K.

1973-01-01

The differential equations governing the transient response of a one-dimensional ablative thermal protection system are presented for thermal nonequilibrium between the pyrolysis gases and the char layer and with finite rate chemical reactions occurring. The system consists of three layers (the char layer, the uncharred layer, and an optical insulation layer) with concentrated heat sinks at the back surface and between the second and third layers. The equations are solved numerically by using a modified implicit finite difference scheme to obtain solutions for the thickness of the charred and uncharred layers, surface recession and pyrolysis rates, solid temperatures, porosity profiles, and profiles of pyrolysis-gas temperature, pressure, composition, and flow rate. Good agreement is obtained between numerical results and exact solutions for a number of simplified cases. The complete numerical analysis is used to obtain solutions for an ablative system subjected to a constant heating environment. Effects of thermal, chemical, and mass transfer processes are shown.

On the limits of numerical astronomical solutions used in paleoclimate studies

NASA Astrophysics Data System (ADS)

Zeebe, Richard E.

2017-04-01

Numerical solutions of the equations of the Solar System estimate Earth's orbital parameters in the past and represent the backbone of cyclostratigraphy and astrochronology, now widely applied in geology and paleoclimatology. Given one numerical realization of a Solar System model (i.e., obtained using one code or integrator package), various parameters determine the properties of the solution and usually limit its validity to a certain time period. Such limitations are denoted here as "internal" and include limitations due to (i) the underlying physics/physical model and (ii) numerics. The physics include initial coordinates and velocities of Solar System bodies, treatment of the Moon and asteroids, the Sun's quadrupole moment, and the intrinsic dynamics of the Solar System itself, i.e., its chaotic nature. Numerical issues include solver algorithm, numerical accuracy (e.g., time step), and round-off errors. At present, internal limitations seem to restrict the validity of astronomical solutions to perhaps the past 50 or 60 myr. However, little is currently known about "external" limitations, that is, how do different numerical realizations compare, say, between different investigators using different codes and integrators? Hitherto only two solutions for Earth's eccentricity appear to be used in paleoclimate studies, provided by two different groups that integrated the full Solar System equations over the past >100 myr (Laskar and coworkers and Varadi et al. 2003). In this contribution, I will present results from new Solar System integrations for Earth's eccentricity obtained using the integrator package HNBody (Rauch and Hamilton 2002). I will discuss the various internal limitations listed above within the framework of the present simulations. I will also compare the results to the existing solutions, the details of which are still being sorted out as several simulations are still running at the time of writing.

Hypersonic shock structure with Burnett terms in the viscous stress and heat flux

NASA Technical Reports Server (NTRS)

Chapman, Dean R.; Fiscko, Kurt A.

1988-01-01

The continuum Navier-Stokes and Burnett equations are solved for one-dimensional shock structure in various monatomic gases. A new numerical method is employed which utilizes the complete time-dependent continuum equations and obtains the steady-state shock structure by allowing the system to relax from arbitrary initial conditions. Included is discussion of numerical difficulties encountered when solving the Burnett equations. Continuum solutions are compared to those obtained utilizing the Direct Simulation Monte Carlo method. Shock solutions are obtained for a hard sphere gas and for argon from Mach 1.3 to Mach 50. Solutions for a Maxwellian gas are obtained from Mach 1.3 to Mach 3.8. It is shown that the Burnett equations yield shock structure solutions in much closer agreement to both Monte Carlo and experimental results than do the Navier-Stokes equations. Shock density thickness, density asymmetry, and density-temperature separation are all more accurately predicted by the Burnett equations than by the Navier-Stokes equations.

A Grobner Basis Solution for Lightning Ground Flash Fraction Retrieval

NASA Technical Reports Server (NTRS)

Solakiewicz, Richard; Attele, Rohan; Koshak, William

2011-01-01

A Bayesian inversion method was previously introduced for retrieving the fraction of ground flashes in a set of flashes observed from a (low earth orbiting or geostationary) satellite lightning imager. The method employed a constrained mixed exponential distribution model to describe the lightning optical measurements. To obtain the optimum model parameters, a scalar function was minimized by a numerical method. In order to improve this optimization, we introduce a Grobner basis solution to obtain analytic representations of the model parameters that serve as a refined initialization scheme to the numerical optimization. Using the Grobner basis, we show that there are exactly 2 solutions involving the first 3 moments of the (exponentially distributed) data. When the mean of the ground flash optical characteristic (e.g., such as the Maximum Group Area, MGA) is larger than that for cloud flashes, then a unique solution can be obtained.

Analytical solutions for benchmarking cold regions subsurface water flow and energy transport models: one-dimensional soil thaw with conduction and advection

USGS Publications Warehouse

Kurylyk, Barret L.; McKenzie, Jeffrey M; MacQuarrie, Kerry T. B.; Voss, Clifford I.

2014-01-01

Numerous cold regions water flow and energy transport models have emerged in recent years. Dissimilarities often exist in their mathematical formulations and/or numerical solution techniques, but few analytical solutions exist for benchmarking flow and energy transport models that include pore water phase change. This paper presents a detailed derivation of the Lunardini solution, an approximate analytical solution for predicting soil thawing subject to conduction, advection, and phase change. Fifteen thawing scenarios are examined by considering differences in porosity, surface temperature, Darcy velocity, and initial temperature. The accuracy of the Lunardini solution is shown to be proportional to the Stefan number. The analytical solution results obtained for soil thawing scenarios with water flow and advection are compared to those obtained from the finite element model SUTRA. Three problems, two involving the Lunardini solution and one involving the classic Neumann solution, are recommended as standard benchmarks for future model development and testing.

On the implementation of an accurate and efficient solver for convection-diffusion equations

NASA Astrophysics Data System (ADS)

Wu, Chin-Tien

In this dissertation, we examine several different aspects of computing the numerical solution of the convection-diffusion equation. The solution of this equation often exhibits sharp gradients due to Dirichlet outflow boundaries or discontinuities in boundary conditions. Because of the singular-perturbed nature of the equation, numerical solutions often have severe oscillations when grid sizes are not small enough to resolve sharp gradients. To overcome such difficulties, the streamline diffusion discretization method can be used to obtain an accurate approximate solution in regions where the solution is smooth. To increase accuracy of the solution in the regions containing layers, adaptive mesh refinement and mesh movement based on a posteriori error estimations can be employed. An error-adapted mesh refinement strategy based on a posteriori error estimations is also proposed to resolve layers. For solving the sparse linear systems that arise from discretization, goemetric multigrid (MG) and algebraic multigrid (AMG) are compared. In addition, both methods are also used as preconditioners for Krylov subspace methods. We derive some convergence results for MG with line Gauss-Seidel smoothers and bilinear interpolation. Finally, while considering adaptive mesh refinement as an integral part of the solution process, it is natural to set a stopping tolerance for the iterative linear solvers on each mesh stage so that the difference between the approximate solution obtained from iterative methods and the finite element solution is bounded by an a posteriori error bound. Here, we present two stopping criteria. The first is based on a residual-type a posteriori error estimator developed by Verfurth. The second is based on an a posteriori error estimator, using local solutions, developed by Kay and Silvester. Our numerical results show the refined mesh obtained from the iterative solution which satisfies the second criteria is similar to the refined mesh obtained from the finite element solution.

Transient Effects in Planar Solidification of Dilute Binary Alloys

NASA Technical Reports Server (NTRS)

Mazuruk, Konstantin; Volz, Martin P.

2008-01-01

The initial transient during planar solidification of dilute binary alloys is studied in the framework of the boundary integral method that leads to the non-linear Volterra integral governing equation. An analytical solution of this equation is obtained for the case of a constant growth rate which constitutes the well-known Tiller's formula for the solute transient. The more physically relevant, constant ramping down temperature case has been studied both numerically and analytically. In particular, an asymptotic analytical solution is obtained for the initial transient behavior. A numerical technique to solve the non-linear Volterra equation is developed and the solution is obtained for a family of the governing parameters. For the rapid solidification condition, growth rate spikes have been observed even for the infinite kinetics model. When recirculating fluid flow is included into the analysis, the spike feature is dramatically diminished. Finally, we have investigated planar solidification with a fluctuating temperature field as a possible mechanism for frequently observed solute trapping bands.

Study of periodic motions of a satellite with a magnetic damper

NASA Technical Reports Server (NTRS)

Sadov, Y. A.; Teterin, A. D.

1979-01-01

The motion of a satellite with a magnetic damper in the plane of a circular polar orbit is studied. The asymptotics of periodic solutions are constructed for a satellite close to axisymmetric and the radius of convergence is evaluated for the power series obtained. In a broad range of values of parameters, a periodic solution is obtained by numerical integration of equations of motion of the satellite. The asymptotics of a bifurcated curve obtained (the curve on which origin of a pair of periodic solutions occurs) in the space of the parameters agrees well with the results of numerical computation with all physical values of these parameters. A breakdown is made of the space of the initial data of phase variables in the field of effect of different types of periodic motion.

Axisymmetric deformation in a micropolar thermoelastic medium under fractional order theory of thermoelasticity

NASA Astrophysics Data System (ADS)

Kumar, Rajneesh; Singh, Kulwinder; Pathania, Devinder Singh

2017-07-01

The purpose of this paper is to study the variations in temperature, radial and normal displacement, normal stress, shear stress and couple stress in a micropolar thermoelastic solid in the context of fractional order theory of thermoelasticity. Eigen value approach together with Laplace and Hankel transforms are employed to obtain the general solution of the problem. The field variables corresponding to different fractional order theories of thermoelasticity have been obtained in the transformed domain. The general solution is applied to an infinite space subjected to a concentrated load at the origin. To obtained solution in the physical domain numerical inversion technique has been applied and numerically computed results are depicted graphically to analyze the effects of fractional order parameter on the field variables.

Spatiotemporal Airy Ince-Gaussian wave packets in strongly nonlocal nonlinear media.

PubMed

Peng, Xi; Zhuang, Jingli; Peng, Yulian; Li, DongDong; Zhang, Liping; Chen, Xingyu; Zhao, Fang; Deng, Dongmei

2018-03-08

The self-accelerating Airy Ince-Gaussian (AiIG) and Airy helical Ince-Gaussian (AihIG) wave packets in strongly nonlocal nonlinear media (SNNM) are obtained by solving the strongly nonlocal nonlinear Schrödinger equation. For the first time, the propagation properties of three dimensional localized AiIG and AihIG breathers and solitons in the SNNM are demonstrated, these spatiotemporal wave packets maintain the self-accelerating and approximately non-dispersion properties in temporal dimension, periodically oscillating (breather state) or steady (soliton state) in spatial dimension. In particular, their numerical experiments of spatial intensity distribution, numerical simulations of spatiotemporal distribution, as well as the transverse energy flow and the angular momentum in SNNM are presented. Typical examples of the obtained solutions are based on the ratio between the input power and the critical power, the ellipticity and the strong nonlocality parameter. The comparisons of analytical solutions with numerical simulations and numerical experiments of the AiIG and AihIG optical solitons show that the numerical results agree well with the analytical solutions in the case of strong nonlocality.

Numerical analysis of singular solutions of two-dimensional problems of asymmetric elasticity

NASA Astrophysics Data System (ADS)

Korepanov, V. V.; Matveenko, V. P.; Fedorov, A. Yu.; Shardakov, I. N.

2013-07-01

An algorithm for the numerical analysis of singular solutions of two-dimensional problems of asymmetric elasticity is considered. The algorithm is based on separation of a power-law dependence from the finite-element solution in a neighborhood of singular points in the domain under study, where singular solutions are possible. The obtained power-law dependencies allow one to conclude whether the stresses have singularities and what the character of these singularities is. The algorithm was tested for problems of classical elasticity by comparing the stress singularity exponents obtained by the proposed method and from known analytic solutions. Problems with various cases of singular points, namely, body surface points at which either the smoothness of the surface is violated, or the type of boundary conditions is changed, or distinct materials are in contact, are considered as applications. The stress singularity exponents obtained by using the models of classical and asymmetric elasticity are compared. It is shown that, in the case of cracks, the stress singularity exponents are the same for the elasticity models under study, but for other cases of singular points, the stress singularity exponents obtained on the basis of asymmetric elasticity have insignificant quantitative distinctions from the solutions of the classical elasticity.

NONLINEAR AND FIBER OPTICS: Self-similar solution obtained by self-focusing of annular laser beams

NASA Astrophysics Data System (ADS)

Azimov, B. S.; Platonenko, Viktor T.; Sagatov, M. M.

1991-03-01

A numerical modeling is reported of steady-state self-focusing of an annular beam with thin "walls." An approximate similar solution is found to describe well the relationships observed in the numerical experiment for a special selection of the input parameters of the beam. This solution is used to estimate the focal length. Such self-similar self-focusing is shown to affect the whole power of the beam.

Analytical Solutions for the Surface States of Bi1-xSbx (0 ≤ x ≲ 0.1)

NASA Astrophysics Data System (ADS)

Fuseya, Yuki; Fukuyama, Hidetoshi

2018-04-01

Analytical solutions for the surface state (SS) of an extended Wolff Hamiltonian, which is a common Hamiltonian for strongly spin-orbit coupled systems, are obtained both for semi-infinite and finite-thickness boundary conditions. For the semi-infinite system, there are two types of SS solutions: (I-a) linearly crossing SSs in the direct bulk band gap, and (I-b) SSs with linear dispersions entering the bulk conduction or valence bands away from the band edge. For the finite-thickness system, a gap opens in the SS of solution I-a. Numerical solutions for the SS are also obtained based on the tight-binding model of Liu and Allen [Phys. Rev. B 52, 1566 (1995)] for Bi1-xSbx (0 ≤ x ≤ 0.1). A perfect correspondence between the analytic and numerical solutions is obtained around the \\bar{M} point including their thickness dependence. This is the first time that the character of the SS numerically obtained is identified with the help of analytical solutions. The size of the gap for I-a SS can be larger than that of bulk band gap even for a "thick" films ( ≲ 200 bilayers ≃ 80 nm) of pure bismuth. Consequently, in such a film of Bi1-xSbx, there is no apparent change in the SSs through the band inversion at x ≃ 0.04, even though the nature of the SS is changed from solution I-a to I-b. Based on our theoretical results, the experimental results on the SS of Bi1-xSbx (0 ≤ x ≲ 0.1) are discussed.

Numerical Simulations of Multidimensional Flows in Presence of either Strong Shocks or Strong Gravitational Fields

NASA Astrophysics Data System (ADS)

Font, J. A.; Ibanez, J. M.; Marti, J. M.

1993-04-01

Some numerical solutions via local characteristic approach have been obtained describing multidimensional flows. These solutions have been used as tests of a two- dimensional code which extends some high-resolution shock-captunng methods, designed recently to solve nonlinear hyperbolic systems of conservation laws. K words: HYDRODYNAMICS - BLACK HOLE - RELATIVITY - SHOCK WAVES

Combined structures-controls optimization of lattice trusses

NASA Technical Reports Server (NTRS)

Balakrishnan, A. V.

1991-01-01

The role that distributed parameter model can play in CSI is demonstrated, in particular in combined structures controls optimization problems of importance in preliminary design. Closed form solutions can be obtained for performance criteria such as rms attitude error, making possible analytical solutions of the optimization problem. This is in contrast to the need for numerical computer solution involving the inversion of large matrices in traditional finite element model (FEM) use. Another advantage of the analytic solution is that it can provide much needed insight into phenomena that can otherwise be obscured or difficult to discern from numerical computer results. As a compromise in level of complexity between a toy lab model and a real space structure, the lattice truss used in the EPS (Earth Pointing Satellite) was chosen. The optimization problem chosen is a generic one: of minimizing the structure mass subject to a specified stability margin and to a specified upper bond on the rms attitude error, using a co-located controller and sensors. Standard FEM treating each bar as a truss element is used, while the continuum model is anisotropic Timoshenko beam model. Performance criteria are derived for each model, except that for the distributed parameter model, explicit closed form solutions was obtained. Numerical results obtained by the two model show complete agreement.

Computation of steady and unsteady quasi-one-dimensional viscous/inviscid interacting internal flows at subsonic, transonic, and supersonic Mach numbers

NASA Technical Reports Server (NTRS)

Swafford, Timothy W.; Huddleston, David H.; Busby, Judy A.; Chesser, B. Lawrence

1992-01-01

Computations of viscous-inviscid interacting internal flowfields are presented for steady and unsteady quasi-one-dimensional (Q1D) test cases. The unsteady Q1D Euler equations are coupled with integral boundary-layer equations for unsteady, two-dimensional (planar or axisymmetric), turbulent flow over impermeable, adiabatic walls. The coupling methodology differs from that used in most techniques reported previously in that the above mentioned equation sets are written as a complete system and solved simultaneously; that is, the coupling is carried out directly through the equations as opposed to coupling the solutions of the different equation sets. Solutions to the coupled system of equations are obtained using both explicit and implicit numerical schemes for steady subsonic, steady transonic, and both steady and unsteady supersonic internal flowfields. Computed solutions are compared with measurements as well as Navier-Stokes and inverse boundary-layer methods. An analysis of the eigenvalues of the coefficient matrix associated with the quasi-linear form of the coupled system of equations indicates the presence of complex eigenvalues for certain flow conditions. It is concluded that although reasonable solutions can be obtained numerically, these complex eigenvalues contribute to the overall difficulty in obtaining numerical solutions to the coupled system of equations.

Second-order numerical methods for multi-term fractional differential equations: Smooth and non-smooth solutions

NASA Astrophysics Data System (ADS)

Zeng, Fanhai; Zhang, Zhongqiang; Karniadakis, George Em

2017-12-01

Starting with the asymptotic expansion of the error equation of the shifted Gr\\"{u}nwald--Letnikov formula, we derive a new modified weighted shifted Gr\\"{u}nwald--Letnikov (WSGL) formula by introducing appropriate correction terms. We then apply one special case of the modified WSGL formula to solve multi-term fractional ordinary and partial differential equations, and we prove the linear stability and second-order convergence for both smooth and non-smooth solutions. We show theoretically and numerically that numerical solutions up to certain accuracy can be obtained with only a few correction terms. Moreover, the correction terms can be tuned according to the fractional derivative orders without explicitly knowing the analytical solutions. Numerical simulations verify the theoretical results and demonstrate that the new formula leads to better performance compared to other known numerical approximations with similar resolution.

Bäcklund transformation, analytic soliton solutions and numerical simulation for a (2+1)-dimensional complex Ginzburg-Landau equation in a nonlinear fiber

NASA Astrophysics Data System (ADS)

Yu, Ming-Xiao; Tian, Bo; Chai, Jun; Yin, Hui-Min; Du, Zhong

2017-10-01

In this paper, we investigate a nonlinear fiber described by a (2+1)-dimensional complex Ginzburg-Landau equation with the chromatic dispersion, optical filtering, nonlinear and linear gain. Bäcklund transformation in the bilinear form is constructed. With the modified bilinear method, analytic soliton solutions are obtained. For the soliton, the amplitude can decrease or increase when the absolute value of the nonlinear or linear gain is enlarged, and the width can be compressed or amplified when the absolute value of the chromatic dispersion or optical filtering is enhanced. We study the stability of the numerical solutions numerically by applying the increasing amplitude, embedding the white noise and adding the Gaussian pulse to the initial values based on the analytic solutions, which shows that the numerical solutions are stable, not influenced by the finite initial perturbations.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Larsen, E.W.

A class of Projected Discrete-Ordinates (PDO) methods is described for obtaining iterative solutions of discrete-ordinates problems with convergence rates comparable to those observed using Diffusion Synthetic Acceleration (DSA). The spatially discretized PDO solutions are generally not equal to the DSA solutions, but unlike DSA, which requires great care in the use of spatial discretizations to preserve stability, the PDO solutions remain stable and rapidly convergent with essentially arbitrary spatial discretizations. Numerical results are presented which illustrate the rapid convergence and the accuracy of solutions obtained using PDO methods with commonplace differencing methods.

Numerical calculations of two dimensional, unsteady transonic flows with circulation

NASA Technical Reports Server (NTRS)

Beam, R. M.; Warming, R. F.

1974-01-01

The feasibility of obtaining two-dimensional, unsteady transonic aerodynamic data by numerically integrating the Euler equations is investigated. An explicit, third-order-accurate, noncentered, finite-difference scheme is used to compute unsteady flows about airfoils. Solutions for lifting and nonlifting airfoils are presented and compared with subsonic linear theory. The applicability and efficiency of the numerical indicial function method are outlined. Numerically computed subsonic and transonic oscillatory aerodynamic coefficients are presented and compared with those obtained from subsonic linear theory and transonic wind-tunnel data.

Kinetic theory analysis of rarefied gas flow through finite length slots

NASA Technical Reports Server (NTRS)

Raghuraman, P.

1972-01-01

An analytic study is made of the flow a rarefied monatomic gas through a two dimensional slot. The parameters of the problem are the ratios of downstream to upstream pressures, the Knudsen number at the high pressure end (based on slot half width) and the length to slot half width ratio. A moment method of solution is used by assuming a discontinuous distribution function consisting of four Maxwellians split equally in angular space. Numerical solutions are obtained for the resulting equations. The characteristics of the transition regime are portrayed. The solutions in the free molecule limit are systematically lower than the results obtained in that limit by more accurate numerical methods.

Spectral methods in general relativity and large Randall-Sundrum II black holes

NASA Astrophysics Data System (ADS)

Abdolrahimi, Shohreh; Cattoën, Céline; Page, Don N.; \\\\; Yaghoobpour-Tari, Shima

2013-06-01

Using a novel numerical spectral method, we have found solutions for large static Randall-Sundrum II (RSII) black holes by perturbing a numerical AdS5-CFT4 solution to the Einstein equation with a negative cosmological constant Λ that is asymptotically conformal to the Schwarzschild metric. We used a numerical spectral method independent of the Ricci-DeTurck-flow method used by Figueras, Lucietti, and Wiseman for a similar numerical solution. We have compared our black-hole solution to the one Figueras and Wiseman have derived by perturbing their numerical AdS5-CFT4 solution, showing that our solution agrees closely with theirs. We have obtained a closed-form approximation to the metric of the black hole on the brane. We have also deduced the new results that to first order in 1/(-ΛM2), the Hawking temperature and entropy of an RSII static black hole have the same values as the Schwarzschild metric with the same mass, but the horizon area is increased by about 4.7/(-Λ).

Eulerian Lagrangian Adaptive Fup Collocation Method for solving the conservative solute transport in heterogeneous porous media

NASA Astrophysics Data System (ADS)

Gotovac, Hrvoje; Srzic, Veljko

2014-05-01

Contaminant transport in natural aquifers is a complex, multiscale process that is frequently studied using different Eulerian, Lagrangian and hybrid numerical methods. Conservative solute transport is typically modeled using the advection-dispersion equation (ADE). Despite the large number of available numerical methods that have been developed to solve it, the accurate numerical solution of the ADE still presents formidable challenges. In particular, current numerical solutions of multidimensional advection-dominated transport in non-uniform velocity fields are affected by one or all of the following problems: numerical dispersion that introduces artificial mixing and dilution, grid orientation effects, unresolved spatial and temporal scales and unphysical numerical oscillations (e.g., Herrera et al, 2009; Bosso et al., 2012). In this work we will present Eulerian Lagrangian Adaptive Fup Collocation Method (ELAFCM) based on Fup basis functions and collocation approach for spatial approximation and explicit stabilized Runge-Kutta-Chebyshev temporal integration (public domain routine SERK2) which is especially well suited for stiff parabolic problems. Spatial adaptive strategy is based on Fup basis functions which are closely related to the wavelets and splines so that they are also compactly supported basis functions; they exactly describe algebraic polynomials and enable a multiresolution adaptive analysis (MRA). MRA is here performed via Fup Collocation Transform (FCT) so that at each time step concentration solution is decomposed using only a few significant Fup basis functions on adaptive collocation grid with appropriate scales (frequencies) and locations, a desired level of accuracy and a near minimum computational cost. FCT adds more collocations points and higher resolution levels only in sensitive zones with sharp concentration gradients, fronts and/or narrow transition zones. According to the our recent achievements there is no need for solving the large linear system on adaptive grid because each Fup coefficient is obtained by predefined formulas equalizing Fup expansion around corresponding collocation point and particular collocation operator based on few surrounding solution values. Furthermore, each Fup coefficient can be obtained independently which is perfectly suited for parallel processing. Adaptive grid in each time step is obtained from solution of the last time step or initial conditions and advective Lagrangian step in the current time step according to the velocity field and continuous streamlines. On the other side, we implement explicit stabilized routine SERK2 for dispersive Eulerian part of solution in the current time step on obtained spatial adaptive grid. Overall adaptive concept does not require the solving of large linear systems for the spatial and temporal approximation of conservative transport. Also, this new Eulerian-Lagrangian-Collocation scheme resolves all mentioned numerical problems due to its adaptive nature and ability to control numerical errors in space and time. Proposed method solves advection in Lagrangian way eliminating problems in Eulerian methods, while optimal collocation grid efficiently describes solution and boundary conditions eliminating usage of large number of particles and other problems in Lagrangian methods. Finally, numerical tests show that this approach enables not only accurate velocity field, but also conservative transport even in highly heterogeneous porous media resolving all spatial and temporal scales of concentration field.

Application of matched asymptotic expansions to lunar and interplanetary trajectories. Volume 1: Technical discussion

NASA Technical Reports Server (NTRS)

Lancaster, J. E.

1973-01-01

Previously published asymptotic solutions for lunar and interplanetary trajectories have been modified and combined to formulate a general analytical solution to the problem on N-bodies. The earlier first-order solutions, derived by the method of matched asymptotic expansions, have been extended to second order for the purpose of obtaining increased accuracy. The derivation of the second-order solution is summarized by showing the essential steps, some in functional form. The general asymptotic solution has been used as a basis for formulating a number of analytical two-point boundary value solutions. These include earth-to-moon, one- and two-impulse moon-to-earth, and interplanetary solutions. The results show that the accuracies of the asymptotic solutions range from an order of magnitude better than conic approximations to that of numerical integration itself. Also, since no iterations are required, the asymptotic boundary value solutions are obtained in a fraction of the time required for comparable numerically integrated solutions. The subject of minimizing the second-order error is discussed, and recommendations made for further work directed toward achieving a uniform accuracy in all applications.

Formulation of a General Technique for Predicting Pneumatic Attenuation Errors in Airborne Pressure Sensing Devices

NASA Technical Reports Server (NTRS)

Whitmore, Stephen A.

1988-01-01

Presented is a mathematical model derived from the Navier-Stokes equations of momentum and continuity, which may be accurately used to predict the behavior of conventionally mounted pneumatic sensing systems subject to arbitrary pressure inputs. Numerical techniques for solving the general model are developed. Both step and frequency response lab tests were performed. These data are compared with solutions of the mathematical model and show excellent agreement. The procedures used to obtain the lab data are described. In-flight step and frequency response data were obtained. Comparisons with numerical solutions of the math model show good agreement. Procedures used to obtain the flight data are described. Difficulties encountered with obtaining the flight data are discussed.

Upwind schemes and bifurcating solutions in real gas computations

NASA Technical Reports Server (NTRS)

Suresh, Ambady; Liou, Meng-Sing

1992-01-01

The area of high speed flow is seeing a renewed interest due to advanced propulsion concepts such as the National Aerospace Plane (NASP), Space Shuttle, and future civil transport concepts. Upwind schemes to solve such flows have become increasingly popular in the last decade due to their excellent shock capturing properties. In the first part of this paper the authors present the extension of the Osher scheme to equilibrium and non-equilibrium gases. For simplicity, the source terms are treated explicitly. Computations based on the above scheme are presented to demonstrate the feasibility, accuracy and efficiency of the proposed scheme. One of the test problems is a Chapman-Jouguet detonation problem for which numerical solutions have been known to bifurcate into spurious weak detonation solutions on coarse grids. Results indicate that the numerical solution obtained depends both on the upwinding scheme used and the limiter employed to obtain second order accuracy. For example, the Osher scheme gives the correct CJ solution when the super-bee limiter is used, but gives the spurious solution when the Van Leer limiter is used. With the Roe scheme the spurious solution is obtained for all limiters.

The dynamics of a delayed predator-prey model with state dependent feedback control

DOE Office of Scientific and Technical Information (OSTI.GOV)

Singh, Anuraj; Gakkhar, Sunita

2011-11-30

A delayed prey-predator model with state-dependent impulses is investigated. The sufficient conditions of existence and stability of semi-trivial solution and positive period-1 solution are obtained by using the Poincare map and analogue of the Poincare Criterion. The qualitative analysis shows that the positive period-one solution bifurcates from the semi-trivial solution through a fold bifurcation. The complex dynamics including chaos is obtained and numerical simulations substantiate the analytical results.

Introduction to Numerical Methods

DOE Office of Scientific and Technical Information (OSTI.GOV)

Schoonover, Joseph A.

2016-06-14

These are slides for a lecture for the Parallel Computing Summer Research Internship at the National Security Education Center. This gives an introduction to numerical methods. Repetitive algorithms are used to obtain approximate solutions to mathematical problems, using sorting, searching, root finding, optimization, interpolation, extrapolation, least squares regresion, Eigenvalue problems, ordinary differential equations, and partial differential equations. Many equations are shown. Discretizations allow us to approximate solutions to mathematical models of physical systems using a repetitive algorithm and introduce errors that can lead to numerical instabilities if we are not careful.

Solutions of conformal Israel-Stewart relativistic viscous fluid dynamics

NASA Astrophysics Data System (ADS)

Marrochio, Hugo; Noronha, Jorge; Denicol, Gabriel S.; Luzum, Matthew; Jeon, Sangyong; Gale, Charles

2015-01-01

We use symmetry arguments developed by Gubser to construct the first radially expanding explicit solutions of the Israel-Stewart formulation of hydrodynamics. Along with a general semi-analytical solution, an exact analytical solution is given which is valid in the cold plasma limit where viscous effects from shear viscosity and the relaxation time coefficient are important. The radially expanding solutions presented in this paper can be used as nontrivial checks of numerical algorithms employed in hydrodynamic simulations of the quark-gluon plasma formed in ultrarelativistic heavy ion collisions. We show this explicitly by comparing such analytic and semi-analytic solutions with the corresponding numerical solutions obtained using the music viscous hydrodynamics simulation code.

Real gas flow fields about three dimensional configurations

NASA Technical Reports Server (NTRS)

Balakrishnan, A.; Lombard, C. K.; Davy, W. C.

1983-01-01

Real gas, inviscid supersonic flow fields over a three-dimensional configuration are determined using a factored implicit algorithm. Air in chemical equilibrium is considered and its local thermodynamic properties are computed by an equilibrium composition method. Numerical solutions are presented for both real and ideal gases at three different Mach numbers and at two different altitudes. Selected results are illustrated by contour plots and are also tabulated for future reference. Results obtained compare well with existing tabulated numerical solutions and hence validate the solution technique.

Numerical analysis of soliton solutions of the modified Korteweg-de Vries-sine-Gordon equation

NASA Astrophysics Data System (ADS)

Popov, S. P.

2015-03-01

Multisoliton solutions of the modified Korteweg-de Vries-sine-Gordon equation (mKdV-SG) are found numerically by applying the quasi-spectral Fourier method and the fourth-order Runge-Kutta method. The accuracy and features of the approach are determined as applied to problems with initial data in the form of various combinations of perturbed soliton distributions. Three-soliton solutions are obtained, and the generation of kinks, breathers, wobblers, perturbed kinks, and nonlinear oscillatory waves is studied.

A family of four stages embedded explicit six-step methods with eliminated phase-lag and its derivatives for the numerical solution of the second order problems

NASA Astrophysics Data System (ADS)

Simos, T. E.

2017-11-01

A family of four stages high algebraic order embedded explicit six-step methods, for the numerical solution of second order initial or boundary-value problems with periodical and/or oscillating solutions, are studied in this paper. The free parameters of the new proposed methods are calculated solving the linear system of equations which is produced by requesting the vanishing of the phase-lag of the methods and the vanishing of the phase-lag's derivatives of the schemes. For the new obtained methods we investigate: • Its local truncation error (LTE) of the methods.• The asymptotic form of the LTE obtained using as model problem the radial Schrödinger equation.• The comparison of the asymptotic forms of LTEs for several methods of the same family. This comparison leads to conclusions on the efficiency of each method of the family.• The stability and the interval of periodicity of the obtained methods of the new family of embedded finite difference pairs.• The applications of the new obtained family of embedded finite difference pairs to the numerical solution of several second order problems like the radial Schrödinger equation, astronomical problems etc. The above applications lead to conclusion on the efficiency of the methods of the new family of embedded finite difference pairs.

Numerical method for solving the nonlinear four-point boundary value problems

NASA Astrophysics Data System (ADS)

Lin, Yingzhen; Lin, Jinnan

2010-12-01

In this paper, a new reproducing kernel space is constructed skillfully in order to solve a class of nonlinear four-point boundary value problems. The exact solution of the linear problem can be expressed in the form of series and the approximate solution of the nonlinear problem is given by the iterative formula. Compared with known investigations, the advantages of our method are that the representation of exact solution is obtained in a new reproducing kernel Hilbert space and accuracy of numerical computation is higher. Meanwhile we present the convergent theorem, complexity analysis and error estimation. The performance of the new method is illustrated with several numerical examples.

Refined numerical solution of the transonic flow past a wedge

NASA Technical Reports Server (NTRS)

Liang, S.-M.; Fung, K.-Y.

1985-01-01

A numerical procedure combining the ideas of solving a modified difference equation and of adaptive mesh refinement is introduced. The numerical solution on a fixed grid is improved by using better approximations of the truncation error computed from local subdomain grid refinements. This technique is used to obtain refined solutions of steady, inviscid, transonic flow past a wedge. The effects of truncation error on the pressure distribution, wave drag, sonic line, and shock position are investigated. By comparing the pressure drag on the wedge and wave drag due to the shocks, a supersonic-to-supersonic shock originating from the wedge shoulder is confirmed.

Thermal Analysis of Antenna Structures. Part 2: Panel Temperature Distribution

NASA Technical Reports Server (NTRS)

Schonfeld, D.; Lansing, F. L.

1983-01-01

This article is the second in a series that analyzes the temperature distribution in microwave antennas. An analytical solution in a series form is obtained for the temperature distribution in a flat plate analogous to an antenna surface panel under arbitrary temperature and boundary conditions. The solution includes the effects of radiation and air convection from the plate. Good agreement is obtained between the numerical and analytical solutions.

The numerical calculation of laminar boundary-layer separation

NASA Technical Reports Server (NTRS)

Klineberg, J. M.; Steger, J. L.

1974-01-01

Iterative finite-difference techniques are developed for integrating the boundary-layer equations, without approximation, through a region of reversed flow. The numerical procedures are used to calculate incompressible laminar separated flows and to investigate the conditions for regular behavior at the point of separation. Regular flows are shown to be characterized by an integrable saddle-type singularity that makes it difficult to obtain numerical solutions which pass continuously into the separated region. The singularity is removed and continuous solutions ensured by specifying the wall shear distribution and computing the pressure gradient as part of the solution. Calculated results are presented for several separated flows and the accuracy of the method is verified. A computer program listing and complete solution case are included.

An exploratory study of a finite difference method for calculating unsteady transonic potential flow

NASA Technical Reports Server (NTRS)

Bennett, R. M.; Bland, S. R.

1979-01-01

A method for calculating transonic flow over steady and oscillating airfoils was developed by Isogai. The full potential equation is solved with a semi-implicit, time-marching, finite difference technique. Steady flow solutions are obtained from time asymptotic solutions for a steady airfoil. Corresponding oscillatory solutions are obtained by initiating an oscillation and marching in time for several cycles until a converged periodic solution is achieved. The method is described in general terms and results for the case of an airfoil with an oscillating flap are presented for Mach numbers 0.500 and 0.875. Although satisfactory results are obtained for some reduced frequencies, it is found that the numerical technique generates spurious oscillations in the indicial response functions and in the variation of the aerodynamic coefficients with reduced frequency. These oscillations are examined with a dynamic data reduction method to evaluate their effects and trends with reduced frequency and Mach number. Further development of the numerical method is needed to eliminate these oscillations.

A numerical and experimental study of three-dimensional liquid sloshing in a rotating spherical container

NASA Technical Reports Server (NTRS)

Chen, Kuo-Huey; Kelecy, Franklyn J.; Pletcher, Richard H.

1992-01-01

A numerical and experimental study of three dimensional liquid sloshing inside a partially-filled spherical container undergoing an orbital rotating motion is described. Solutions of the unsteady, three-dimensional Navier-Stokes equations for the case of a gradual spin-up from rest are compared with experimental data obtained using a rotating test rig fitted with two liquid-filled spherical tanks. Data gathered from several experiments are reduced in terms of a dimensionless free surface height for comparison with transient results from the numerical simulations. The numerical solutions are found to compare favorably with the experimental data.

A new theoretical basis for numerical simulations of nonlinear acoustic fields

NASA Astrophysics Data System (ADS)

Wójcik, Janusz

2000-07-01

Nonlinear acoustic equations can be considerably simplified. The presented model retains the accuracy of a more complex description of nonlinearity and a uniform description of near and far fields (in contrast to the KZK equation). A method has been presented for obtaining solutions of Kuznetsov's equation from the solutions of the model under consideration. Results of numerical calculations, including comparative ones, are presented.

Existence and numerical simulation of periodic traveling wave solutions to the Casimir equation for the Ito system

NASA Astrophysics Data System (ADS)

Abbasbandy, S.; Van Gorder, R. A.; Hajiketabi, M.; Mesrizadeh, M.

2015-10-01

We consider traveling wave solutions to the Casimir equation for the Ito system (a two-field extension of the KdV equation). These traveling waves are governed by a nonlinear initial value problem with an interesting nonlinearity (which actually amplifies in magnitude as the size of the solution becomes small). The nonlinear problem is parameterized by two initial constant values, and we demonstrate that the existence of solutions is strongly tied to these parameter values. For our interests, we are concerned with positive, bounded, periodic wave solutions. We are able to classify parameter regimes which admit such solutions in full generality, thereby obtaining a nice existence result. Using the existence result, we are then able to numerically simulate the positive, bounded, periodic solutions. We elect to employ a group preserving scheme in order to numerically study these solutions, and an outline of this approach is provided. The numerical simulations serve to illustrate the properties of these solutions predicted analytically through the existence result. Physically, these results demonstrate the existence of a type of space-periodic structure in the Casimir equation for the Ito model, which propagates as a traveling wave.

FOCUSING OF HIGH POWER ULTRASOUND BEAMS AND LIMITING VALUES OF SHOCK WAVE PARAMETERS

PubMed Central

Bessonova, O.V.; Khokhlova, V.A.; Bailey, M.R.; Canney, M.S.; Crum, L.A.

2009-01-01

In this work, the influence of nonlinear and diffraction effects on amplification factors of focused ultrasound systems is investigated. The limiting values of acoustic field parameters obtained by focusing of high power ultrasound are studied. The Khokhlov-Zabolotskaya-Kuznetsov (KZK) equation was used for the numerical modeling. Solutions for the nonlinear acoustic field were obtained at output levels corresponding to both pre- and post- shock formation conditions in the focal area of the beam in a weakly dissipative medium. Numerical solutions were compared with experimental data as well as with known analytic predictions. PMID:20161349

FOCUSING OF HIGH POWER ULTRASOUND BEAMS AND LIMITING VALUES OF SHOCK WAVE PARAMETERS.

PubMed

Bessonova, O V; Khokhlova, V A; Bailey, M R; Canney, M S; Crum, L A

2009-07-21

In this work, the influence of nonlinear and diffraction effects on amplification factors of focused ultrasound systems is investigated. The limiting values of acoustic field parameters obtained by focusing of high power ultrasound are studied. The Khokhlov-Zabolotskaya-Kuznetsov (KZK) equation was used for the numerical modeling. Solutions for the nonlinear acoustic field were obtained at output levels corresponding to both pre- and post- shock formation conditions in the focal area of the beam in a weakly dissipative medium. Numerical solutions were compared with experimental data as well as with known analytic predictions.

Focusing of high power ultrasound beams and limiting values of shock wave parameters

NASA Astrophysics Data System (ADS)

Bessonova, O. V.; Khokhlova, V. A.; Bailey, M. R.; Canney, M. S.; Crum, L. A.

2009-10-01

In this work, the influence of nonlinear and diffraction effects on amplification factors of focused ultrasound systems is investigated. The limiting values of acoustic field parameters obtained by focusing of high power ultrasound are studied. The Khokhlov-Zabolotskaya-Kuznetsov (KZK) equation was used for the numerical modeling. Solutions for the nonlinear acoustic field were obtained at output levels corresponding to both pre- and post-shock formation conditions in the focal area of the beam in a weakly dissipative medium. Numerical solutions were compared with experimental data as well as with known analytic predictions.

Numerical solution of modified differential equations based on symmetry preservation.

PubMed

Ozbenli, Ersin; Vedula, Prakash

2017-12-01

In this paper, we propose a method to construct invariant finite-difference schemes for solution of partial differential equations (PDEs) via consideration of modified forms of the underlying PDEs. The invariant schemes, which preserve Lie symmetries, are obtained based on the method of equivariant moving frames. While it is often difficult to construct invariant numerical schemes for PDEs due to complicated symmetry groups associated with cumbersome discrete variable transformations, we note that symmetries associated with more convenient transformations can often be obtained by appropriately modifying the original PDEs. In some cases, modifications to the original PDEs are also found to be useful in order to avoid trivial solutions that might arise from particular selections of moving frames. In our proposed method, modified forms of PDEs can be obtained either by addition of perturbation terms to the original PDEs or through defect correction procedures. These additional terms, whose primary purpose is to enable symmetries with more convenient transformations, are then removed from the system by considering moving frames for which these specific terms go to zero. Further, we explore selection of appropriate moving frames that result in improvement in accuracy of invariant numerical schemes based on modified PDEs. The proposed method is tested using the linear advection equation (in one- and two-dimensions) and the inviscid Burgers' equation. Results obtained for these tests cases indicate that numerical schemes derived from the proposed method perform significantly better than existing schemes not only by virtue of improvement in numerical accuracy but also due to preservation of qualitative properties or symmetries of the underlying differential equations.

Jacobi-Gauss-Lobatto collocation method for the numerical solution of 1+1 nonlinear Schrödinger equations

NASA Astrophysics Data System (ADS)

Doha, E. H.; Bhrawy, A. H.; Abdelkawy, M. A.; Van Gorder, Robert A.

2014-03-01

A Jacobi-Gauss-Lobatto collocation (J-GL-C) method, used in combination with the implicit Runge-Kutta method of fourth order, is proposed as a numerical algorithm for the approximation of solutions to nonlinear Schrödinger equations (NLSE) with initial-boundary data in 1+1 dimensions. Our procedure is implemented in two successive steps. In the first one, the J-GL-C is employed for approximating the functional dependence on the spatial variable, using (N-1) nodes of the Jacobi-Gauss-Lobatto interpolation which depends upon two general Jacobi parameters. The resulting equations together with the two-point boundary conditions induce a system of 2(N-1) first-order ordinary differential equations (ODEs) in time. In the second step, the implicit Runge-Kutta method of fourth order is applied to solve this temporal system. The proposed J-GL-C method, used in combination with the implicit Runge-Kutta method of fourth order, is employed to obtain highly accurate numerical approximations to four types of NLSE, including the attractive and repulsive NLSE and a Gross-Pitaevskii equation with space-periodic potential. The numerical results obtained by this algorithm have been compared with various exact solutions in order to demonstrate the accuracy and efficiency of the proposed method. Indeed, for relatively few nodes used, the absolute error in our numerical solutions is sufficiently small.

Some Finite Difference Solutions of the Laminar Compressible Boundary Layer Showing the Effects of Upstream Transpiration Cooling

NASA Technical Reports Server (NTRS)

Howe, John T.

1959-01-01

Three numerical solutions of the partial differential equations describing the compressible laminar boundary layer are obtained by the finite difference method described in reports by I. Flugge-Lotz, D.C. Baxter, and this author. The solutions apply to steady-state supersonic flow without pressure gradient, over a cold wall and over an adiabatic wall, both having transpiration cooling upstream, and over an adiabatic wall with upstream cooling but without upstream transpiration. It is shown that for a given upstream wall temperature, upstream transpiration cooling affords much better protection to the adiabatic solid wall than does upstream cooling without transpiration. The results of the numerical solutions are compared with those of approximate solutions. The thermal results of the finite difference solution lie between the results of Rubesin and Inouye, and those of Libby and Pallone. When the skin-friction results of one finite difference solution are used in the thermal analysis of Rubesin and Inouye, improved agreement between the thermal results of the two methods of solution is obtained.

Time-Accurate Solutions of Incompressible Navier-Stokes Equations for Potential Turbopump Applications

NASA Technical Reports Server (NTRS)

Kiris, Cetin; Kwak, Dochan

2001-01-01

Two numerical procedures, one based on artificial compressibility method and the other pressure projection method, are outlined for obtaining time-accurate solutions of the incompressible Navier-Stokes equations. The performance of the two method are compared by obtaining unsteady solutions for the evolution of twin vortices behind a at plate. Calculated results are compared with experimental and other numerical results. For an un- steady ow which requires small physical time step, pressure projection method was found to be computationally efficient since it does not require any subiterations procedure. It was observed that the artificial compressibility method requires a fast convergence scheme at each physical time step in order to satisfy incompressibility condition. This was obtained by using a GMRES-ILU(0) solver in our computations. When a line-relaxation scheme was used, the time accuracy was degraded and time-accurate computations became very expensive.

Study on a Multi-Frequency Homotopy Analysis Method for Period-Doubling Solutions of Nonlinear Systems

NASA Astrophysics Data System (ADS)

Fu, H. X.; Qian, Y. H.

In this paper, a modification of homotopy analysis method (HAM) is applied to study the two-degree-of-freedom coupled Duffing system. Firstly, the process of calculating the two-degree-of-freedom coupled Duffing system is presented. Secondly, the single periodic solutions and double periodic solutions are obtained by solving the constructed nonlinear algebraic equations. Finally, comparing the periodic solutions obtained by the multi-frequency homotopy analysis method (MFHAM) and the fourth-order Runge-Kutta method, it is found that the approximate solution agrees well with the numerical solution.

Assessment of a high-resolution central scheme for the solution of the relativistic hydrodynamics equations

NASA Astrophysics Data System (ADS)

Lucas-Serrano, A.; Font, J. A.; Ibáñez, J. M.; Martí, J. M.

2004-12-01

We assess the suitability of a recent high-resolution central scheme developed by \\cite{kurganov} for the solution of the relativistic hydrodynamic equations. The novelty of this approach relies on the absence of Riemann solvers in the solution procedure. The computations we present are performed in one and two spatial dimensions in Minkowski spacetime. Standard numerical experiments such as shock tubes and the relativistic flat-faced step test are performed. As an astrophysical application the article includes two-dimensional simulations of the propagation of relativistic jets using both Cartesian and cylindrical coordinates. The simulations reported clearly show the capabilities of the numerical scheme of yielding satisfactory results, with an accuracy comparable to that obtained by the so-called high-resolution shock-capturing schemes based upon Riemann solvers (Godunov-type schemes), even well inside the ultrarelativistic regime. Such a central scheme can be straightforwardly applied to hyperbolic systems of conservation laws for which the characteristic structure is not explicitly known, or in cases where a numerical computation of the exact solution of the Riemann problem is prohibitively expensive. Finally, we present comparisons with results obtained using various Godunov-type schemes as well as with those obtained using other high-resolution central schemes which have recently been reported in the literature.

A two-dimensional numerical study of the flow inside the combustion chambers of a motored rotary engine

NASA Technical Reports Server (NTRS)

Shih, T. I. P.; Yang, S. L.; Schock, H. J.

1986-01-01

A numerical study was performed to investigate the unsteady, multidimensional flow inside the combustion chambers of an idealized, two-dimensional, rotary engine under motored conditions. The numerical study was based on the time-dependent, two-dimensional, density-weighted, ensemble-averaged conservation equations of mass, species, momentum, and total energy valid for two-component ideal gas mixtures. The ensemble-averaged conservation equations were closed by a K-epsilon model of turbulence. This K-epsilon model of turbulence was modified to account for some of the effects of compressibility, streamline curvature, low-Reynolds number, and preferential stress dissipation. Numerical solutions to the conservation equations were obtained by the highly efficient implicit-factored method of Beam and Warming. The grid system needed to obtain solutions were generated by an algebraic grid generation technique based on transfinite interpolation. Results of the numerical study are presented in graphical form illustrating the flow patterns during intake, compression, gaseous fuel injection, expansion, and exhaust.

A two-dimensional numerical study of the flow inside the combustion chamber of a motored rotary engine

NASA Technical Reports Server (NTRS)

Shih, T. I-P.; Yang, S. L.; Schock, H. J.

1986-01-01

A numerical study was performed to investigate the unsteady, multidimensional flow inside the combustion chambers of an idealized, two-dimensional, rotary engine under motored conditions. The numerical study was based on the time-dependent, two-dimensional, density-weighted, ensemble-averaged conservation equations of mass, species, momentum, and total energy valid for two-component ideal gas mixtures. The ensemble-averaged conservation equations were closed by a K-epsilon model of turbulence. This K-epsilon model of turbulence was modified to account for some of the effects of compressibility, streamline curvature, low-Reynolds number, and preferential stress dissipation. Numerical solutions to the conservation equations were obtained by the highly efficient implicit-factored method of Beam and Warming. The grid system needed to obtain solutions were generated by an algebraic grid generation technique based on transfinite interpolation. Results of the numerical study are presented in graphical form illustrating the flow patterns during intake, compression, gaseous fuel injection, expansion, and exhaust.

Transonic Navier-Stokes solutions of three-dimensional afterbody flows

NASA Technical Reports Server (NTRS)

Compton, William B., III; Thomas, James L.; Abeyounis, William K.; Mason, Mary L.

1989-01-01

The performance of a three-dimensional Navier-Stokes solution technique in predicting the transonic flow past a nonaxisymmetric nozzle was investigated. The investigation was conducted at free-stream Mach numbers ranging from 0.60 to 0.94 and an angle of attack of 0 degrees. The numerical solution procedure employs the three-dimensional, unsteady, Reynolds-averaged Navier-Stokes equations written in strong conservation form, a thin layer assumption, and the Baldwin-Lomax turbulence model. The equations are solved by using the finite-volume principle in conjunction with an approximately factored upwind-biased numerical algorithm. In the numerical procedure, the jet exhaust is represented by a solid sting. Wind-tunnel data with the jet exhaust simulated by high pressure air were also obtained to compare with the numerical calculations.

Numerical solution of the exterior oblique derivative BVP using the direct BEM formulation

NASA Astrophysics Data System (ADS)

Čunderlík, Róbert; Špir, Róbert; Mikula, Karol

2016-04-01

The fixed gravimetric boundary value problem (FGBVP) represents an exterior oblique derivative problem for the Laplace equation. A direct formulation of the boundary element method (BEM) for the Laplace equation leads to a boundary integral equation (BIE) where a harmonic function is represented as a superposition of the single-layer and double-layer potential. Such a potential representation is applied to obtain a numerical solution of FGBVP. The oblique derivative problem is treated by a decomposition of the gradient of the unknown disturbing potential into its normal and tangential components. Our numerical scheme uses the collocation with linear basis functions. It involves a triangulated discretization of the Earth's surface as our computational domain considering its complicated topography. To achieve high-resolution numerical solutions, parallel implementations using the MPI subroutines as well as an iterative elimination of far zones' contributions are performed. Numerical experiments present a reconstruction of a harmonic function above the Earth's topography given by the spherical harmonic approach, namely by the EGM2008 geopotential model up to degree 2160. The SRTM30 global topography model is used to approximate the Earth's surface by the triangulated discretization. The obtained BEM solution with the resolution 0.05 deg (12,960,002 nodes) is compared with EGM2008. The standard deviation of residuals 5.6 cm indicates a good agreement. The largest residuals are obviously in high mountainous regions. They are negative reaching up to -0.7 m in Himalayas and about -0.3 m in Andes and Rocky Mountains. A local refinement in the area of Slovakia confirms an improvement of the numerical solution in this mountainous region despite of the fact that the Earth's topography is here considered in more details.

Charged Analogues of Henning Knutsen Type Solutions in General Relativity

NASA Astrophysics Data System (ADS)

Gupta, Y. K.; Kumar, Sachin; Pratibha

2011-11-01

In the present article, we have found charged analogues of Henning Knutsen's interior solutions which join smoothly to the Reissner-Nordstrom metric at the pressure free interface. The solutions are singularity free and analyzed numerically with respect to pressure, energy-density and charge-density in details. The solutions so obtained also present the generalization of A.L. Mehra's solutions.

Numerical simulation of the transonic flow past the blunted wedge in the diverging channel

NASA Astrophysics Data System (ADS)

Ryabinin, Anatoly

2018-05-01

Positions of shock waves in the 2D channel with a blunted wedge are studied numerically. Solutions of the Euler equations are obtained with finite-volume solver SU2 for 15 variants of channel geometry. Numerical simulations reveal a considerable hysteresis in the shock wave position versus the supersonic Mach number given at the inlet. In the certain range of inlet Mach number, there are asymmetrical solutions of the equations. Small change in the geometry of the channel leads to shift of boundaries of the hysteresis range.

A multistage time-stepping scheme for the Navier-Stokes equations

NASA Technical Reports Server (NTRS)

Swanson, R. C.; Turkel, E.

1985-01-01

A class of explicit multistage time-stepping schemes is used to construct an algorithm for solving the compressible Navier-Stokes equations. Flexibility in treating arbitrary geometries is obtained with a finite-volume formulation. Numerical efficiency is achieved by employing techniques for accelerating convergence to steady state. Computer processing is enhanced through vectorization of the algorithm. The scheme is evaluated by solving laminar and turbulent flows over a flat plate and an NACA 0012 airfoil. Numerical results are compared with theoretical solutions or other numerical solutions and/or experimental data.

Small-x asymptotics of the quark helicity distribution: Analytic results

DOE PAGES

Kovchegov, Yuri V.; Pitonyak, Daniel; Sievert, Matthew D.

2017-06-15

In this Letter, we analytically solve the evolution equations for the small-x asymptotic behavior of the (flavor singlet) quark helicity distribution in the large- N c limit. Here, these evolution equations form a set of coupled integro-differential equations, which previously could only be solved numerically. This approximate numerical solution, however, revealed simplifying properties of the small-x asymptotics, which we exploit here to obtain an analytic solution.

Numerical solution of the incompressible Navier-Stokes equations. Ph.D. Thesis - Stanford Univ., Mar. 1989

NASA Technical Reports Server (NTRS)

Rogers, Stuart E.

1990-01-01

The current work is initiated in an effort to obtain an efficient, accurate, and robust algorithm for the numerical solution of the incompressible Navier-Stokes equations in two- and three-dimensional generalized curvilinear coordinates for both steady-state and time-dependent flow problems. This is accomplished with the use of the method of artificial compressibility and a high-order flux-difference splitting technique for the differencing of the convective terms. Time accuracy is obtained in the numerical solutions by subiterating the equations in psuedo-time for each physical time step. The system of equations is solved with a line-relaxation scheme which allows the use of very large pseudo-time steps leading to fast convergence for steady-state problems as well as for the subiterations of time-dependent problems. Numerous laminar test flow problems are computed and presented with a comparison against analytically known solutions or experimental results. These include the flow in a driven cavity, the flow over a backward-facing step, the steady and unsteady flow over a circular cylinder, flow over an oscillating plate, flow through a one-dimensional inviscid channel with oscillating back pressure, the steady-state flow through a square duct with a 90 degree bend, and the flow through an artificial heart configuration with moving boundaries. An adequate comparison with the analytical or experimental results is obtained in all cases. Numerical comparisons of the upwind differencing with central differencing plus artificial dissipation indicates that the upwind differencing provides a much more robust algorithm, which requires significantly less computing time. The time-dependent problems require on the order of 10 to 20 subiterations, indicating that the elliptical nature of the problem does require a substantial amount of computing effort.

Stable Numerical Approach for Fractional Delay Differential Equations

NASA Astrophysics Data System (ADS)

Singh, Harendra; Pandey, Rajesh K.; Baleanu, D.

2017-12-01

In this paper, we present a new stable numerical approach based on the operational matrix of integration of Jacobi polynomials for solving fractional delay differential equations (FDDEs). The operational matrix approach converts the FDDE into a system of linear equations, and hence the numerical solution is obtained by solving the linear system. The error analysis of the proposed method is also established. Further, a comparative study of the approximate solutions is provided for the test examples of the FDDE by varying the values of the parameters in the Jacobi polynomials. As in special case, the Jacobi polynomials reduce to the well-known polynomials such as (1) Legendre polynomial, (2) Chebyshev polynomial of second kind, (3) Chebyshev polynomial of third and (4) Chebyshev polynomial of fourth kind respectively. Maximum absolute error and root mean square error are calculated for the illustrated examples and presented in form of tables for the comparison purpose. Numerical stability of the presented method with respect to all four kind of polynomials are discussed. Further, the obtained numerical results are compared with some known methods from the literature and it is observed that obtained results from the proposed method is better than these methods.

Graph theory approach to the eigenvalue problem of large space structures

NASA Technical Reports Server (NTRS)

Reddy, A. S. S. R.; Bainum, P. M.

1981-01-01

Graph theory is used to obtain numerical solutions to eigenvalue problems of large space structures (LSS) characterized by a state vector of large dimensions. The LSS are considered as large, flexible systems requiring both orientation and surface shape control. Graphic interpretation of the determinant of a matrix is employed to reduce a higher dimensional matrix into combinations of smaller dimensional sub-matrices. The reduction is implemented by means of a Boolean equivalent of the original matrices formulated to obtain smaller dimensional equivalents of the original numerical matrix. Computation time becomes less and more accurate solutions are possible. An example is provided in the form of a free-free square plate. Linearized system equations and numerical values of a stiffness matrix are presented, featuring a state vector with 16 components.

An efficient technique for the numerical solution of the bidomain equations.

PubMed

Whiteley, Jonathan P

2008-08-01

Computing the numerical solution of the bidomain equations is widely accepted to be a significant computational challenge. In this study we extend a previously published semi-implicit numerical scheme with good stability properties that has been used to solve the bidomain equations (Whiteley, J.P. IEEE Trans. Biomed. Eng. 53:2139-2147, 2006). A new, efficient numerical scheme is developed which utilizes the observation that the only component of the ionic current that must be calculated on a fine spatial mesh and updated frequently is the fast sodium current. Other components of the ionic current may be calculated on a coarser mesh and updated less frequently, and then interpolated onto the finer mesh. Use of this technique to calculate the transmembrane potential and extracellular potential induces very little error in the solution. For the simulations presented in this study an increase in computational efficiency of over two orders of magnitude over standard numerical techniques is obtained.

Hybrid asymptotic-numerical approach for estimating first-passage-time densities of the two-dimensional narrow capture problem.

PubMed

Lindsay, A E; Spoonmore, R T; Tzou, J C

2016-10-01

A hybrid asymptotic-numerical method is presented for obtaining an asymptotic estimate for the full probability distribution of capture times of a random walker by multiple small traps located inside a bounded two-dimensional domain with a reflecting boundary. As motivation for this study, we calculate the variance in the capture time of a random walker by a single interior trap and determine this quantity to be comparable in magnitude to the mean. This implies that the mean is not necessarily reflective of typical capture times and that the full density must be determined. To solve the underlying diffusion equation, the method of Laplace transforms is used to obtain an elliptic problem of modified Helmholtz type. In the limit of vanishing trap sizes, each trap is represented as a Dirac point source that permits the solution of the transform equation to be represented as a superposition of Helmholtz Green's functions. Using this solution, we construct asymptotic short-time solutions of the first-passage-time density, which captures peaks associated with rapid capture by the absorbing traps. When numerical evaluation of the Helmholtz Green's function is employed followed by numerical inversion of the Laplace transform, the method reproduces the density for larger times. We demonstrate the accuracy of our solution technique with a comparison to statistics obtained from a time-dependent solution of the diffusion equation and discrete particle simulations. In particular, we demonstrate that the method is capable of capturing the multimodal behavior in the capture time density that arises when the traps are strategically arranged. The hybrid method presented can be applied to scenarios involving both arbitrary domains and trap shapes.

Nonlinear earthquake analysis of reinforced concrete frames with fiber and Bernoulli-Euler beam-column element.

PubMed

Karaton, Muhammet

2014-01-01

A beam-column element based on the Euler-Bernoulli beam theory is researched for nonlinear dynamic analysis of reinforced concrete (RC) structural element. Stiffness matrix of this element is obtained by using rigidity method. A solution technique that included nonlinear dynamic substructure procedure is developed for dynamic analyses of RC frames. A predicted-corrected form of the Bossak-α method is applied for dynamic integration scheme. A comparison of experimental data of a RC column element with numerical results, obtained from proposed solution technique, is studied for verification the numerical solutions. Furthermore, nonlinear cyclic analysis results of a portal reinforced concrete frame are achieved for comparing the proposed solution technique with Fibre element, based on flexibility method. However, seismic damage analyses of an 8-story RC frame structure with soft-story are investigated for cases of lumped/distributed mass and load. Damage region, propagation, and intensities according to both approaches are researched.

Finite analytic numerical solution of heat transfer and flow past a square channel cavity

NASA Technical Reports Server (NTRS)

Chen, C.-J.; Obasih, K.

1982-01-01

A numerical solution of flow and heat transfer characteristics is obtained by the finite analytic method for a two dimensional laminar channel flow over a two-dimensional square cavity. The finite analytic method utilizes the local analytic solution in a small element of the problem region to form the algebraic equation relating an interior nodal value with its surrounding nodal values. Stable and rapidly converged solutions were obtained for Reynolds numbers ranging to 1000 and Prandtl number to 10. Streamfunction, vorticity and temperature profiles are solved. Local and mean Nusselt number are given. It is found that the separation streamlines between the cavity and channel flow are concave into the cavity at low Reynolds number and convex at high Reynolds number (Re greater than 100) and for square cavity the mean Nusselt number may be approximately correlated with Peclet number as Nu(m) = 0.365 Pe exp 0.2.

Design of two-dimensional channels with prescribed velocity distributions along the channel walls

NASA Technical Reports Server (NTRS)

Stanitz, John D

1953-01-01

A general method of design is developed for two-dimensional unbranched channels with prescribed velocities as a function of arc length along the channel walls. The method is developed for both compressible and incompressible, irrotational, nonviscous flow and applies to the design of elbows, diffusers, nozzles, and so forth. In part I solutions are obtained by relaxation methods; in part II solutions are obtained by a Green's function. Five numerical examples are given in part I including three elbow designs with the same prescribed velocity as a function of arc length along the channel walls but with incompressible, linearized compressible, and compressible flow. One numerical example is presented in part II for an accelerating elbow with linearized compressible flow, and the time required for the solution by a Green's function in part II was considerably less than the time required for the same solution by relaxation methods in part I.

Diffusion Influenced Adsorption Kinetics.

PubMed

Miura, Toshiaki; Seki, Kazuhiko

2015-08-27

When the kinetics of adsorption is influenced by the diffusive flow of solutes, the solute concentration at the surface is influenced by the surface coverage of solutes, which is given by the Langmuir-Hinshelwood adsorption equation. The diffusion equation with the boundary condition given by the Langmuir-Hinshelwood adsorption equation leads to the nonlinear integro-differential equation for the surface coverage. In this paper, we solved the nonlinear integro-differential equation using the Grünwald-Letnikov formula developed to solve fractional kinetics. Guided by the numerical results, analytical expressions for the upper and lower bounds of the exact numerical results were obtained. The upper and lower bounds were close to the exact numerical results in the diffusion- and reaction-controlled limits, respectively. We examined the validity of the two simple analytical expressions obtained in the diffusion-controlled limit. The results were generalized to include the effect of dispersive diffusion. We also investigated the effect of molecular rearrangement of anisotropic molecules on surface coverage.

Large Black Holes in the Randall-Sundrum II Model

NASA Astrophysics Data System (ADS)

Yaghoobpour Tari, Shima

The Einstein equation with a negative cosmological constant ! in the five dimensions for the Randall-Sundrum II model, which includes a black hole, has been solved numerically. We have constructed an AdS5-CFT 4 solution numerically, using a spectral method to minimize the integral of the square of the error of the Einstein equation, with 210 parameters to be determined by optimization. This metric is conformal to the Schwarzschild metric at an AdS5 boundary with an infinite scale factor. So, we consider this solution as an infinite-mass black hole solution. We have rewritten the infinite-mass black hole in the Fefferman-Graham form and obtained the numerical components of the CFT energy-momentum tensor. Using them, we have perturbed the metric to relocate the brane from infinity and obtained a large static black hole solution for the Randall- Sundrum II model. The changes of mass, entropy, temperature and area of the large black hole from the Schwarzschild metric are studied up to the first order for the perturbation parameter 1/(-Λ5M 2). The Hawking temperature and entropy for our large black hole have the same values as the Schwarzschild metric with the same mass, but the horizon area is increased by about 4.7/(-Λ5). Figueras, Lucietti, and Wiseman found an AdS5-CFT4 solution using an independent and different method from us, called the Ricci-DeTurck-flow method. Then, Figueras and Wiseman perturbed this solution in a same way as we have done and obtained the solution for the large black hole in the Randall-Sundrum II model. These two numerical solutions are the first mathematical proofs for having a large black hole in the Randall-Sundrum II. We have compared their results with ours for the CFT energy-momentum tensor components and the perturbed metric. We have shown that the results are closely in agreement, which can be considered as evidence that the solution for the large black hole in the Randall-Sundrum II model exists.

The L sub 1 finite element method for pure convection problems

NASA Technical Reports Server (NTRS)

Jiang, Bo-Nan

1991-01-01

The least squares (L sub 2) finite element method is introduced for 2-D steady state pure convection problems with smooth solutions. It is proven that the L sub 2 method has the same stability estimate as the original equation, i.e., the L sub 2 method has better control of the streamline derivative. Numerical convergence rates are given to show that the L sub 2 method is almost optimal. This L sub 2 method was then used as a framework to develop an iteratively reweighted L sub 2 finite element method to obtain a least absolute residual (L sub 1) solution for problems with discontinuous solutions. This L sub 1 finite element method produces a nonoscillatory, nondiffusive and highly accurate numerical solution that has a sharp discontinuity in one element on both coarse and fine meshes. A robust reweighting strategy was also devised to obtain the L sub 1 solution in a few iterations. A number of examples solved by using triangle and bilinear elements are presented.

Application of fractional derivative with exponential law to bi-fractional-order wave equation with frictional memory kernel

NASA Astrophysics Data System (ADS)

Cuahutenango-Barro, B.; Taneco-Hernández, M. A.; Gómez-Aguilar, J. F.

2017-12-01

Analytical solutions of the wave equation with bi-fractional-order and frictional memory kernel of Mittag-Leffler type are obtained via Caputo-Fabrizio fractional derivative in the Liouville-Caputo sense. Through the method of separation of variables and Laplace transform method we derive closed-form solutions and establish fundamental solutions. Special cases with homogeneous Dirichlet boundary conditions and nonhomogeneous initial conditions, as well as for the external force are considered. Numerical simulations of the special solutions were done and novel behaviors are obtained.

Exact analytic solution for the spin-up maneuver of an axially symmetric spacecraft

NASA Astrophysics Data System (ADS)

Ventura, Jacopo; Romano, Marcello

2014-11-01

The problem of spinning-up an axially symmetric spacecraft subjected to an external torque constant in magnitude and parallel to the symmetry axis is considered. The existing exact analytic solution for an axially symmetric body is applied for the first time to this problem. The proposed solution is valid for any initial conditions of attitude and angular velocity and for any length of time and rotation amplitude. Furthermore, the proposed solution can be numerically evaluated up to any desired level of accuracy. Numerical experiments and comparison with an existing approximated solution and with the integration of the equations of motion are reported in the paper. Finally, a new approximated solution obtained from the exact one is introduced in this paper.

Homotopy perturbation method: a versatile tool to evaluate linear and nonlinear fuzzy Volterra integral equations of the second kind.

PubMed

Narayanamoorthy, S; Sathiyapriya, S P

2016-01-01

In this article, we focus on linear and nonlinear fuzzy Volterra integral equations of the second kind and we propose a numerical scheme using homotopy perturbation method (HPM) to obtain fuzzy approximate solutions to them. To facilitate the benefits of this proposal, an algorithmic form of the HPM is also designed to handle the same. In order to illustrate the potentiality of the approach, two test problems are offered and the obtained numerical results are compared with the existing exact solutions and are depicted in terms of plots to reveal its precision and reliability.

Fracture analysis of a transversely isotropic high temperature superconductor strip based on real fundamental solutions

NASA Astrophysics Data System (ADS)

Gao, Zhiwen; Zhou, Youhe

2015-04-01

Real fundamental solution for fracture problem of transversely isotropic high temperature superconductor (HTS) strip is obtained. The superconductor E-J constitutive law is characterized by the Bean model where the critical current density is independent of the flux density. Fracture analysis is performed by the methods of singular integral equations which are solved numerically by Gauss-Lobatto-Chybeshev (GSL) collocation method. To guarantee a satisfactory accuracy, the convergence behavior of the kernel function is investigated. Numerical results of fracture parameters are obtained and the effects of the geometric characteristics, applied magnetic field and critical current density on the stress intensity factors (SIF) are discussed.

Dispersive shock waves in the Kadomtsev-Petviashvili and two dimensional Benjamin-Ono equations

NASA Astrophysics Data System (ADS)

Ablowitz, Mark J.; Demirci, Ali; Ma, Yi-Ping

2016-10-01

Dispersive shock waves (DSWs) in the Kadomtsev-Petviashvili (KP) equation and two dimensional Benjamin-Ono (2DBO) equation are considered using step like initial data along a parabolic front. Employing a parabolic similarity reduction exactly reduces the study of such DSWs in two space one time (2 + 1) dimensions to finding DSW solutions of (1 + 1) dimensional equations. With this ansatz, the KP and 2DBO equations can be exactly reduced to the cylindrical Korteweg-de Vries (cKdV) and cylindrical Benjamin-Ono (cBO) equations, respectively. Whitham modulation equations which describe DSW evolution in the cKdV and cBO equations are derived and Riemann type variables are introduced. DSWs obtained from the numerical solutions of the corresponding Whitham systems and direct numerical simulations of the cKdV and cBO equations are compared with very good agreement obtained. In turn, DSWs obtained from direct numerical simulations of the KP and 2DBO equations are compared with the cKdV and cBO equations, again with good agreement. It is concluded that the (2 + 1) DSW behavior along self similar parabolic fronts can be effectively described by the DSW solutions of the reduced (1 + 1) dimensional equations.

A three-dimensional, compressible, laminar boundary-layer method for general fuselages. Volume 1: Numerical method

NASA Technical Reports Server (NTRS)

Wie, Yong-Sun

1990-01-01

A procedure for calculating 3-D, compressible laminar boundary layer flow on general fuselage shapes is described. The boundary layer solutions can be obtained in either nonorthogonal 'body oriented' coordinates or orthogonal streamline coordinates. The numerical procedure is 'second order' accurate, efficient and independent of the cross flow velocity direction. Numerical results are presented for several test cases, including a sharp cone, an ellipsoid of revolution, and a general aircraft fuselage at angle of attack. Comparisons are made between numerical results obtained using nonorthogonal curvilinear 'body oriented' coordinates and streamline coordinates.

Einstein gravity with torsion induced by the scalar field

NASA Astrophysics Data System (ADS)

Özçelik, H. T.; Kaya, R.; Hortaçsu, M.

2018-06-01

We couple a conformal scalar field in (2+1) dimensions to Einstein gravity with torsion. The field equations are obtained by a variational principle. We could not solve the Einstein and Cartan equations analytically. These equations are solved numerically with 4th order Runge-Kutta method. From the numerical solution, we make an ansatz for the rotation parameter in the proposed metric, which gives an analytical solution for the scalar field for asymptotic regions.

Comparison of viscous-shock-layer solutions by time-asymptotic and steady-state methods. [flow distribution around a Jupiter entry probe

NASA Technical Reports Server (NTRS)

Gupta, R. N.; Moss, J. N.; Simmonds, A. L.

1982-01-01

Two flow-field codes employing the time- and space-marching numerical techniques were evaluated. Both methods were used to analyze the flow field around a massively blown Jupiter entry probe under perfect-gas conditions. In order to obtain a direct point-by-point comparison, the computations were made by using identical grids and turbulence models. For the same degree of accuracy, the space-marching scheme takes much less time as compared to the time-marching method and would appear to provide accurate results for the problems with nonequilibrium chemistry, free from the effect of local differences in time on the final solution which is inherent in time-marching methods. With the time-marching method, however, the solutions are obtainable for the realistic entry probe shapes with massive or uniform surface blowing rates; whereas, with the space-marching technique, it is difficult to obtain converged solutions for such flow conditions. The choice of the numerical method is, therefore, problem dependent. Both methods give equally good results for the cases where results are compared with experimental data.

A new Jacobi spectral collocation method for solving 1+1 fractional Schrödinger equations and fractional coupled Schrödinger systems

NASA Astrophysics Data System (ADS)

Bhrawy, A. H.; Doha, E. H.; Ezz-Eldien, S. S.; Van Gorder, Robert A.

2014-12-01

The Jacobi spectral collocation method (JSCM) is constructed and used in combination with the operational matrix of fractional derivatives (described in the Caputo sense) for the numerical solution of the time-fractional Schrödinger equation (T-FSE) and the space-fractional Schrödinger equation (S-FSE). The main characteristic behind this approach is that it reduces such problems to those of solving a system of algebraic equations, which greatly simplifies the solution process. In addition, the presented approach is also applied to solve the time-fractional coupled Schrödinger system (T-FCSS). In order to demonstrate the validity and accuracy of the numerical scheme proposed, several numerical examples with their approximate solutions are presented with comparisons between our numerical results and those obtained by other methods.

A block iterative finite element algorithm for numerical solution of the steady-state, compressible Navier-Stokes equations

NASA Technical Reports Server (NTRS)

Cooke, C. H.

1976-01-01

An iterative method for numerically solving the time independent Navier-Stokes equations for viscous compressible flows is presented. The method is based upon partial application of the Gauss-Seidel principle in block form to the systems of nonlinear algebraic equations which arise in construction of finite element (Galerkin) models approximating solutions of fluid dynamic problems. The C deg-cubic element on triangles is employed for function approximation. Computational results for a free shear flow at Re = 1,000 indicate significant achievement of economy in iterative convergence rate over finite element and finite difference models which employ the customary time dependent equations and asymptotic time marching procedure to steady solution. Numerical results are in excellent agreement with those obtained for the same test problem employing time marching finite element and finite difference solution techniques.

Supersonic flow of chemically reacting gas-particle mixtures. Volume 1: A theoretical analysis and development of the numerical solution

NASA Technical Reports Server (NTRS)

Penny, M. M.; Smith, S. D.; Anderson, P. G.; Sulyma, P. R.; Pearson, M. L.

1976-01-01

A numerical solution for chemically reacting supersonic gas-particle flows in rocket nozzles and exhaust plumes was described. The gas-particle flow solution is fully coupled in that the effects of particle drag and heat transfer between the gas and particle phases are treated. Gas and particles exchange momentum via the drag exerted on the gas by the particles. Energy is exchanged between the phases via heat transfer (convection and/or radiation). Thermochemistry calculations (chemical equilibrium, frozen or chemical kinetics) were shown to be uncoupled from the flow solution and, as such, can be solved separately. The solution to the set of governing equations is obtained by utilizing the method of characteristics. The equations cast in characteristic form are shown to be formally the same for ideal, frozen, chemical equilibrium and chemical non-equilibrium reacting gas mixtures. The particle distribution is represented in the numerical solution by a finite distribution of particle sizes.

Theoretical analysis of incompressible flow through a radial-inlet centrifugal impeller at various weight flows

NASA Technical Reports Server (NTRS)

Kramer, James J; Prian, Vasily D; Wu, Chung-Hua

1956-01-01

A method for the solution of the incompressible nonviscous flow through a centrifugal impeller, including the inlet region, is presented. Several numerical solutions are obtained for four weight flows through an impeller at one operating speed. These solutions are refined in the leading-edge region. The results are presented in a series of figures showing streamlines and relative velocity contours. A comparison is made with the results obtained by using a rapid approximate method of analysis.

Numerical simulation of the flow about the F-18 HARV at high angle of attack

NASA Technical Reports Server (NTRS)

Murman, Scott M.

1994-01-01

As part of NASA's High Alpha Technology Program, research has been aimed at developing and extending numerical methods to accurately predict the high Reynolds number flow about the NASA F-18 High Alpha Research Vehicle (HARV) at large angles of attack. The HARV aircraft is equipped with a bidirectional thrust vectoring unit which enables stable, controlled flight through 70 deg angle of attack. Currently, high-fidelity numerical solutions for the flow about the HARV have been obtained at alpha = 30 deg, and validated against flight-test data. It is planned to simulate the flow about the HARV through alpha = 60 deg, and obtain solutions of the same quality as those at the lower angles of attack. This report presents the status of work aimed at extending the HARV computations to the extreme angle of attack range.

Numerical quadrature methods for integrals of singular periodic functions and their application to singular and weakly singular integral equations

NASA Technical Reports Server (NTRS)

Sidi, A.; Israeli, M.

1986-01-01

High accuracy numerical quadrature methods for integrals of singular periodic functions are proposed. These methods are based on the appropriate Euler-Maclaurin expansions of trapezoidal rule approximations and their extrapolations. They are used to obtain accurate quadrature methods for the solution of singular and weakly singular Fredholm integral equations. Such periodic equations are used in the solution of planar elliptic boundary value problems, elasticity, potential theory, conformal mapping, boundary element methods, free surface flows, etc. The use of the quadrature methods is demonstrated with numerical examples.

Numerical method and FORTRAN program for the solution of an axisymmetric electrostatic collector design problem

NASA Technical Reports Server (NTRS)

Reese, O. W.

1972-01-01

The numerical calculation is described of the steady-state flow of electrons in an axisymmetric, spherical, electrostatic collector for a range of boundary conditions. The trajectory equations of motion are solved alternately with Poisson's equation for the potential field until convergence is achieved. A direct (noniterative) numerical technique is used to obtain the solution to Poisson's equation. Space charge effects are included for initial current densities as large as 100 A/sq cm. Ways of dealing successfully with the difficulties associated with these high densities are discussed. A description of the mathematical model, a discussion of numerical techniques, results from two typical runs, and the FORTRAN computer program are included.

Generalized bipartite quantum state discrimination problems with sequential measurements

NASA Astrophysics Data System (ADS)

Nakahira, Kenji; Kato, Kentaro; Usuda, Tsuyoshi Sasaki

2018-02-01

We investigate an optimization problem of finding quantum sequential measurements, which forms a wide class of state discrimination problems with the restriction that only local operations and one-way classical communication are allowed. Sequential measurements from Alice to Bob on a bipartite system are considered. Using the fact that the optimization problem can be formulated as a problem with only Alice's measurement and is convex programming, we derive its dual problem and necessary and sufficient conditions for an optimal solution. Our results are applicable to various practical optimization criteria, including the Bayes criterion, the Neyman-Pearson criterion, and the minimax criterion. In the setting of the problem of finding an optimal global measurement, its dual problem and necessary and sufficient conditions for an optimal solution have been widely used to obtain analytical and numerical expressions for optimal solutions. Similarly, our results are useful to obtain analytical and numerical expressions for optimal sequential measurements. Examples in which our results can be used to obtain an analytical expression for an optimal sequential measurement are provided.

Study of viscous flow about airfoils by the integro-differential method

NASA Technical Reports Server (NTRS)

Wu, J. C.; Sampath, S.

1975-01-01

An integro-differential method was used for numerically solving unsteady incompressible viscous flow problems. A computer program was prepared to solve the problem of an impulsively started 9% thick symmetric Joukowski airfoil at an angle of attack of 15 deg and a Reynolds number of 1000. Some of the results obtained for this problem were discussed and compared with related work completed previously. Two numerical procedures were used, an Alternating Direction Implicit (ADI) method and a Successive Line Relaxation (SLR) method. Generally, the ADI solution agrees well with the SLR solution and with previous results are stations away from the trailing edge. At the trailing edge station, the ADI solution differs substantially from previous results, while the vorticity profiles obtained from the SLR method there are in good qualitative agreement with previous results.

On the anomaly of velocity-pressure decoupling in collocated mesh solutions

NASA Technical Reports Server (NTRS)

Kim, Sang-Wook; Vanoverbeke, Thomas

1991-01-01

The use of various pressure correction algorithms originally developed for fully staggered meshes can yield a velocity-pressure decoupled solution for collocated meshes. The mechanism that causes velocity-pressure decoupling is identified. It is shown that the use of a partial differential equation for the incremental pressure eliminates such a mechanism and yields a velocity-pressure coupled solution. Example flows considered are a three dimensional lid-driven cavity flow and a laminar flow through a 90 deg bend square duct. Numerical results obtained using the collocated mesh are in good agreement with the measured data and other numerical results.

Flow through three-dimensional arrangements of cylinders with alternating streamwise planar tilt

NASA Astrophysics Data System (ADS)

Sahraoui, M.; Marshall, H.; Kaviany, M.

1993-09-01

In this report, fluid flow through a three-dimensional model for the fibrous filters is examined. In this model, the three-dimensional Stokes equation with the appropriate periodic boundary conditions is solved using the finite volume method. In addition to the numerical solution, we attempt to model this flow analytically by using the two-dimensional extended analytic solution in each of the unit cells of the three-dimensional structure. Particle trajectories computed using the superimposed analytic solution of the flow field are closed to those computed using the numerical solution of the flow field. The numerical results show that the pressure drop is not affected significantly by the relative angle of rotation of the cylinders for the high porosity used in this study (epsilon = 0.8 and epsilon = 0.95). The numerical solution and the superimposed analytic solution are also compared in terms of the particle capture efficiency. The results show that the efficiency predictions using the two methods are within 10% for St = 0.01 and 5% for St = 100. As the the porosity decreases, the three-dimensional effect becomes more significant and a difference of 35% is obtained for epsilon = 0.8.

An interative solution of an integral equation for radiative transfer by using variational technique

NASA Technical Reports Server (NTRS)

Yoshikawa, K. K.

1973-01-01

An effective iterative technique is introduced to solve a nonlinear integral equation frequently associated with radiative transfer problems. The problem is formulated in such a way that each step of an iterative sequence requires the solution of a linear integral equation. The advantage of a previously introduced variational technique which utilizes a stepwise constant trial function is exploited to cope with the nonlinear problem. The method is simple and straightforward. Rapid convergence is obtained by employing a linear interpolation of the iterative solutions. Using absorption coefficients of the Milne-Eddington type, which are applicable to some planetary atmospheric radiation problems. Solutions are found in terms of temperature and radiative flux. These solutions are presented numerically and show excellent agreement with other numerical solutions.

System Simulation by Recursive Feedback: Coupling a Set of Stand-Alone Subsystem Simulations

NASA Technical Reports Server (NTRS)

Nixon, D. D.

2001-01-01

Conventional construction of digital dynamic system simulations often involves collecting differential equations that model each subsystem, arran g them to a standard form, and obtaining their numerical gin solution as a single coupled, total-system simultaneous set. Simulation by numerical coupling of independent stand-alone subsimulations is a fundamentally different approach that is attractive because, among other things, the architecture naturally facilitates high fidelity, broad scope, and discipline independence. Recursive feedback is defined and discussed as a candidate approach to multidiscipline dynamic system simulation by numerical coupling of self-contained, single-discipline subsystem simulations. A satellite motion example containing three subsystems (orbit dynamics, attitude dynamics, and aerodynamics) has been defined and constructed using this approach. Conventional solution methods are used in the subsystem simulations. Distributed and centralized implementations of coupling have been considered. Numerical results are evaluated by direct comparison with a standard total-system, simultaneous-solution approach.

Nonideal isentropic gas flow through converging-diverging nozzles

NASA Technical Reports Server (NTRS)

Bober, W.; Chow, W. L.

1990-01-01

A method for treating nonideal gas flows through converging-diverging nozzles is described. The method incorporates the Redlich-Kwong equation of state. The Runge-Kutta method is used to obtain a solution. Numerical results were obtained for methane gas. Typical plots of pressure, temperature, and area ratios as functions of Mach number are given. From the plots, it can be seen that there exists a range of reservoir conditions that require the gas to be treated as nonideal if an accurate solution is to be obtained.

GHM method for obtaining rationalsolutions of nonlinear differential equations.

PubMed

Vazquez-Leal, Hector; Sarmiento-Reyes, Arturo

2015-01-01

In this paper, we propose the application of the general homotopy method (GHM) to obtain rational solutions of nonlinear differential equations. It delivers a high precision representation of the nonlinear differential equation using a few linear algebraic terms. In order to assess the benefits of this proposal, three nonlinear problems are solved and compared against other semi-analytic methods or numerical methods. The obtained results show that GHM is a powerful tool, capable to generate highly accurate rational solutions. AMS subject classification 34L30.

The exact solution of the monoenergetic transport equation for critical cylinders

NASA Technical Reports Server (NTRS)

Westfall, R. M.; Metcalf, D. R.

1972-01-01

An analytic solution for the critical, monoenergetic, bare, infinite cylinder is presented. The solution is obtained by modifying a previous development based on a neutron density transform and Case's singular eigenfunction method. Numerical results for critical radii and the neutron density as a function of position are included and compared with the results of other methods.

Numerical Solution of the Problem of the Expansion of the Universe in the Schwarzschild Metric

NASA Astrophysics Data System (ADS)

Vasenin, I. M.; Goiko, V. L.

2018-02-01

The statement and solution of the problem of the expansion of the Universe in nonstationary spherically-symmetric coordinates in the Schwarzschild metric are considered without pressure taken into account. The observational data of astronomers investigating the rates of recession of distant stars are explained on the basis of the obtained solutions.

On the Possibilities of Predicting Geomagnetic Secular Variation with Geodynamo Modeling

NASA Technical Reports Server (NTRS)

Kuang, Wei-Jia; Tangborn, Andrew; Sabaka, Terrance

2004-01-01

We use our MoSST core dynamics model and geomagnetic field at the core-mantle boundary (CMB) continued downward from surface observations to investigate possibilities of geomagnetic data assimilation, so that model results and current geomagnetic observations can be used to predict geomagnetic secular variation in future. As the first attempt, we apply data insertion technique to examine evolution of the model solution that is modified by geomagnetic input. Our study demonstrate that, with a single data insertion, large-scale poloidal magnetic field obtained from subsequent numerical simulation evolves similarly to the observed geomagnetic variation, regardless of the initial choice of the model solution (so long it is a well developed numerical solution). The model solution diverges on the time scales on the order of 60 years, similar to the time scales of the torsional oscillations in the Earth's core. Our numerical test shows that geomagnetic data assimilation is promising with our MoSST model.

Exact solutions of the population balance equation including particle transport, using group analysis

NASA Astrophysics Data System (ADS)

Lin, Fubiao; Meleshko, Sergey V.; Flood, Adrian E.

2018-06-01

The population balance equation (PBE) has received an unprecedented amount of attention in recent years from both academics and industrial practitioners because of its long history, widespread use in engineering, and applicability to a wide variety of particulate and discrete-phase processes. However it is typically impossible to obtain analytical solutions, although in almost every case a numerical solution of the PBEs can be obtained. In this article, the symmetries of PBEs with homogeneous coagulation kernels involving aggregation, breakage and growth processes and particle transport in one dimension are found by direct solving the determining equations. Using the optimal system of one and two-dimensional subalgebras, all invariant solutions and reduced equations are obtained. In particular, an explicit analytical physical solution is also presented.

Nonlinear Earthquake Analysis of Reinforced Concrete Frames with Fiber and Bernoulli-Euler Beam-Column Element

PubMed Central

Karaton, Muhammet

2014-01-01

A beam-column element based on the Euler-Bernoulli beam theory is researched for nonlinear dynamic analysis of reinforced concrete (RC) structural element. Stiffness matrix of this element is obtained by using rigidity method. A solution technique that included nonlinear dynamic substructure procedure is developed for dynamic analyses of RC frames. A predicted-corrected form of the Bossak-α method is applied for dynamic integration scheme. A comparison of experimental data of a RC column element with numerical results, obtained from proposed solution technique, is studied for verification the numerical solutions. Furthermore, nonlinear cyclic analysis results of a portal reinforced concrete frame are achieved for comparing the proposed solution technique with Fibre element, based on flexibility method. However, seismic damage analyses of an 8-story RC frame structure with soft-story are investigated for cases of lumped/distributed mass and load. Damage region, propagation, and intensities according to both approaches are researched. PMID:24578667

Approximate Solution Methods for Spectral Radiative Transfer in High Refractive Index Layers

NASA Technical Reports Server (NTRS)

Siegel, R.; Spuckler, C. M.

1994-01-01

Some ceramic materials for high temperature applications are partially transparent for radiative transfer. The refractive indices of these materials can be substantially greater than one which influences internal radiative emission and reflections. Heat transfer behavior of single and laminated layers has been obtained in the literature by numerical solutions of the radiative transfer equations coupled with heat conduction and heating at the boundaries by convection and radiation. Two-flux and diffusion methods are investigated here to obtain approximate solutions using a simpler formulation than required for exact numerical solutions. Isotropic scattering is included. The two-flux method for a single layer yields excellent results for gray and two band spectral calculations. The diffusion method yields a good approximation for spectral behavior in laminated multiple layers if the overall optical thickness is larger than about ten. A hybrid spectral model is developed using the two-flux method in the optically thin bands, and radiative diffusion in bands that are optically thick.

Mechanics of additively manufactured porous biomaterials based on the rhombicuboctahedron unit cell.

PubMed

Hedayati, R; Sadighi, M; Mohammadi-Aghdam, M; Zadpoor, A A

2016-01-01

Thanks to recent developments in additive manufacturing techniques, it is now possible to fabricate porous biomaterials with arbitrarily complex micro-architectures. Micro-architectures of such biomaterials determine their physical and biological properties, meaning that one could potentially improve the performance of such biomaterials through rational design of micro-architecture. The relationship between the micro-architecture of porous biomaterials and their physical and biological properties has therefore received increasing attention recently. In this paper, we studied the mechanical properties of porous biomaterials made from a relatively unexplored unit cell, namely rhombicuboctahedron. We derived analytical relationships that relate the micro-architecture of such porous biomaterials, i.e. the dimensions of the rhombicuboctahedron unit cell, to their elastic modulus, Poisson's ratio, and yield stress. Finite element models were also developed to validate the analytical solutions. Analytical and numerical results were compared with experimental data from one of our recent studies. It was found that analytical solutions and numerical results show a very good agreement particularly for smaller values of apparent density. The elastic moduli predicted by analytical and numerical models were in very good agreement with experimental observations too. While in excellent agreement with each other, analytical and numerical models somewhat over-predicted the yield stress of the porous structures as compared to experimental data. As the ratio of the vertical struts to the inclined struts, α, approaches zero and infinity, the rhombicuboctahedron unit cell respectively approaches the octahedron (or truncated cube) and cube unit cells. For those limits, the analytical solutions presented here were found to approach the analytic solutions obtained for the octahedron, truncated cube, and cube unit cells, meaning that the presented solutions are generalizations of the analytical solutions obtained for several other types of porous biomaterials. Copyright © 2015 Elsevier Ltd. All rights reserved.

Second order finite-difference ghost-point multigrid methods for elliptic problems with discontinuous coefficients on an arbitrary interface

NASA Astrophysics Data System (ADS)

Coco, Armando; Russo, Giovanni

2018-05-01

In this paper we propose a second-order accurate numerical method to solve elliptic problems with discontinuous coefficients (with general non-homogeneous jumps in the solution and its gradient) in 2D and 3D. The method consists of a finite-difference method on a Cartesian grid in which complex geometries (boundaries and interfaces) are embedded, and is second order accurate in the solution and the gradient itself. In order to avoid the drop in accuracy caused by the discontinuity of the coefficients across the interface, two numerical values are assigned on grid points that are close to the interface: a real value, that represents the numerical solution on that grid point, and a ghost value, that represents the numerical solution extrapolated from the other side of the interface, obtained by enforcing the assigned non-homogeneous jump conditions on the solution and its flux. The method is also extended to the case of matrix coefficient. The linear system arising from the discretization is solved by an efficient multigrid approach. Unlike the 1D case, grid points are not necessarily aligned with the normal derivative and therefore suitable stencils must be chosen to discretize interface conditions in order to achieve second order accuracy in the solution and its gradient. A proper treatment of the interface conditions will allow the multigrid to attain the optimal convergence factor, comparable with the one obtained by Local Fourier Analysis for rectangular domains. The method is robust enough to handle large jump in the coefficients: order of accuracy, monotonicity of the errors and good convergence factor are maintained by the scheme.

Steady states and outbreaks of two-phase nonlinear age-structured model of population dynamics with discrete time delay.

PubMed

Akimenko, Vitalii; Anguelov, Roumen

2017-12-01

In this paper we study the nonlinear age-structured model of a polycyclic two-phase population dynamics including delayed effect of population density growth on the mortality. Both phases are modelled as a system of initial boundary values problem for semi-linear transport equation with delay and initial problem for nonlinear delay ODE. The obtained system is studied both theoretically and numerically. Three different regimes of population dynamics for asymptotically stable states of autonomous systems are obtained in numerical experiments for the different initial values of population density. The quasi-periodical travelling wave solutions are studied numerically for the autonomous system with the different values of time delays and for the system with oscillating death rate and birth modulus. In both cases it is observed three types of travelling wave solutions: harmonic oscillations, pulse sequence and single pulse.

Large eddy simulation of incompressible turbulent channel flow

NASA Technical Reports Server (NTRS)

Moin, P.; Reynolds, W. C.; Ferziger, J. H.

1978-01-01

The three-dimensional, time-dependent primitive equations of motion were numerically integrated for the case of turbulent channel flow. A partially implicit numerical method was developed. An important feature of this scheme is that the equation of continuity is solved directly. The residual field motions were simulated through an eddy viscosity model, while the large-scale field was obtained directly from the solution of the governing equations. An important portion of the initial velocity field was obtained from the solution of the linearized Navier-Stokes equations. The pseudospectral method was used for numerical differentiation in the horizontal directions, and second-order finite-difference schemes were used in the direction normal to the walls. The large eddy simulation technique is capable of reproducing some of the important features of wall-bounded turbulent flows. The resolvable portions of the root-mean square wall pressure fluctuations, pressure velocity-gradient correlations, and velocity pressure-gradient correlations are documented.

Simplified method for numerical modeling of fiber lasers.

PubMed

Shtyrina, O V; Yarutkina, I A; Fedoruk, M P

2014-12-29

A simplified numerical approach to modeling of dissipative dispersion-managed fiber lasers is examined. We present a new numerical iteration algorithm for finding the periodic solutions of the system of nonlinear ordinary differential equations describing the intra-cavity dynamics of the dissipative soliton characteristics in dispersion-managed fiber lasers. We demonstrate that results obtained using simplified model are in good agreement with full numerical modeling based on the corresponding partial differential equations.

Numerical Model of Multiple Scattering and Emission from Layering Snowpack for Microwave Remote Sensing

NASA Astrophysics Data System (ADS)

Jin, Y.; Liang, Z.

2002-12-01

The vector radiative transfer (VRT) equation is an integral-deferential equation to describe multiple scattering, absorption and transmission of four Stokes parameters in random scatter media. From the integral formal solution of VRT equation, the lower order solutions, such as the first-order scattering for a layer medium or the second order scattering for a half space, can be obtained. The lower order solutions are usually good at low frequency when high-order scattering is negligible. It won't be feasible to continue iteration for obtaining high order scattering solution because too many folds integration would be involved. In the space-borne microwave remote sensing, for example, the DMSP (Defense Meterological Satellite Program) SSM/I (Special Sensor Microwave/Imager) employed seven channels of 19, 22, 37 and 85GHz. Multiple scattering from the terrain surfaces such as snowpack cannot be neglected at these channels. The discrete ordinate and eigen-analysis method has been studied to take into account for multiple scattering and applied to remote sensing of atmospheric precipitation, snowpack etc. Snowpack was modeled as a layer of dense spherical particles, and the VRT for a layer of uniformly dense spherical particles has been numerically studied by the discrete ordinate method. However, due to surface melting and refrozen crusts, the snowpack undergoes stratifying to form inhomegeneous profiles of the ice grain size, fractional volume and physical temperature etc. It becomes necessary to study multiple scattering and emission from stratified snowpack of dense ice grains. But, the discrete ordinate and eigen-analysis method cannot be simply applied to multi-layers model, because numerically solving a set of multi-equations of VRT is difficult. Stratifying the inhomogeneous media into multi-slabs and employing the first order Mueller matrix of each thin slab, this paper developed an iterative method to derive high orders scattering solutions of whole scatter media. High order scattering and emission from inhomogeneous stratifying media of dense spherical particles are numerically obtained. The brightness temperature at low frequency such as 5.3 GHz without high order scattering and at SSM/I channels with high order scattering are obtained. This approach is also compared with the conventional discrete ordinate method for an uniform layer model. Numerical simulation for inhomogeneous snowpack is also compared with the measurements of microwave remote sensing.

Numerical Issues for Circulation Control Calculations

NASA Technical Reports Server (NTRS)

Swanson, Roy C., Jr.; Rumsey, Christopher L.

2006-01-01

Steady-state and time-accurate two-dimensional solutions of the compressible Reynolds-averaged Navier- Stokes equations are obtained for flow over the Lockheed circulation control (CC) airfoil and the General Aviation CC (GACC) airfoil. Numerical issues in computing circulation control flows such as the effects of grid resolution, boundary and initial conditions, and unsteadiness are addressed. For the Lockheed CC airfoil computed solutions are compared with detailed experimental data, which include velocity and Reynolds stress profiles. Three turbulence models, having either one or two transport equations, are considered. Solutions are obtained on a sequence of meshes, with mesh refinement primarily concentrated on the airfoil circular trailing edge. Several effects related to mesh refinement are identified. For example, sometimes sufficient mesh resolution can exclude nonphysical solutions, which can occur in CC airfoil calculations. Also, sensitivities of the turbulence models with mesh refinement are discussed. In the case of the GACC airfoil the focus is on the difference between steady-state and time-accurate solutions. A specific objective is to determine if there is self-excited vortex shedding from the jet slot lip.

Numerical Solution of the Flow of a Perfect Gas Over A Circular Cylinder at Infinite Mach Number

NASA Technical Reports Server (NTRS)

Hamaker, Frank M.

1959-01-01

A solution for the two-dimensional flow of an inviscid perfect gas over a circular cylinder at infinite Mach number is obtained by numerical methods of analysis. Nonisentropic conditions of curved shock waves and vorticity are included in the solution. The analysis is divided into two distinct regions, the subsonic region which is analyzed by the relaxation method of Southwell and the supersonic region which was treated by the method of characteristics. Both these methods of analysis are inapplicable on the sonic line which is therefore considered separately. The shapes of the sonic line and the shock wave are obtained by iteration techniques. The striking result of the solution is the strong curvature of the sonic line and of the other lines of constant Mach number. Because of this the influence of the supersonic flow on the sonic line is negligible. On comparison with Newtonian flow methods, it is found that the approximate methods show a larger variation of surface pressure than is given by the present solution.

Nonlinear dynamics and numerical uncertainties in CFD

NASA Technical Reports Server (NTRS)

Yee, H. C.; Sweby, P. K.

1996-01-01

The application of nonlinear dynamics to improve the understanding of numerical uncertainties in computational fluid dynamics (CFD) is reviewed. Elementary examples in the use of dynamics to explain the nonlinear phenomena and spurious behavior that occur in numerics are given. The role of dynamics in the understanding of long time behavior of numerical integrations and the nonlinear stability, convergence, and reliability of using time-marching, approaches for obtaining steady-state numerical solutions in CFD is explained. The study is complemented with spurious behavior observed in CFD computations.

Numerical solution of the unsteady diffusion-convection-reaction equation based on improved spectral Galerkin method

NASA Astrophysics Data System (ADS)

Zhong, Jiaqi; Zeng, Cheng; Yuan, Yupeng; Zhang, Yuzhe; Zhang, Ye

2018-04-01

The aim of this paper is to present an explicit numerical algorithm based on improved spectral Galerkin method for solving the unsteady diffusion-convection-reaction equation. The principal characteristics of this approach give the explicit eigenvalues and eigenvectors based on the time-space separation method and boundary condition analysis. With the help of Fourier series and Galerkin truncation, we can obtain the finite-dimensional ordinary differential equations which facilitate the system analysis and controller design. By comparing with the finite element method, the numerical solutions are demonstrated via two examples. It is shown that the proposed method is effective.

Numerical Simulation of the Flow over a Segment-Conical Body on the Basis of Reynolds Equations

NASA Astrophysics Data System (ADS)

Egorov, I. V.; Novikov, A. V.; Palchekovskaya, N. V.

2018-01-01

Numerical simulation was used to study the 3D supersonic flow over a segment-conical body similar in shape to the ExoMars space vehicle. The nonmonotone behavior of the normal force acting on the body placed in a supersonic gas flow was analyzed depending on the angle of attack. The simulation was based on the numerical solution of the unsteady Reynolds-averaged Navier-Stokes equations with a two-parameter differential turbulence model. The solution of the problem was obtained using the in-house solver HSFlow with an efficient parallel algorithm intended for multiprocessor super computers.

An unstaggered central scheme on nonuniform grids for the simulation of a compressible two-phase flow model

DOE Office of Scientific and Technical Information (OSTI.GOV)

Touma, Rony; Zeidan, Dia

In this paper we extend a central finite volume method on nonuniform grids to the case of drift-flux two-phase flow problems. The numerical base scheme is an unstaggered, non oscillatory, second-order accurate finite volume scheme that evolves a piecewise linear numerical solution on a single grid and uses dual cells intermediately while updating the numerical solution to avoid the resolution of the Riemann problems arising at the cell interfaces. We then apply the numerical scheme and solve a classical drift-flux problem. The obtained results are in good agreement with corresponding ones appearing in the recent literature, thus confirming the potentialmore » of the proposed scheme.« less

Comments on real tachyon vacuum solution without square roots

NASA Astrophysics Data System (ADS)

Arroyo, E. Aldo

2018-01-01

We analyze the consistency of a recently proposed real tachyon vacuum solution without square roots in open bosonic string field theory. We show that the equation of motion contracted with the solution itself is satisfied. Additionally, by expanding the solution in the basis of the curly ℒ0 and the traditional L 0 eigenstates, we evaluate numerically the vacuum energy and obtain a result in agreement with Sen's conjecture.

Approximated analytical solution to an Ebola optimal control problem

NASA Astrophysics Data System (ADS)

Hincapié-Palacio, Doracelly; Ospina, Juan; Torres, Delfim F. M.

2016-11-01

An analytical expression for the optimal control of an Ebola problem is obtained. The analytical solution is found as a first-order approximation to the Pontryagin Maximum Principle via the Euler-Lagrange equation. An implementation of the method is given using the computer algebra system Maple. Our analytical solutions confirm the results recently reported in the literature using numerical methods.

A preconditioned formulation of the Cauchy-Riemann equations

NASA Technical Reports Server (NTRS)

Phillips, T. N.

1983-01-01

A preconditioning of the Cauchy-Riemann equations which results in a second-order system is described. This system is shown to have a unique solution if the boundary conditions are chosen carefully. This choice of boundary condition enables the solution of the first-order system to be retrieved. A numerical solution of the preconditioned equations is obtained by the multigrid method.

Essentially nonoscillatory postprocessing filtering methods

NASA Technical Reports Server (NTRS)

Lafon, F.; Osher, S.

1992-01-01

High order accurate centered flux approximations used in the computation of numerical solutions to nonlinear partial differential equations produce large oscillations in regions of sharp transitions. Here, we present a new class of filtering methods denoted by Essentially Nonoscillatory Least Squares (ENOLS), which constructs an upgraded filtered solution that is close to the physically correct weak solution of the original evolution equation. Our method relies on the evaluation of a least squares polynomial approximation to oscillatory data using a set of points which is determined via the ENO network. Numerical results are given in one and two space dimensions for both scalar and systems of hyperbolic conservation laws. Computational running time, efficiency, and robustness of method are illustrated in various examples such as Riemann initial data for both Burgers' and Euler's equations of gas dynamics. In all standard cases, the filtered solution appears to converge numerically to the correct solution of the original problem. Some interesting results based on nonstandard central difference schemes, which exactly preserve entropy, and have been recently shown generally not to be weakly convergent to a solution of the conservation law, are also obtained using our filters.

ADM For Solving Linear Second-Order Fredholm Integro-Differential Equations

NASA Astrophysics Data System (ADS)

Karim, Mohd F.; Mohamad, Mahathir; Saifullah Rusiman, Mohd; Che-Him, Norziha; Roslan, Rozaini; Khalid, Kamil

2018-04-01

In this paper, we apply Adomian Decomposition Method (ADM) as numerically analyse linear second-order Fredholm Integro-differential Equations. The approximate solutions of the problems are calculated by Maple package. Some numerical examples have been considered to illustrate the ADM for solving this equation. The results are compared with the existing exact solution. Thus, the Adomian decomposition method can be the best alternative method for solving linear second-order Fredholm Integro-Differential equation. It converges to the exact solution quickly and in the same time reduces computational work for solving the equation. The result obtained by ADM shows the ability and efficiency for solving these equations.

Numerical calculation of flow fields about rectangular wings of finite thickness in supersonic flow. Ph.D. Thesis

NASA Technical Reports Server (NTRS)

Vogel, J. M.

1973-01-01

The calculation of the outer inviscid flow about a rectangular wing moving at supersonic speeds is reported. The inviscid equations of motion governing the flow generated by the wing form a set of hyperbolic differential equations. The flow field about the rectangular wing is separated into three regions consisting of the forebody, the afterbody, and the wing wake. Solutions for the forebody are obtained using conical flow techniques while the afterbody and the wing wake regions are treated as initial value problems. The numerical solutions are compared in the two dimensional regions with known exact solutions.

Singular perturbation solutions of steady-state Poisson-Nernst-Planck systems.

PubMed

Wang, Xiang-Sheng; He, Dongdong; Wylie, Jonathan J; Huang, Huaxiong

2014-02-01

We study the Poisson-Nernst-Planck (PNP) system with an arbitrary number of ion species with arbitrary valences in the absence of fixed charges. Assuming point charges and that the Debye length is small relative to the domain size, we derive an asymptotic formula for the steady-state solution by matching outer and boundary layer solutions. The case of two ionic species has been extensively studied, the uniqueness of the solution has been proved, and an explicit expression for the solution has been obtained. However, the case of three or more ions has received significantly less attention. Previous work has indicated that the solution may be nonunique and that even obtaining numerical solutions is a difficult task since one must solve complicated systems of nonlinear equations. By adopting a methodology that preserves the symmetries of the PNP system, we show that determining the outer solution effectively reduces to solving a single scalar transcendental equation. Due to the simple form of the transcendental equation, it can be solved numerically in a straightforward manner. Our methodology thus provides a standard procedure for solving the PNP system and we illustrate this by solving some practical examples. Despite the fact that for three ions, previous studies have indicated that multiple solutions may exist, we show that all except for one of these solutions are unphysical and thereby prove the existence and uniqueness for the three-ion case.

Transitions of stationary to pulsating solutions in the complex cubic-quintic Ginzburg-Landau equation under the influence of nonlinear gain and higher-order effects

NASA Astrophysics Data System (ADS)

Uzunov, Ivan M.; Georgiev, Zhivko D.; Arabadzhiev, Todor N.

2018-05-01

In this paper we study the transitions of stationary to pulsating solutions in the complex cubic-quintic Ginzburg-Landau equation (CCQGLE) under the influence of nonlinear gain, its saturation, and higher-order effects: self-steepening, third-order of dispersion, and intrapulse Raman scattering in the anomalous dispersion region. The variation method and the method of moments are applied in order to obtain the dynamic models with finite degrees of freedom for the description of stationary and pulsating solutions. Having applied the first model and its bifurcation analysis we have discovered the existence of families of subcritical Poincaré-Andronov-Hopf bifurcations due to the intrapulse Raman scattering, as well as some small nonlinear gain and the saturation of the nonlinear gain. A phenomenon of nonlinear stability has been studied and it has been shown that long living pulsating solutions with relatively small fluctuations of amplitude and frequencies exist at the bifurcation point. The numerical analysis of the second model has revealed the existence of Poincaré-Andronov-Hopf bifurcations of Raman dissipative soliton under the influence of the self-steepening effect and large nonlinear gain. All our theoretical predictions have been confirmed by the direct numerical solution of the full perturbed CCQGLE. The detailed comparison between the results obtained by both dynamic models and the direct numerical solution of the perturbed CCQGLE has proved the applicability of the proposed models in the investigation of the solutions of the perturbed CCQGLE.

Drifting solutions with elliptic symmetry for the compressible Navier-Stokes equations with density-dependent viscosity

DOE Office of Scientific and Technical Information (OSTI.GOV)

An, Hongli, E-mail: kaixinguoan@163.com; Yuen, Manwai, E-mail: nevetsyuen@hotmail.com

2014-05-15

In this paper, we investigate the analytical solutions of the compressible Navier-Stokes equations with dependent-density viscosity. By using the characteristic method, we successfully obtain a class of drifting solutions with elliptic symmetry for the Navier-Stokes model wherein the velocity components are governed by a generalized Emden dynamical system. In particular, when the viscosity variables are taken the same as Yuen [M. W. Yuen, “Analytical solutions to the Navier-Stokes equations,” J. Math. Phys. 49, 113102 (2008)], our solutions constitute a generalization of that obtained by Yuen. Interestingly, numerical simulations show that the analytical solutions can be used to explain the driftingmore » phenomena of the propagation wave like Tsunamis in oceans.« less

A Generalization of the Karush-Kuhn-Tucker Theorem for Approximate Solutions of Mathematical Programming Problems Based on Quadratic Approximation

NASA Astrophysics Data System (ADS)

Voloshinov, V. V.

2018-03-01

In computations related to mathematical programming problems, one often has to consider approximate, rather than exact, solutions satisfying the constraints of the problem and the optimality criterion with a certain error. For determining stopping rules for iterative procedures, in the stability analysis of solutions with respect to errors in the initial data, etc., a justified characteristic of such solutions that is independent of the numerical method used to obtain them is needed. A necessary δ-optimality condition in the smooth mathematical programming problem that generalizes the Karush-Kuhn-Tucker theorem for the case of approximate solutions is obtained. The Lagrange multipliers corresponding to the approximate solution are determined by solving an approximating quadratic programming problem.

Complexiton and solitary wave solutions of the coupled nonlinear Maccari’s system using two integration schemes

NASA Astrophysics Data System (ADS)

Inc, Mustafa; Aliyu, Aliyu Isa; Yusuf, Abdullahi; Baleanu, Dumitru; Nuray, Elif

2018-01-01

In this paper, we consider a coupled nonlinear Maccari’s system (CNMS) which describes the motion of isolated waves localized in a small part of space. There are some integration tools that are adopted to retrieve the solitary wave solutions. They are the modified F-Expansion and the generalized projective Riccati equation methods. Topological, non-topological, complexiton, singular and trigonometric function solutions are derived. A comparison between the results in this paper and the well-known results in the literature is also given. The derived structures of the obtained solutions offer a rich platform to study the nonlinear CNMS. Numerical simulation of the obtained solutions are presented with interesting figures showing the physical meaning of the solutions.

Exact analytical solutions of continuity equation for electron beams precipitating in Coulomb collisions

DOE Office of Scientific and Technical Information (OSTI.GOV)

Dobranskis, R. R.; Zharkova, V. V., E-mail: valentina.zharkova@northumbria.ac.uk

2014-06-10

The original continuity equation (CE) used for the interpretation of the power law energy spectra of beam electrons in flares was written and solved for an electron beam flux while ignoring an additional free term with an electron density. In order to remedy this omission, the original CE for electron flux, considering beam's energy losses in Coulomb collisions, was first differentiated by the two independent variables: depth and energy leading to partial differential equation for an electron beam density instead of flux with the additional free term. The analytical solution of this partial differential continuity equation (PDCE) is obtained bymore » using the method of characteristics. This solution is further used to derive analytical expressions for mean electron spectra for Coulomb collisions and to carry out numeric calculations of hard X-ray (HXR) photon spectra for beams with different parameters. The solutions revealed a significant departure of electron densities at lower energies from the original results derived from the CE for the flux obtained for Coulomb collisions. This departure is caused by the additional exponential term that appeared in the updated solutions for electron differential density leading to its faster decrease at lower energies (below 100 keV) with every precipitation depth similar to the results obtained with numerical Fokker-Planck solutions. The effects of these updated solutions for electron densities on mean electron spectra and HXR photon spectra are also discussed.« less

On the Unreasonable Effectiveness of post-Newtonian Theory in Gravitational-Wave Physics

ScienceCinema

Will, Clifford M.

2017-12-22

The first indirect detection of gravitational waves involved a binary system of neutron stars.  In the future, the first direct detection may also involve binary systems -- inspiralling and merging binary neutron stars or black holes. This means that it is essential to understand in full detail the two-body system in general relativity, a notoriously difficult problem with a long history. Post-Newtonian approximation methods are thought to work only under slow motion and weak field conditions, while numerical solutions of Einstein's equations are thought to be limited to the final merger phase.  Recent results have shown that post-Newtonian approximations seem to remain unreasonably valid well into the relativistic regime, while advances in numerical relativity now permit solutions for numerous orbits before merger.  It is now possible to envision linking post-Newtonian theory and numerical relativity to obtain a complete "solution" of the general relativistic two-body problem.  These solutions will play a central role in detecting and understanding gravitational wave signals received by interferometric observatories on Earth and in space.

Numerical method for the solution of large systems of differential equations of the boundary layer type

NASA Technical Reports Server (NTRS)

Green, M. J.; Nachtsheim, P. R.

1972-01-01

A numerical method for the solution of large systems of nonlinear differential equations of the boundary-layer type is described. The method is a modification of the technique for satisfying asymptotic boundary conditions. The present method employs inverse interpolation instead of the Newton method to adjust the initial conditions of the related initial-value problem. This eliminates the so-called perturbation equations. The elimination of the perturbation equations not only reduces the user's preliminary work in the application of the method, but also reduces the number of time-consuming initial-value problems to be numerically solved at each iteration. For further ease of application, the solution of the overdetermined system for the unknown initial conditions is obtained automatically by applying Golub's linear least-squares algorithm. The relative ease of application of the proposed numerical method increases directly as the order of the differential-equation system increases. Hence, the method is especially attractive for the solution of large-order systems. After the method is described, it is applied to a fifth-order problem from boundary-layer theory.

Direct numerical solution of the Ornstein-Zernike integral equation and spatial distribution of water around hydrophobic molecules

NASA Astrophysics Data System (ADS)

Ikeguchi, Mitsunori; Doi, Junta

1995-09-01

The Ornstein-Zernike integral equation (OZ equation) has been used to evaluate the distribution function of solvents around solutes, but its numerical solution is difficult for molecules with a complicated shape. This paper proposes a numerical method to directly solve the OZ equation by introducing the 3D lattice. The method employs no approximation the reference interaction site model (RISM) equation employed. The method enables one to obtain the spatial distribution of spherical solvents around solutes with an arbitrary shape. Numerical accuracy is sufficient when the grid-spacing is less than 0.5 Å for solvent water. The spatial water distribution around a propane molecule is demonstrated as an example of a nonspherical hydrophobic molecule using iso-value surfaces. The water model proposed by Pratt and Chandler is used. The distribution agrees with the molecular dynamics simulation. The distribution increases offshore molecular concavities. The spatial distribution of water around 5α-cholest-2-ene (C27H46) is visualized using computer graphics techniques and a similar trend is observed.

Application of CFD to aerothermal heating problems

NASA Technical Reports Server (NTRS)

Macaraeg, M. G.

1986-01-01

Numerical solutions of the compressible Navier-Stokes equations by an alternating direction implicit scheme, applied to two experimental investigations are presented. The first is cooling by injection of a gas jet through the nose of an ogive-cone, and the second is the aerothermal environment in the gap formed by the wing and elevon section of a test model of the space shuttle. The simulations demonstrate that accurate pressure calculations are easily obtained on a coarse grid, while convergence is obtained after the residual reduces by four orders of magnitude. Accurate heating rates, however, require a fine grid solution, with convergence requiring at least a reduction of six orders of magnitude in the residual. The effect of artificial dissipation on numerical results is also assessed.

Soliton and kink jams in traffic flow with open boundaries.

PubMed

Muramatsu, M; Nagatani, T

1999-07-01

Soliton density wave is investigated numerically and analytically in the optimal velocity model (a car-following model) of a one-dimensional traffic flow with open boundaries. Soliton density wave is distinguished from the kink density wave. It is shown that the soliton density wave appears only at the threshold of occurrence of traffic jams. The Korteweg-de Vries (KdV) equation is derived from the optimal velocity model by the use of the nonlinear analysis. It is found that the traffic soliton appears only near the neutral stability line. The soliton solution is analytically obtained from the perturbed KdV equation. It is shown that the soliton solution obtained from the nonlinear analysis is consistent with that of the numerical simulation.

Controllable parabolic-cylinder optical rogue wave.

PubMed

Zhong, Wei-Ping; Chen, Lang; Belić, Milivoj; Petrović, Nikola

2014-10-01

We demonstrate controllable parabolic-cylinder optical rogue waves in certain inhomogeneous media. An analytical rogue wave solution of the generalized nonlinear Schrödinger equation with spatially modulated coefficients and an external potential in the form of modulated quadratic potential is obtained by the similarity transformation. Numerical simulations are performed for comparison with the analytical solutions and to confirm the stability of the rogue wave solution obtained. These optical rogue waves are built by the products of parabolic-cylinder functions and the basic rogue wave solution of the standard nonlinear Schrödinger equation. Such rogue waves may appear in different forms, as the hump and paw profiles.

Hopf solitons in the Nicole model

DOE Office of Scientific and Technical Information (OSTI.GOV)

Gillard, Mike; Sutcliffe, Paul

2010-12-15

The Nicole model is a conformal field theory in a three-dimensional space. It has topological soliton solutions classified by the integer-valued Hopf charge, and all currently known solitons are axially symmetric. A volume-preserving flow is used to construct soliton solutions numerically for all Hopf charges from 1 to 8. It is found that the known axially symmetric solutions are unstable for Hopf charges greater than 2 and new lower energy solutions are obtained that include knots and links. A comparison with the Skyrme-Faddeev model suggests many universal features, though there are some differences in the link types obtained in themore » two theories.« less

On the solution of the complex eikonal equation in acoustic VTI media: A perturbation plus optimization scheme

NASA Astrophysics Data System (ADS)

Huang, Xingguo; Sun, Jianguo; Greenhalgh, Stewart

2018-04-01

We present methods for obtaining numerical and analytic solutions of the complex eikonal equation in inhomogeneous acoustic VTI media (transversely isotropic media with a vertical symmetry axis). The key and novel point of the method for obtaining numerical solutions is to transform the problem of solving the highly nonlinear acoustic VTI eikonal equation into one of solving the relatively simple eikonal equation for the background (isotropic) medium and a system of linear partial differential equations. Specifically, to obtain the real and imaginary parts of the complex traveltime in inhomogeneous acoustic VTI media, we generalize a perturbation theory, which was developed earlier for solving the conventional real eikonal equation in inhomogeneous anisotropic media, to the complex eikonal equation in such media. After the perturbation analysis, we obtain two types of equations. One is the complex eikonal equation for the background medium and the other is a system of linearized partial differential equations for the coefficients of the corresponding complex traveltime formulas. To solve the complex eikonal equation for the background medium, we employ an optimization scheme that we developed for solving the complex eikonal equation in isotropic media. Then, to solve the system of linearized partial differential equations for the coefficients of the complex traveltime formulas, we use the finite difference method based on the fast marching strategy. Furthermore, by applying the complex source point method and the paraxial approximation, we develop the analytic solutions of the complex eikonal equation in acoustic VTI media, both for the isotropic and elliptical anisotropic background medium. Our numerical results demonstrate the effectiveness of our derivations and illustrate the influence of the beam widths and the anisotropic parameters on the complex traveltimes.

Fully pseudospectral solution of the conformally invariant wave equation near the cylinder at spacelike infinity. III: nonspherical Schwarzschild waves and singularities at null infinity

NASA Astrophysics Data System (ADS)

Frauendiener, Jörg; Hennig, Jörg

2018-03-01

We extend earlier numerical and analytical considerations of the conformally invariant wave equation on a Schwarzschild background from the case of spherically symmetric solutions, discussed in Frauendiener and Hennig (2017 Class. Quantum Grav. 34 045005), to the case of general, nonsymmetric solutions. A key element of our approach is the modern standard representation of spacelike infinity as a cylinder. With a decomposition into spherical harmonics, we reduce the four-dimensional wave equation to a family of two-dimensional equations. These equations can be used to study the behaviour at the cylinder, where the solutions turn out to have, in general, logarithmic singularities at infinitely many orders. We derive regularity conditions that may be imposed on the initial data, in order to avoid the first singular terms. We then demonstrate that the fully pseudospectral time evolution scheme can be applied to this problem leading to a highly accurate numerical reconstruction of the nonsymmetric solutions. We are particularly interested in the behaviour of the solutions at future null infinity, and we numerically show that the singularities spread to null infinity from the critical set, where the cylinder approaches null infinity. The observed numerical behaviour is consistent with similar logarithmic singularities found analytically on the critical set. Finally, we demonstrate that even solutions with singularities at low orders can be obtained with high accuracy by virtue of a coordinate transformation that converts solutions with logarithmic singularities into smooth solutions.

Eikonal solutions to optical model coupled-channel equations

NASA Technical Reports Server (NTRS)

Cucinotta, Francis A.; Khandelwal, Govind S.; Maung, Khin M.; Townsend, Lawrence W.; Wilson, John W.

1988-01-01

Methods of solution are presented for the Eikonal form of the nucleus-nucleus coupled-channel scattering amplitudes. Analytic solutions are obtained for the second-order optical potential for elastic scattering. A numerical comparison is made between the first and second order optical model solutions for elastic and inelastic scattering of H-1 and He-4 on C-12. The effects of bound-state excitations on total and reaction cross sections are also estimated.

Study of stability of the difference scheme for the model problem of the gaslift process

NASA Astrophysics Data System (ADS)

Temirbekov, Nurlan; Turarov, Amankeldy

2017-09-01

The paper studies a model of the gaslift process where the motion in a gas-lift well is described by partial differential equations. The system describing the studied process consists of equations of motion, continuity, equations of thermodynamic state, and hydraulic resistance. A two-layer finite-difference Lax-Vendroff scheme is constructed for the numerical solution of the problem. The stability of the difference scheme for the model problem is investigated using the method of a priori estimates, the order of approximation is investigated, the algorithm for numerical implementation of the gaslift process model is given, and the graphs are presented. The development and investigation of difference schemes for the numerical solution of systems of equations of gas dynamics makes it possible to obtain simultaneously exact and monotonic solutions.

Accurate spectral solutions for the parabolic and elliptic partial differential equations by the ultraspherical tau method

NASA Astrophysics Data System (ADS)

Doha, E. H.; Abd-Elhameed, W. M.

2005-09-01

We present a double ultraspherical spectral methods that allow the efficient approximate solution for the parabolic partial differential equations in a square subject to the most general inhomogeneous mixed boundary conditions. The differential equations with their boundary and initial conditions are reduced to systems of ordinary differential equations for the time-dependent expansion coefficients. These systems are greatly simplified by using tensor matrix algebra, and are solved by using the step-by-step method. Numerical applications of how to use these methods are described. Numerical results obtained compare favorably with those of the analytical solutions. Accurate double ultraspherical spectral approximations for Poisson's and Helmholtz's equations are also noted. Numerical experiments show that spectral approximation based on Chebyshev polynomials of the first kind is not always better than others based on ultraspherical polynomials.

Numerical and Experimental Studies of the Natural Convection Flow Within a Horizontal Cylinder Subjected to a Uniformly Cold Wall Boundary Condition. Ph.D. Thesis - Va. Poly. Inst. and State Univ.

NASA Technical Reports Server (NTRS)

Stewart, R. B.

1972-01-01

Numberical solutions are obtained for the quasi-compressible Navier-Stokes equations governing the time dependent natural convection flow within a horizontal cylinder. The early time flow development and wall heat transfer is obtained after imposing a uniformly cold wall boundary condition on the cylinder. Solutions are also obtained for the case of a time varying cold wall boundary condition. Windware explicit differ-encing is used for the numerical solutions. The viscous truncation error associated with this scheme is controlled so that first order accuracy is maintained in time and space. The results encompass a range of Grashof numbers from 8.34 times 10,000 to 7 times 10 to the 7th power which is within the laminar flow regime for gravitationally driven fluid flows. Experiments within a small scale instrumented horizontal cylinder revealed the time development of the temperature distribution across the boundary layer and also the decay of wall heat transfer with time.

Singular boundary method for global gravity field modelling

NASA Astrophysics Data System (ADS)

Cunderlik, Robert

2014-05-01

The singular boundary method (SBM) and method of fundamental solutions (MFS) are meshless boundary collocation techniques that use the fundamental solution of a governing partial differential equation (e.g. the Laplace equation) as their basis functions. They have been developed to avoid singular numerical integration as well as mesh generation in the traditional boundary element method (BEM). SBM have been proposed to overcome a main drawback of MFS - its controversial fictitious boundary outside the domain. The key idea of SBM is to introduce a concept of the origin intensity factors that isolate singularities of the fundamental solution and its derivatives using some appropriate regularization techniques. Consequently, the source points can be placed directly on the real boundary and coincide with the collocation nodes. In this study we deal with SBM applied for high-resolution global gravity field modelling. The first numerical experiment presents a numerical solution to the fixed gravimetric boundary value problem. The achieved results are compared with the numerical solutions obtained by MFS or the direct BEM indicating efficiency of all methods. In the second numerical experiments, SBM is used to derive the geopotential and its first derivatives from the Tzz components of the gravity disturbing tensor observed by the GOCE satellite mission. A determination of the origin intensity factors allows to evaluate the disturbing potential and gravity disturbances directly on the Earth's surface where the source points are located. To achieve high-resolution numerical solutions, the large-scale parallel computations are performed on the cluster with 1TB of the distributed memory and an iterative elimination of far zones' contributions is applied.

Finite element solution of optimal control problems with inequality constraints

NASA Technical Reports Server (NTRS)

Bless, Robert R.; Hodges, Dewey H.

1990-01-01

A finite-element method based on a weak Hamiltonian form of the necessary conditions is summarized for optimal control problems. Very crude shape functions (so simple that element numerical quadrature is not necessary) can be used to develop an efficient procedure for obtaining candidate solutions (i.e., those which satisfy all the necessary conditions) even for highly nonlinear problems. An extension of the formulation allowing for discontinuities in the states and derivatives of the states is given. A theory that includes control inequality constraints is fully developed. An advanced launch vehicle (ALV) model is presented. The model involves staging and control constraints, thus demonstrating the full power of the weak formulation to date. Numerical results are presented along with total elapsed computer time required to obtain the results. The speed and accuracy in obtaining the results make this method a strong candidate for a real-time guidance algorithm.

A variational numerical method based on finite elements for the nonlinear solution characteristics of the periodically forced Chen system

NASA Astrophysics Data System (ADS)

Khan, Sabeel M.; Sunny, D. A.; Aqeel, M.

2017-09-01

Nonlinear dynamical systems and their solutions are very sensitive to initial conditions and therefore need to be approximated carefully. In this article, we present and analyze nonlinear solution characteristics of the periodically forced Chen system with the application of a variational method based on the concept of finite time-elements. Our approach is based on the discretization of physical time space into finite elements where each time-element is mapped to a natural time space. The solution of the system is then determined in natural time space using a set of suitable basis functions. The numerical algorithm is presented and implemented to compute and analyze nonlinear behavior at different time-step sizes. The obtained results show an excellent agreement with the classical RK-4 and RK-5 methods. The accuracy and convergence of the method is shown by comparing numerically computed results with the exact solution for a test problem. The presented method has shown a great potential in dealing with the solutions of nonlinear dynamical systems and thus can be utilized in delineating different features and characteristics of their solutions.

Planet-disc interactions with Discontinuous Galerkin Methods using GPUs

NASA Astrophysics Data System (ADS)

Velasco Romero, David A.; Veiga, Maria Han; Teyssier, Romain; Masset, Frédéric S.

2018-05-01

We present a two-dimensional Cartesian code based on high order discontinuous Galerkin methods, implemented to run in parallel over multiple GPUs. A simple planet-disc setup is used to compare the behaviour of our code against the behaviour found using the FARGO3D code with a polar mesh. We make use of the time dependence of the torque exerted by the disc on the planet as a mean to quantify the numerical viscosity of the code. We find that the numerical viscosity of the Keplerian flow can be as low as a few 10-8r2Ω, r and Ω being respectively the local orbital radius and frequency, for fifth order schemes and resolution of ˜10-2r. Although for a single disc problem a solution of low numerical viscosity can be obtained at lower computational cost with FARGO3D (which is nearly an order of magnitude faster than a fifth order method), discontinuous Galerkin methods appear promising to obtain solutions of low numerical viscosity in more complex situations where the flow cannot be captured on a polar or spherical mesh concentric with the disc.

Recovery of time-dependent volatility in option pricing model

NASA Astrophysics Data System (ADS)

Deng, Zui-Cha; Hon, Y. C.; Isakov, V.

2016-11-01

In this paper we investigate an inverse problem of determining the time-dependent volatility from observed market prices of options with different strikes. Due to the non linearity and sparsity of observations, an analytical solution to the problem is generally not available. Numerical approximation is also difficult to obtain using most of the existing numerical algorithms. Based on our recent theoretical results, we apply the linearisation technique to convert the problem into an inverse source problem from which recovery of the unknown volatility function can be achieved. Two kinds of strategies, namely, the integral equation method and the Landweber iterations, are adopted to obtain the stable numerical solution to the inverse problem. Both theoretical analysis and numerical examples confirm that the proposed approaches are effective. The work described in this paper was partially supported by a grant from the Research Grant Council of the Hong Kong Special Administrative Region (Project No. CityU 101112) and grants from the NNSF of China (Nos. 11261029, 11461039), and NSF grants DMS 10-08902 and 15-14886 and by Emylou Keith and Betty Dutcher Distinguished Professorship at the Wichita State University (USA).

A remark on fractional differential equation involving I-function

NASA Astrophysics Data System (ADS)

Mishra, Jyoti

2018-02-01

The present paper deals with the solution of the fractional differential equation using the Laplace transform operator and its corresponding properties in the fractional calculus; we derive an exact solution of a complex fractional differential equation involving a special function known as I-function. The analysis of the some fractional integral with two parameters is presented using the suggested Theorem 1. In addition, some very useful corollaries are established and their proofs presented in detail. Some obtained exact solutions are depicted to see the effect of each fractional order. Owing to the wider applicability of the I-function, we can conclude that, the obtained results in our work generalize numerous well-known results obtained by specializing the parameters.

Periodic, complexiton solutions and stability for a (2+1)-dimensional variable-coefficient Gross-Pitaevskii equation in the Bose-Einstein condensation

NASA Astrophysics Data System (ADS)

Yin, Hui-Min; Tian, Bo; Zhao, Xin-Chao

2018-06-01

This paper presents an investigation of a (2 + 1)-dimensional variable-coefficient Gross-Pitaevskii equation in the Bose-Einstein condensation. Periodic and complexiton solutions are obtained. Solitons solutions are also gotten through the periodic solutions. Numerical solutions via the split step method are stable. Effects of the weak and strong modulation instability on the solitons are shown: the weak modulation instability permits an observable soliton, and the strong one overwhelms its development.

Optimal impulsive manoeuvres and aerodynamic braking

NASA Technical Reports Server (NTRS)

Jezewski, D. J.

1985-01-01

A method developed for obtaining solutions to the aerodynamic braking problem, using impulses in the exoatmospheric phases is discussed. The solution combines primer vector theory and the results of a suboptimal atmospheric guidance program. For a specified initial and final orbit, the solution determines: (1) the minimum impulsive cost using a maximum of four impulses, (2) the optimal atmospheric entry and exit-state vectors subject to equality and inequality constraints, and (3) the optimal coast times. Numerical solutions which illustrate the characteristics of the solution are presented.

Magnetohydrodynamic viscous flow over a nonlinearly moving surface: Closed-form solutions

NASA Astrophysics Data System (ADS)

Fang, Tiegang

2014-05-01

In this paper, the magnetohydrodynamic (MHD) flow over a nonlinearly (power-law velocity) moving surface is investigated analytically and solutions are presented for a few special conditions. The solutions are obtained in closed forms with hyperbolic functions. The effects of the magnetic, the wall moving, and the mass transpiration parameters are discussed. These solutions are important to show the flow physics as well as to be used as bench mark problems for numerical validation and development of new solution schemes.

A new approach to exact optical soliton solutions for the nonlinear Schrödinger equation

NASA Astrophysics Data System (ADS)

Morales-Delgado, V. F.; Gómez-Aguilar, J. F.; Baleanu, Dumitru

2018-05-01

By using the modified homotopy analysis transform method, we construct the analytical solutions of the space-time generalized nonlinear Schrödinger equation involving a new fractional conformable derivative in the Liouville-Caputo sense and the fractional-order derivative with the Mittag-Leffler law. Employing theoretical parameters, we present some numerical simulations and compare the solutions obtained.

A new perspective for quintic B-spline based Crank-Nicolson-differential quadrature method algorithm for numerical solutions of the nonlinear Schrödinger equation

NASA Astrophysics Data System (ADS)

Başhan, Ali; Uçar, Yusuf; Murat Yağmurlu, N.; Esen, Alaattin

2018-01-01

In the present paper, a Crank-Nicolson-differential quadrature method (CN-DQM) based on utilizing quintic B-splines as a tool has been carried out to obtain the numerical solutions for the nonlinear Schrödinger (NLS) equation. For this purpose, first of all, the Schrödinger equation has been converted into coupled real value differential equations and then they have been discretized using both the forward difference formula and the Crank-Nicolson method. After that, Rubin and Graves linearization techniques have been utilized and the differential quadrature method has been applied to obtain an algebraic equation system. Next, in order to be able to test the efficiency of the newly applied method, the error norms, L2 and L_{∞}, as well as the two lowest invariants, I1 and I2, have been computed. Besides those, the relative changes in those invariants have been presented. Finally, the newly obtained numerical results have been compared with some of those available in the literature for similar parameters. This comparison clearly indicates that the currently utilized method, namely CN-DQM, is an effective and efficient numerical scheme and allows us to propose to solve a wide range of nonlinear equations.

A Numerical Study of Coupled Non-Linear Equations of Thermo-Viscous Fluid Flow in Cylindrical Geometry

NASA Astrophysics Data System (ADS)

Pothanna, N.; Aparna, P.; Gorla, R. S. R.

2017-12-01

In this paper we present numerical solutions to coupled non-linear governing equations of thermo-viscous fluid flow in cylindrical geometry using MATHEMATICA software solver. The numerical results are presented in terms of velocity, temperature and pressure distribution for various values of the material parameters such as the thermo-mechanical stress coefficient, thermal conductivity coefficient, Reiner Rivlin cross viscosity coefficient and the Prandtl number in the form of tables and graphs. Also, the solutions to governing equations for slow steady motion of a fluid have been obtained numerically and compared with the existing analytical results and are found to be in excellent agreement. The results of the present study will hopefully enable a better understanding applications of the flow under consideration.

A spectral tau algorithm based on Jacobi operational matrix for numerical solution of time fractional diffusion-wave equations

NASA Astrophysics Data System (ADS)

Bhrawy, A. H.; Doha, E. H.; Baleanu, D.; Ezz-Eldien, S. S.

2015-07-01

In this paper, an efficient and accurate spectral numerical method is presented for solving second-, fourth-order fractional diffusion-wave equations and fractional wave equations with damping. The proposed method is based on Jacobi tau spectral procedure together with the Jacobi operational matrix for fractional integrals, described in the Riemann-Liouville sense. The main characteristic behind this approach is to reduce such problems to those of solving systems of algebraic equations in the unknown expansion coefficients of the sought-for spectral approximations. The validity and effectiveness of the method are demonstrated by solving five numerical examples. Numerical examples are presented in the form of tables and graphs to make comparisons with the results obtained by other methods and with the exact solutions more easier.

Numerical solution of equations governing longitudinal suspension line wave motion during the parachute unfurling process. Ph.D. Thesis - George Washington Univ., Washington, D. C.

NASA Technical Reports Server (NTRS)

Poole, L. R.

1973-01-01

Equations are presented which govern the dynamics of the lines-first parachute unfurling process, including wave motion in the parachute suspension lines. Techniques are developed for obtaining numerical solutions to the governing equations. Histories of tension at test data, and generally good agreement is observed. Errors in computed results are attributed to several areas of uncertainty, the most significant being a poorly defined boundary condition on the wave motion at the vehicle-suspension line boundary.

Numerical simulations of the flow with the prescribed displacement of the airfoil and comparison with experiment

NASA Astrophysics Data System (ADS)

Řidký, V.; Šidlof, P.; Vlček, V.

2013-04-01

The work is devoted to comparing measured data with the results of numerical simulations. As mathematical model was used mathematical model whitout turbulence for incompressible flow In the experiment was observed the behavior of designed NACA0015 airfoil in airflow. For the numerical solution was used OpenFOAM computational package, this is open-source software based on finite volume method. In the numerical solution is prescribed displacement of the airfoil, which corresponds to the experiment. The velocity at a point close to the airfoil surface is compared with the experimental data obtained from interferographic measurements of the velocity field. Numerical solution is computed on a 3D mesh composed of about 1 million ortogonal hexahedron elements. The time step is limited by the Courant number. Parallel computations are run on supercomputers of the CIV at Technical University in Prague (HAL and FOX) and on a computer cluster of the Faculty of Mechatronics of Liberec (HYDRA). Run time is fixed at five periods, the results from the fifth periods and average value for all periods are then be compared with experiment.

Analytical approximation and numerical simulations for periodic travelling water waves

NASA Astrophysics Data System (ADS)

Kalimeris, Konstantinos

2017-12-01

We present recent analytical and numerical results for two-dimensional periodic travelling water waves with constant vorticity. The analytical approach is based on novel asymptotic expansions. We obtain numerical results in two different ways: the first is based on the solution of a constrained optimization problem, and the second is realized as a numerical continuation algorithm. Both methods are applied on some examples of non-constant vorticity. This article is part of the theme issue 'Nonlinear water waves'.

Numerical Hydrodynamics in Special Relativity.

PubMed

Martí, J M; Müller, E

1999-01-01

This review is concerned with a discussion of numerical methods for the solution of the equations of special relativistic hydrodynamics (SRHD). Particular emphasis is put on a comprehensive review of the application of high-resolution shock-capturing methods in SRHD. Results obtained with different numerical SRHD methods are compared, and two astrophysical applications of SRHD flows are discussed. An evaluation of the various numerical methods is given and future developments are analyzed. Supplementary material is available for this article at 10.12942/lrr-1999-3.

A solution to the Navier-Stokes equations based upon the Newton Kantorovich method

NASA Technical Reports Server (NTRS)

Davis, J. E.; Gabrielsen, R. E.; Mehta, U. B.

1977-01-01

An implicit finite difference scheme based on the Newton-Kantorovich technique was developed for the numerical solution of the nonsteady, incompressible, two-dimensional Navier-Stokes equations in conservation-law form. The algorithm was second-order-time accurate, noniterative with regard to the nonlinear terms in the vorticity transport equation except at the earliest few time steps, and spatially factored. Numerical results were obtained with the technique for a circular cylinder at Reynolds number 15. Results indicate that the technique is in excellent agreement with other numerical techniques for all geometries and Reynolds numbers investigated, and indicates a potential for significant reduction in computation time over current iterative techniques.

New distributed activation energy model: numerical solution and application to pyrolysis kinetics of some types of biomass.

PubMed

Cai, Junmeng; Liu, Ronghou

2008-05-01

In the present paper, a new distributed activation energy model has been developed, considering the reaction order and the dependence of frequency factor on temperature. The proposed DAEM cannot be solved directly in a closed from, thus a method was used to obtain the numerical solution of the new DAEM equation. Two numerical examples to illustrate the proposed method were presented. The traditional DAEM and new DAEM have been used to simulate the pyrolytic process of some types of biomass. The new DAEM fitted the experimental data much better than the traditional DAEM as the dependence of the frequency factor on temperature was taken into account.

Analysis of Piezoelectric Actuator for Vibration Control of Composite plate

NASA Astrophysics Data System (ADS)

Gomaa, Ahmed R.; Hai, Huang

2017-07-01

Vibration analysis is studied numerically in this paper for a simply supported composite plate subjected to external loadings. Vibrations are controlled by using piezoelectric patches. Finite element method (ANSYS) is used for obtaining finite element model of the smart plate structure, a layered composite plate is manufactured experimentally and tested to obtain the structure mechanical properties. Different piezoelectric patch areas and different applied gain voltage effects on vibration attenuation is studied. The numerical solution is compared with the experimental work, a good agreement achieved.

Electromagnetic field scattering by a triangular aperture.

PubMed

Harrison, R E; Hyman, E

1979-03-15

The multiple Laplace transform has been applied to analysis and computation of scattering by a double triangular aperture. Results are obtained which match far-field intensity distributions observed in experiments. Arbitrary polarization components, as well as in-phase and quadrature-phase components, may be determined, in the transform domain, as a continuous function of distance from near to far-field for any orientation, aperture, and transformable waveform. Numerical results are obtained by application of numerical multiple inversions of the fully transformed solution.

Surface electrical properties experiment. Part 2: Theory of radio-frequency interferometry in geophysical subsurface probing

NASA Technical Reports Server (NTRS)

Kong, J. A.; Tsang, L.

1974-01-01

The radiation fields due to a horizontal electric dipole laid on the surface of a stratified medium were calculated using a geometrical optics approximation, a modal approach, and direct numerical integration. The solutions were obtained from the reflection coefficient formulation and written in integral forms. The calculated interference patterns are compared in terms of the usefulness of the methods used to obtain them. Scattering effects are also discussed and all numerical results for anisotropic and isotropic cases are presented.

Radiation of sound from unflanged cylindrical ducts

NASA Technical Reports Server (NTRS)

Hartharan, S. L.; Bayliss, A.

1983-01-01

Calculations of sound radiated from unflanged cylindrical ducts are presented. The numerical simulation models the problem of an aero-engine inlet. The time dependent linearized Euler equations are solved from a state of rest until a harmonic solution is attained. A fourth order accurate finite difference scheme is used and solutions are obtained from a fully vectorized Cyber-203 computer program. Cases of both plane waves and spin modes are treated. Spin modes model the sound generated by a turbofan engine. Boundary conditions for both plane waves and spin modes are treated. Solutions obtained are compared with experiments conducted at NASA Langley Research Center.

Fully Nonlinear Modeling and Analysis of Precision Membranes

NASA Technical Reports Server (NTRS)

Pai, P. Frank; Young, Leyland G.

2003-01-01

High precision membranes are used in many current space applications. This paper presents a fully nonlinear membrane theory with forward and inverse analyses of high precision membrane structures. The fully nonlinear membrane theory is derived from Jaumann strains and stresses, exact coordinate transformations, the concept of local relative displacements, and orthogonal virtual rotations. In this theory, energy and Newtonian formulations are fully correlated, and every structural term can be interpreted in terms of vectors. Fully nonlinear ordinary differential equations (ODES) governing the large static deformations of known axisymmetric membranes under known axisymmetric loading (i.e., forward problems) are presented as first-order ODES, and a method for obtaining numerically exact solutions using the multiple shooting procedure is shown. A method for obtaining the undeformed geometry of any axisymmetric membrane with a known inflated geometry and a known internal pressure (i.e., inverse problems) is also derived. Numerical results from forward analysis are verified using results in the literature, and results from inverse analysis are verified using known exact solutions and solutions from the forward analysis. Results show that the membrane theory and the proposed numerical methods for solving nonlinear forward and inverse membrane problems are accurate.

Perturbation solutions of combustion instability problems

NASA Technical Reports Server (NTRS)

Googerdy, A.; Peddieson, J., Jr.; Ventrice, M.

1979-01-01

A method involving approximate modal analysis using the Galerkin method followed by an approximate solution of the resulting modal-amplitude equations by the two-variable perturbation method (method of multiple scales) is applied to two problems of pressure-sensitive nonlinear combustion instability in liquid-fuel rocket motors. One problem exhibits self-coupled instability while the other exhibits mode-coupled instability. In both cases it is possible to carry out the entire linear stability analysis and significant portions of the nonlinear stability analysis in closed form. In the problem of self-coupled instability the nonlinear stability boundary and approximate forms of the limit-cycle amplitudes and growth and decay rates are determined in closed form while the exact limit-cycle amplitudes and growth and decay rates are found numerically. In the problem of mode-coupled instability the limit-cycle amplitudes are found in closed form while the growth and decay rates are found numerically. The behavior of the solutions found by the perturbation method are in agreement with solutions obtained using complex numerical methods.

Numerical solution to the oblique derivative boundary value problem on non-uniform grids above the Earth topography

NASA Astrophysics Data System (ADS)

Medl'a, Matej; Mikula, Karol; Čunderlík, Róbert; Macák, Marek

2018-01-01

The paper presents a numerical solution of the oblique derivative boundary value problem on and above the Earth's topography using the finite volume method (FVM). It introduces a novel method for constructing non-uniform hexahedron 3D grids above the Earth's surface. It is based on an evolution of a surface, which approximates the Earth's topography, by mean curvature. To obtain optimal shapes of non-uniform 3D grid, the proposed evolution is accompanied by a tangential redistribution of grid nodes. Afterwards, the Laplace equation is discretized using FVM developed for such a non-uniform grid. The oblique derivative boundary condition is treated as a stationary advection equation, and we derive a new upwind type discretization suitable for non-uniform 3D grids. The discretization of the Laplace equation together with the discretization of the oblique derivative boundary condition leads to a linear system of equations. The solution of this system gives the disturbing potential in the whole computational domain including the Earth's surface. Numerical experiments aim to show properties and demonstrate efficiency of the developed FVM approach. The first experiments study an experimental order of convergence of the method. Then, a reconstruction of the harmonic function on the Earth's topography, which is generated from the EGM2008 or EIGEN-6C4 global geopotential model, is presented. The obtained FVM solutions show that refining of the computational grid leads to more precise results. The last experiment deals with local gravity field modelling in Slovakia using terrestrial gravity data. The GNSS-levelling test shows accuracy of the obtained local quasigeoid model.

Computer simulation of concentrated solid solution strengthening

NASA Technical Reports Server (NTRS)

Kuo, C. T. K.; Arsenault, R. J.

1976-01-01

The interaction forces between a straight edge dislocation moving through a three-dimensional block containing a random array of solute atoms were determined. The yield stress at 0 K was obtained by determining the average maximum solute-dislocation interaction force that is encountered by edge dislocation, and an expression relating the yield stress to the length of the dislocation and the solute concentration is provided. The magnitude of the solid solution strengthening due to solute atoms can be determined directly from the numerical results, provided the dislocation line length that moves as a unit is specified.

The New Method of Tsunami Source Reconstruction With r-Solution Inversion Method

NASA Astrophysics Data System (ADS)

Voronina, T. A.; Romanenko, A. A.

2016-12-01

Application of the r-solution method to reconstructing the initial tsunami waveform is discussed. This methodology is based on the inversion of remote measurements of water-level data. The wave propagation is considered within the scope of a linear shallow-water theory. The ill-posed inverse problem in question is regularized by means of a least square inversion using the truncated Singular Value Decomposition method. As a result of the numerical process, an r-solution is obtained. The method proposed allows one to control the instability of a numerical solution and to obtain an acceptable result in spite of ill posedness of the problem. Implementation of this methodology to reconstructing of the initial waveform to 2013 Solomon Islands tsunami validates the theoretical conclusion for synthetic data and a model tsunami source: the inversion result strongly depends on data noisiness, the azimuthal and temporal coverage of recording stations with respect to the source area. Furthermore, it is possible to make a preliminary selection of the most informative set of the available recording stations used in the inversion process.

Numerical Solution of the Extended Nernst-Planck Model.

PubMed

Samson; Marchand

1999-07-01

The main features of a numerical model aiming at predicting the drift of ions in an electrolytic solution upon a chemical potential gradient are presented. The mechanisms of ionic diffusion are described by solving the extended Nernst-Planck system of equations. The electrical coupling between the various ionic fluxes is accounted for by the Poisson equation. Furthermore, chemical activity effects are considered in the model. The whole system of nonlinear equations is solved using the finite-element method. Results yielded by the model for simple test cases are compared to those obtained using an analytical solution. Applications of the model to more complex problems are also presented and discussed. Copyright 1999 Academic Press.

Variational and numerical analysis of a quasistatic viscoelastic problem with normal compliance, friction and damage

NASA Astrophysics Data System (ADS)

Han, Weimin; Shillor, Meir; Sofonea, Mircea

2001-12-01

We consider a model for quasistatic frictional contact between a viscoelastic body and a foundation. The material constitutive relation is assumed to be nonlinear. The mechanical damage of the material, caused by excessive stress or strain, is described by the damage function, the evolution of which is determined by a parabolic inclusion. The contact is modeled with the normal compliance condition and the associated version of Coulomb's law of dry friction. We derive a variational formulation for the problem and prove the existence of its unique weak solution. We then study a fully discrete scheme for the numerical solutions of the problem and obtain error estimates on the approximate solutions.

Low-density lipoprotein transport through an arterial wall under hyperthermia and hypertension conditions--An analytical solution.

PubMed

Iasiello, Marcello; Vafai, Kambiz; Andreozzi, Assunta; Bianco, Nicola

2016-01-25

An analytical solution for Low-Density Lipoprotein transport through an arterial wall under hyperthermia conditions is established in this work. A four-layer model is used to characterize the arterial wall. Transport governing equations are obtained as a combination between Staverman-Kedem-Katchalsky membrane equations and volume-averaged porous media equations. Temperature and solute transport fields are coupled by means of Ludwig-Soret effect. Results are in excellent agreement with numerical and analytical literature data under isothermal conditions, and with numerical literature data for the hyperthermia case. Effects of hypertension combined with hyperthermia, are also analyzed in this work. Copyright © 2015 Elsevier Ltd. All rights reserved.

Neoclassical transport including collisional nonlinearity.

PubMed

Candy, J; Belli, E A

2011-06-10

In the standard δf theory of neoclassical transport, the zeroth-order (Maxwellian) solution is obtained analytically via the solution of a nonlinear equation. The first-order correction δf is subsequently computed as the solution of a linear, inhomogeneous equation that includes the linearized Fokker-Planck collision operator. This equation admits analytic solutions only in extreme asymptotic limits (banana, plateau, Pfirsch-Schlüter), and so must be solved numerically for realistic plasma parameters. Recently, numerical codes have appeared which attempt to compute the total distribution f more accurately than in the standard ordering by retaining some nonlinear terms related to finite-orbit width, while simultaneously reusing some form of the linearized collision operator. In this work we show that higher-order corrections to the distribution function may be unphysical if collisional nonlinearities are ignored.

On Chorin's Method for Stationary Solutions of the Oberbeck-Boussinesq Equation

NASA Astrophysics Data System (ADS)

Kagei, Yoshiyuki; Nishida, Takaaki

2017-06-01

Stability of stationary solutions of the Oberbeck-Boussinesq system (OB) and the corresponding artificial compressible system is considered. The latter system is obtained by adding the time derivative of the pressure with small parameter ɛ > 0 to the continuity equation of (OB), which was proposed by A. Chorin to find stationary solutions of (OB) numerically. Both systems have the same sets of stationary solutions and the system (OB) is obtained from the artificial compressible one as the limit ɛ \\to 0 which is a singular limit. It is proved that if a stationary solution of the artificial compressible system is stable for sufficiently small ɛ > 0, then it is also stable as a solution of (OB). The converse is proved provided that the velocity field of the stationary solution satisfies some smallness condition.

Numerical solution to the glancing sidewall oblique shock wave/turbulent boundary layer interaction in three dimension

NASA Technical Reports Server (NTRS)

Anderson, B. H.; Benson, T. J.

1983-01-01

A supersonic three-dimensional viscous forward-marching computer design code called PEPSIS is used to obtain a numerical solution of the three-dimensional problem of the interaction of a glancing sidewall oblique shock wave and a turbulent boundary layer. Very good results are obtained for a test case that was run to investigate the use of the wall-function boundary-condition approximation for a highly complex three-dimensional shock-boundary layer interaction. Two additional test cases (coarse mesh and medium mesh) are run to examine the question of near-wall resolution when no-slip boundary conditions are applied. A comparison with experimental data shows that the PEPSIS code gives excellent results in general and is practical for three-dimensional supersonic inlet calculations.

Numerical solution to the glancing sidewall oblique shock wave/turbulent boundary layer interaction in three-dimension

NASA Technical Reports Server (NTRS)

Anderson, B. H.; Benson, T. J.

1983-01-01

A supersonic three-dimensional viscous forward-marching computer design code called PEPSIS is used to obtain a numerical solution of the three-dimensional problem of the interaction of a glancing sidewall oblique shock wave and a turbulent boundary layer. Very good results are obtained for a test case that was run to investigate the use of the wall-function boundary-condition approximation for a highly complex three-dimensional shock-boundary layer interaction. Two additional test cases (coarse mesh and medium mesh) are run to examine the question of near-wall resolution when no-slip boundary conditions are applied. A comparison with experimental data shows that the PEPSIS code gives excellent results in general and is practical for three-dimensional supersonic inlet calculations.

Aeroacoustic Simulation of Nose Landing Gear on Adaptive Unstructured Grids With FUN3D

NASA Technical Reports Server (NTRS)

Vatsa, Veer N.; Khorrami, Mehdi R.; Park, Michael A.; Lockard, David P.

2013-01-01

Numerical simulations have been performed for a partially-dressed, cavity-closed nose landing gear configuration that was tested in NASA Langley s closed-wall Basic Aerodynamic Research Tunnel (BART) and in the University of Florida's open-jet acoustic facility known as the UFAFF. The unstructured-grid flow solver FUN3D, developed at NASA Langley Research center, is used to compute the unsteady flow field for this configuration. Starting with a coarse grid, a series of successively finer grids were generated using the adaptive gridding methodology available in the FUN3D code. A hybrid Reynolds-averaged Navier-Stokes/large eddy simulation (RANS/LES) turbulence model is used for these computations. Time-averaged and instantaneous solutions obtained on these grids are compared with the measured data. In general, the correlation with the experimental data improves with grid refinement. A similar trend is observed for sound pressure levels obtained by using these CFD solutions as input to a FfowcsWilliams-Hawkings noise propagation code to compute the farfield noise levels. In general, the numerical solutions obtained on adapted grids compare well with the hand-tuned enriched fine grid solutions and experimental data. In addition, the grid adaption strategy discussed here simplifies the grid generation process, and results in improved computational efficiency of CFD simulations.

Three-dimensional MHD (magnetohydrodynamic) flows in rectangular ducts of liquid-metal-cooled blankets

DOE Office of Scientific and Technical Information (OSTI.GOV)

Hua, T.Q.; Walker, J.S.; Picologlou, B.F.

1988-07-01

Magnetohydrodynamic flows of liquid metals in rectangular ducts with thin conducting walls in the presence of strong nonuniform transverse magnetic fields are examined. The interaction parameter and Hartmann number are assumed to be large, whereas the magnetic Reynolds number is assumed to be small. Under these assumptions, viscous and inertial effects are confined in very thin boundary layers adjacent to the walls. A significant fraction of the fluid flow is concentrated in the boundary layers adjacent to the side walls which are parallel to the magnetic field. This paper describes the analysis and numerical methods for obtaining 3-D solutions formore » flow parameters outside these layers, without solving explicitly for the layers themselves. Numerical solutions are presented for cases which are relevant to the flows of liquid metals in fusion reactor blankets. Experimental results obtained from the ALEX experiments at Argonne National Laboratory are used to validate the numerical code. In general, the agreement is excellent. 5 refs., 14 figs.« less

Boundary-layer effects in composite laminates: Free-edge stress singularities, part 6

NASA Technical Reports Server (NTRS)

Wanag, S. S.; Choi, I.

1981-01-01

A rigorous mathematical model was obtained for the boundary-layer free-edge stress singularity in angleplied and crossplied fiber composite laminates. The solution was obtained using a method consisting of complex-variable stress function potentials and eigenfunction expansions. The required order of the boundary-layer stress singularity is determined by solving the transcendental characteristic equation obtained from the homogeneous solution of the partial differential equations. Numerical results obtained show that the boundary-layer stress singularity depends only upon material elastic constants and fiber orientation of the adjacent plies. For angleplied and crossplied laminates the order of the singularity is weak in general.

Return Trajectory of the SpaceShipTwo Spacecraft--Numerical Solution

ERIC Educational Resources Information Center

Slegr, J.; Kraus, I.

2012-01-01

SpaceShipTwo is a private spaceplane project which is intended for space tourism. Very few details about its construction and flight characteristics are available for the public, but with proper numerical methods some interesting results can be obtained using secondary school mathematics. An exercise about SpaceShipTwo can be used as a…

N-person differential games. Part 2: The penalty method

NASA Technical Reports Server (NTRS)

Chen, G.; Mills, W. H.; Zheng, Q.; Shaw, W. H.

1983-01-01

The equilibrium strategy for N-person differential games can be found by studying a min-max problem subject to differential systems constraints. The differential constraints are penalized and finite elements are used to compute numerical solutions. Convergence proof and error estimates are given. Numerical results are also included and compared with those obtained by the dual method.

Numerical uncertainty in computational engineering and physics

DOE Office of Scientific and Technical Information (OSTI.GOV)

Hemez, Francois M

2009-01-01

Obtaining a solution that approximates ordinary or partial differential equations on a computational mesh or grid does not necessarily mean that the solution is accurate or even 'correct'. Unfortunately assessing the quality of discrete solutions by questioning the role played by spatial and temporal discretizations generally comes as a distant third to test-analysis comparison and model calibration. This publication is contributed to raise awareness of the fact that discrete solutions introduce numerical uncertainty. This uncertainty may, in some cases, overwhelm in complexity and magnitude other sources of uncertainty that include experimental variability, parametric uncertainty and modeling assumptions. The concepts ofmore » consistency, convergence and truncation error are overviewed to explain the articulation between the exact solution of continuous equations, the solution of modified equations and discrete solutions computed by a code. The current state-of-the-practice of code and solution verification activities is discussed. An example in the discipline of hydro-dynamics illustrates the significant effect that meshing can have on the quality of code predictions. A simple method is proposed to derive bounds of solution uncertainty in cases where the exact solution of the continuous equations, or its modified equations, is unknown. It is argued that numerical uncertainty originating from mesh discretization should always be quantified and accounted for in the overall uncertainty 'budget' that supports decision-making for applications in computational physics and engineering.« less

Slow Crack Growth Analysis of Brittle Materials with Finite Thickness Subjected to Constant Stress-Rate Flexural Loading

NASA Technical Reports Server (NTRS)

Chio, S. R.; Gyekenyesi, J. P.

1999-01-01

A two-dimensional, numerical analysis of slow crack growth (SCG) was performed for brittle materials with finite thickness subjected to constant stress-rate ("dynamic fatigue") loading in flexure. The numerical solution showed that the conventional, simple, one-dimensional analytical solution can be used with a maximum error of about 5% in determining the SCG parameters of a brittle material with the conditions of a normalized thickness (a ratio of specimen thickness to initial crack size) T > 3.3 and of a SCG parameter n > 10. The change in crack shape from semicircular to elliptical configurations was significant particularly at both low stress rate and low T, attributed to predominant difference in stress intensity factor along the crack front. The numerical solution of SCG parameters was supported within the experimental range by the data obtained from constant stress-rate flexural testing for soda-lime glass microslides at ambient temperature.

Fluid Flow and Solidification Under Combined Action of Magnetic Fields and Microgravity

NASA Technical Reports Server (NTRS)

Li, B. Q.; Shu, Y.; Li, K.; deGroh, H. C.

2002-01-01

Mathematical models, both 2-D and 3-D, are developed to represent g-jitter induced fluid flows and their effects on solidification under combined action of magnetic fields and microgravity. The numerical model development is based on the finite element solution of governing equations describing the transient g-jitter driven fluid flows, heat transfer and solutal transport during crystal growth with and without an applied magnetic field in space vehicles. To validate the model predictions, a ground-based g-jitter simulator is developed using the oscillating wall temperatures where timely oscillating fluid flows are measured using a laser PIV system. The measurements are compared well with numerical results obtained from the numerical models. Results show that a combined action derived from magnetic damping and microgravity can be an effective means to control the melt flow and solutal transport in space single crystal growth systems.

Finite element computation of a viscous compressible free shear flow governed by the time dependent Navier-Stokes equations

NASA Technical Reports Server (NTRS)

Cooke, C. H.; Blanchard, D. K.

1975-01-01

A finite element algorithm for solution of fluid flow problems characterized by the two-dimensional compressible Navier-Stokes equations was developed. The program is intended for viscous compressible high speed flow; hence, primitive variables are utilized. The physical solution was approximated by trial functions which at a fixed time are piecewise cubic on triangular elements. The Galerkin technique was employed to determine the finite-element model equations. A leapfrog time integration is used for marching asymptotically from initial to steady state, with iterated integrals evaluated by numerical quadratures. The nonsymmetric linear systems of equations governing time transition from step-to-step are solved using a rather economical block iterative triangular decomposition scheme. The concept was applied to the numerical computation of a free shear flow. Numerical results of the finite-element method are in excellent agreement with those obtained from a finite difference solution of the same problem.

An adaptive gridless methodology in one dimension

DOE Office of Scientific and Technical Information (OSTI.GOV)

Snyder, N.T.; Hailey, C.E.

1996-09-01

Gridless numerical analysis offers great potential for accurately solving for flow about complex geometries or moving boundary problems. Because gridless methods do not require point connection, the mesh cannot twist or distort. The gridless method utilizes a Taylor series about each point to obtain the unknown derivative terms from the current field variable estimates. The governing equation is then numerically integrated to determine the field variables for the next iteration. Effects of point spacing and Taylor series order on accuracy are studied, and they follow similar trends of traditional numerical techniques. Introducing adaption by point movement using a spring analogymore » allows the solution method to track a moving boundary. The adaptive gridless method models linear, nonlinear, steady, and transient problems. Comparison with known analytic solutions is given for these examples. Although point movement adaption does not provide a significant increase in accuracy, it helps capture important features and provides an improved solution.« less

Adiabatic model of field reversal by fast ions in an axisymmetric open trap

DOE Office of Scientific and Technical Information (OSTI.GOV)

Tsidulko, Yu. A., E-mail: tsidulko@mail.ru

2016-06-15

A model of field reversal by fast ions has been developed under the assumption of preservation of fast-ion adiabatic invariants. Analytical solutions obtained in the approximation of a narrow fast-ion layer and numerical solutions to the evolutionary problem are presented. The solutions demonstrate the process of formation of a field reversed configuration with parameters close to those of the planned experiment.

Localized solutions of Lugiato-Lefever equations with focused pump.

PubMed

Cardoso, Wesley B; Salasnich, Luca; Malomed, Boris A

2017-12-04

Lugiato-Lefever (LL) equations in one and two dimensions (1D and 2D) accurately describe the dynamics of optical fields in pumped lossy cavities with the intrinsic Kerr nonlinearity. The external pump is usually assumed to be uniform, but it can be made tightly focused too-in particular, for building small pixels. We obtain solutions of the LL equations, with both the focusing and defocusing intrinsic nonlinearity, for 1D and 2D confined modes supported by the localized pump. In the 1D setting, we first develop a simple perturbation theory, based in the sech ansatz, in the case of weak pump and loss. Then, a family of exact analytical solutions for spatially confined modes is produced for the pump focused in the form of a delta-function, with a nonlinear loss (two-photon absorption) added to the LL model. Numerical findings demonstrate that these exact solutions are stable, both dynamically and structurally (the latter means that stable numerical solutions close to the exact ones are found when a specific condition, necessary for the existence of the analytical solution, does not hold). In 2D, vast families of stable confined modes are produced by means of a variational approximation and full numerical simulations.

Analytical and numerical solutions for heat transfer and effective thermal conductivity of cracked media

NASA Astrophysics Data System (ADS)

Tran, A. B.; Vu, M. N.; Nguyen, S. T.; Dong, T. Q.; Le-Nguyen, K.

2018-02-01

This paper presents analytical solutions to heat transfer problems around a crack and derive an adaptive model for effective thermal conductivity of cracked materials based on singular integral equation approach. Potential solution of heat diffusion through two-dimensional cracked media, where crack filled by air behaves as insulator to heat flow, is obtained in a singular integral equation form. It is demonstrated that the temperature field can be described as a function of temperature and rate of heat flow on the boundary and the temperature jump across the cracks. Numerical resolution of this boundary integral equation allows determining heat conduction and effective thermal conductivity of cracked media. Moreover, writing this boundary integral equation for an infinite medium embedding a single crack under a far-field condition allows deriving the closed-form solution of temperature discontinuity on the crack and particularly the closed-form solution of temperature field around the crack. These formulas are then used to establish analytical effective medium estimates. Finally, the comparison between the developed numerical and analytical solutions allows developing an adaptive model for effective thermal conductivity of cracked media. This model takes into account both the interaction between cracks and the percolation threshold.

Approximate bound-state solutions of the Dirac equation for the generalized yukawa potential plus the generalized tensor interaction

NASA Astrophysics Data System (ADS)

Ikot, Akpan N.; Maghsoodi, Elham; Hassanabadi, Hassan; Obu, Joseph A.

2014-05-01

In this paper, we obtain the approximate analytical bound-state solutions of the Dirac particle with the generalized Yukawa potential within the framework of spin and pseudospin symmetries for the arbitrary к state with a generalized tensor interaction. The generalized parametric Nikiforov-Uvarov method is used to obtain the energy eigenvalues and the corresponding wave functions in closed form. We also report some numerical results and present figures to show the effect of the tensor interaction.

Numerical solution of linear and nonlinear Fredholm integral equations by using weighted mean-value theorem.

PubMed

Altürk, Ahmet

2016-01-01

Mean value theorems for both derivatives and integrals are very useful tools in mathematics. They can be used to obtain very important inequalities and to prove basic theorems of mathematical analysis. In this article, a semi-analytical method that is based on weighted mean-value theorem for obtaining solutions for a wide class of Fredholm integral equations of the second kind is introduced. Illustrative examples are provided to show the significant advantage of the proposed method over some existing techniques.

Numerical solution of problems concerning the thermal convection of a variable-viscosity liquid

NASA Astrophysics Data System (ADS)

Zherebiatev, I. F.; Lukianov, A. T.; Podkopaev, Iu. L.

A stabilizing-correction scheme is constructed for integrating the fourth-order equation describing the dynamics of a viscous incompressible liquid. As an example, a solution is obtained to the problem of the solidification of a liquid in a rectangular region with allowance for convective energy transfer in the liquid phase as well as temperature-dependent changes of viscosity. It is noted that the proposed method can be used to study steady-state problems of thermal convection in ingots obtained through continuous casting.

The method of lines in analyzing solids containing cracks

NASA Technical Reports Server (NTRS)

Gyekenyesi, John P.

1990-01-01

A semi-numerical method is reviewed for solving a set of coupled partial differential equations subject to mixed and possibly coupled boundary conditions. The line method of analysis is applied to the Navier-Cauchy equations of elastic and elastoplastic equilibrium to calculate the displacement distributions in various, simple geometry bodies containing cracks. The application of this method to the appropriate field equations leads to coupled sets of simultaneous ordinary differential equations whose solutions are obtained along sets of lines in a discretized region. When decoupling of the equations and their boundary conditions is not possible, the use of a successive approximation procedure permits the analytical solution of the resulting ordinary differential equations. The use of this method is illustrated by reviewing and presenting selected solutions of mixed boundary value problems in three dimensional fracture mechanics. These solutions are of great importance in fracture toughness testing, where accurate stress and displacement distributions are required for the calculation of certain fracture parameters. Computations obtained for typical flawed specimens include that for elastic as well as elastoplastic response. Problems in both Cartesian and cylindrical coordinate systems are included. Results are summarized for a finite geometry rectangular bar with a central through-the-thickness or rectangular surface crack under remote uniaxial tension. In addition, stress and displacement distributions are reviewed for finite circular bars with embedded penny-shaped cracks, and rods with external annular or ring cracks under opening mode tension. The results obtained show that the method of lines presents a systematic approach to the solution of some three-dimensional mechanics problems with arbitrary boundary conditions. The advantage of this method over other numerical solutions is that good results are obtained even from the use of a relatively coarse grid.

The exact solution of shear-lag problems in flat panels and box beams assumed rigid in the transverse direction

NASA Technical Reports Server (NTRS)

Hildebrand, Francis B

1943-01-01

A mathematical procedure is herein developed for obtaining exact solutions of shear-lag problems in flat panels and box beams: the method is based on the assumption that the amount of stretching of the sheets in the direction perpendicular to the direction of essential normal stresses is negligible. Explicit solutions, including the treatment of cut-outs, are given for several cases and numerical results are presented in graphic and tabular form. The general theory is presented in a from which further solutions can be readily obtained. The extension of the theory to cover certain cases of non-uniform cross section is indicated. Although the solutions are obtained in terms of infinite series, the present developments differ from those previously given in that, in practical cases, the series usually converge so rapidly that sufficient accuracy is afforded by a small number of terms. Comparisons are made in several cases between the present results and the corresponding solutions obtained by approximate procedures devised by Reissner and by Kuhn and Chiarito.

The precise time-dependent solution of the Fokker–Planck equation with anomalous diffusion

DOE Office of Scientific and Technical Information (OSTI.GOV)

Guo, Ran; Du, Jiulin, E-mail: jiulindu@aliyun.com

2015-08-15

We study the time behavior of the Fokker–Planck equation in Zwanzig’s rule (the backward-Ito’s rule) based on the Langevin equation of Brownian motion with an anomalous diffusion in a complex medium. The diffusion coefficient is a function in momentum space and follows a generalized fluctuation–dissipation relation. We obtain the precise time-dependent analytical solution of the Fokker–Planck equation and at long time the solution approaches to a stationary power-law distribution in nonextensive statistics. As a test, numerically we have demonstrated the accuracy and validity of the time-dependent solution. - Highlights: • The precise time-dependent solution of the Fokker–Planck equation with anomalousmore » diffusion is found. • The anomalous diffusion satisfies a generalized fluctuation–dissipation relation. • At long time the time-dependent solution approaches to a power-law distribution in nonextensive statistics. • Numerically we have demonstrated the accuracy and validity of the time-dependent solution.« less

A scheme to calculate higher-order homogenization as applied to micro-acoustic boundary value problems

NASA Astrophysics Data System (ADS)

Vagh, Hardik A.; Baghai-Wadji, Alireza

2008-12-01

Current technological challenges in materials science and high-tech device industry require the solution of boundary value problems (BVPs) involving regions of various scales, e.g. multiple thin layers, fibre-reinforced composites, and nano/micro pores. In most cases straightforward application of standard variational techniques to BVPs of practical relevance necessarily leads to unsatisfactorily ill-conditioned analytical and/or numerical results. To remedy the computational challenges associated with sub-sectional heterogeneities various sophisticated homogenization techniques need to be employed. Homogenization refers to the systematic process of smoothing out the sub-structural heterogeneities, leading to the determination of effective constitutive coefficients. Ordinarily, homogenization involves a sophisticated averaging and asymptotic order analysis to obtain solutions. In the majority of the cases only zero-order terms are constructed due to the complexity of the processes involved. In this paper we propose a constructive scheme for obtaining homogenized solutions involving higher order terms, and thus, guaranteeing higher accuracy and greater robustness of the numerical results. We present

Analytical and numerical solutions for mass diffusion in a composite cylindrical body

NASA Astrophysics Data System (ADS)

Kumar, A.

1980-12-01

The analytical and numerical solution techniques were investigated to study moisture diffusion problems in cylindrical bodies that are assumed to be composed of a finite number of layers of different materials. A generalized diffusion model for an n-layer cylindrical body with discontinuous moisture content at the interfaces was developed and the formal solutions were obtained. The model is to be used for describing mass transfer rates of any composite body, such as an ear of corn which could be assumed of consisting two different layers: the inner core represents the woody cob and the outer cylinder represents the kernel layer. Data describing the fully exposed drying characteristics of ear corn at high air velocity were obtained under different drying conditions. Ear corns were modeled as homogeneous bodies since composite model did not improve the fit substantially. A computer program using multidimensional optimization technique showed that diffusivity was an exponential function of moisture content and an arrhenius function of temperature of drying air.

Magnetic Field in a Screw Flow with Fluctuations

NASA Astrophysics Data System (ADS)

Titov, V. V.; Stepanov, R. A.; Sokoloff, D. D.

2018-04-01

We consider the influence of fluctuations in a screw flow of a conducting liquid on the effect of magnetic field self-excitation; the solution of this problem is important for experimental realization of a turbulent dynamo. We propose a theoretical approach based on the solution of averaged equations obtained in the limit of a short correlation time. The applicability of this approach is confirmed by direct numerical simulation of the initial equations. We demonstrate the influence of the correlation of fluctuations on the dynamo effect threshold. It is shown that the solution of the mean-field equations differs from the solution based on direct numerical simulation for a finite correlation time. The advantages and disadvantages of the two approaches are estimates, as well as the importance of the discovered difference in the context of problems of magnetic field self-excitation. The influence of helicity and intermittency on the type of the solution is considered.

Relative motion of orbiting satellites

NASA Technical Reports Server (NTRS)

Eades, J. B., Jr.

1972-01-01

The relative motion problem is analyzed, as a linearized case, and as a numerically determined solution to provide a time history of the geometries representing the motion state. The displacement history and the hodographs for families of solutions are provided, analytically and graphically, to serve as an aid to understanding this problem area. Linearized solutions to relative motion problems of orbiting particles are presented for the impulsive and fixed thrust cases. Second order solutions are described to enhance the accuracy of prediction. A method was developed to obtain accurate, numerical solutions to the intercept and rendezvous problem; and, special situations are examined. A particular problem related to relative motions, where the motion traces develop a cusp, is examined in detail. This phenomenon is found to be dependent on a particular relationship between orbital eccentricity and the inclination between orbital planes. These conditions are determined, and, example situations are presented and discussed.

Numerical solutions for heat flow in adhesive lap joints

NASA Technical Reports Server (NTRS)

Howell, P. A.; Winfree, William P.

1992-01-01

The present formulation for the modeling of heat transfer in thin, adhesively bonded lap joints precludes difficulties associated with large aspect ratio grids required by standard FEM formulations. This quasi-static formulation also reduces the problem dimensionality (by one), thereby minimizing computational requirements. The solutions obtained are found to be in good agreement with both analytical solutions and solutions from standard FEM programs. The approach is noted to yield a more accurate representation of heat-flux changes between layers due to a disbond.

Two-dimensional solitary waves and periodic waves on coupled nonlinear electrical transmission lines

NASA Astrophysics Data System (ADS)

Wang, Heng; Zheng, Shuhua

2017-06-01

By using the dynamical system approach, the exact travelling wave solutions for a system of coupled nonlinear electrical transmission lines are studied. Based on this method, the bifurcations of phase portraits of a dynamical system are given. The two-dimensional solitary wave solutions and periodic wave solutions on coupled nonlinear transmission lines are obtained. With the aid of Maple, the numerical simulations are conducted for solitary wave solutions and periodic wave solutions to the model equation. The results presented in this paper improve upon previous studies.

Using Predictor-Corrector Methods in Numerical Solutions to Mathematical Problems of Motion

ERIC Educational Resources Information Center

Lewis, Jerome

2005-01-01

In this paper, the author looks at some classic problems in mathematics that involve motion in the plane. Many case problems like these are difficult and beyond the mathematical skills of most undergraduates, but computational approaches often require less insight into the subtleties of the problems and can be used to obtain reliable solutions.…

Upscaling of Solute Transport in Heterogeneous Media with Non-uniform Flow and Dispersion Fields

DOE Office of Scientific and Technical Information (OSTI.GOV)

Xu, Zhijie; Meakin, Paul

2013-10-01

An analytical and computational model for non-reactive solute transport in periodic heterogeneous media with arbitrary non-uniform flow and dispersion fields within the unit cell of length ε is described. The model lumps the effect of non-uniform flow and dispersion into an effective advection velocity Ve and an effective dispersion coefficient De. It is shown that both Ve and De are scale-dependent (dependent on the length scale of the microscopic heterogeneity, ε), dependent on the Péclet number Pe, and on a dimensionless parameter α that represents the effects of microscopic heterogeneity. The parameter α, confined to the range of [-0.5, 0.5]more » for the numerical example presented, depends on the flow direction and non-uniform flow and dispersion fields. Effective advection velocity Ve and dispersion coefficient De can be derived for any given flow and dispersion fields, and . Homogenized solutions describing the macroscopic variations can be obtained from the effective model. Solutions with sub-unit-cell accuracy can be constructed by homogenized solutions and its spatial derivatives. A numerical implementation of the model compared with direct numerical solutions using a fine grid, demonstrated that the new method was in good agreement with direct solutions, but with significant computational savings.« less

Dispersive models describing mosquitoes’ population dynamics

NASA Astrophysics Data System (ADS)

Yamashita, W. M. S.; Takahashi, L. T.; Chapiro, G.

2016-08-01

The global incidences of dengue and, more recently, zica virus have increased the interest in studying and understanding the mosquito population dynamics. Understanding this dynamics is important for public health in countries where climatic and environmental conditions are favorable for the propagation of these diseases. This work is based on the study of nonlinear mathematical models dealing with the life cycle of the dengue mosquito using partial differential equations. We investigate the existence of traveling wave solutions using semi-analytical method combining dynamical systems techniques and numerical integration. Obtained solutions are validated through numerical simulations using finite difference schemes.

Accurate boundary conditions for exterior problems in gas dynamics

NASA Technical Reports Server (NTRS)

Hagstrom, Thomas; Hariharan, S. I.

1988-01-01

The numerical solution of exterior problems is typically accomplished by introducing an artificial, far field boundary and solving the equations on a truncated domain. For hyperbolic systems, boundary conditions at this boundary are often derived by imposing a principle of no reflection. However, waves with spherical symmetry in gas dynamics satisfy equations where incoming and outgoing Riemann variables are coupled. This suggests that natural reflections may be important. A reflecting boundary condition is proposed based on an asymptotic solution of the far field equations. Nonlinear energy estimates are obtained for the truncated problem and numerical experiments presented to validate the theory.

Accurate boundary conditions for exterior problems in gas dynamics

NASA Technical Reports Server (NTRS)

Hagstrom, Thomas; Hariharan, S. I.

1988-01-01

The numerical solution of exterior problems is typically accomplished by introducing an artificial, far-field boundary and solving the equations on a truncated domain. For hyperbolic systems, boundary conditions at this boundary are often derived by imposing a principle of no reflection. However, waves with spherical symmetry in gas dynamics satisfy equations where incoming and outgoing Riemann variables are coupled. This suggests that natural reflections may be important. A reflecting boundary condition is proposed based on an asymptotic solution of the far-field equations. Nonlinear energy estimates are obtained for the truncated problem and numerical experiments presented to validate the theory.

A network thermodynamic method for numerical solution of the Nernst-Planck and Poisson equation system with application to ionic transport through membranes.

PubMed

Horno, J; González-Caballero, F; González-Fernández, C F

1990-01-01

Simple techniques of network thermodynamics are used to obtain the numerical solution of the Nernst-Planck and Poisson equation system. A network model for a particular physical situation, namely ionic transport through a thin membrane with simultaneous diffusion, convection and electric current, is proposed. Concentration and electric field profiles across the membrane, as well as diffusion potential, have been simulated using the electric circuit simulation program, SPICE. The method is quite general and extremely efficient, permitting treatments of multi-ion systems whatever the boundary and experimental conditions may be.

Numeric Solutions of Dirac-Gursey Spinor Field Equation Under External Gaussian White Noise

NASA Astrophysics Data System (ADS)

Aydogmus, Fatma

2016-06-01

In this paper, we consider the Dirac-Gursey spinor field equation that has particle-like solutions derived classical field equations so-called instantons, formed by using Heisenberg ansatz, under the effect of an additional Gaussian white noise term. Our purpose is to understand how the behavior of spinor-type excited instantons in four dimensions can be affected by noise. Thus, we simulate the phase portraits and Poincaré sections of the obtained system numerically both with and without noise. Recurrence plots are also given for more detailed information regarding the system.

Fourth order scheme for wavelet based solution of Black-Scholes equation

NASA Astrophysics Data System (ADS)

Finěk, Václav

2017-12-01

The present paper is devoted to the numerical solution of the Black-Scholes equation for pricing European options. We apply the Crank-Nicolson scheme with Richardson extrapolation for time discretization and Hermite cubic spline wavelets with four vanishing moments for space discretization. This scheme is the fourth order accurate both in time and in space. Computational results indicate that the Crank-Nicolson scheme with Richardson extrapolation significantly decreases the amount of computational work. We also numerically show that optimal convergence rate for the used scheme is obtained without using startup procedure despite the data irregularities in the model.

Fuzzy multi objective transportation problem – evolutionary algorithm approach

NASA Astrophysics Data System (ADS)

Karthy, T.; Ganesan, K.

2018-04-01

This paper deals with fuzzy multi objective transportation problem. An fuzzy optimal compromise solution is obtained by using Fuzzy Genetic Algorithm. A numerical example is provided to illustrate the methodology.

Some Aspects of Nonlinear Dynamics and CFD

NASA Technical Reports Server (NTRS)

Yee, Helen C.; Merriam, Marshal (Technical Monitor)

1996-01-01

The application of nonlinear dynamics to improve the understanding of numerical uncertainties in computational fluid dynamics (CFD) is reviewed. Elementary examples in the use of dynamics to explain the nonlinear phenomena and spurious behavior that occur in numerics are given. The role of dynamics in the understanding of long time behavior of numerical integrations and the nonlinear stability, convergence, and reliability of using time-marching approaches for obtaining steady-state numerical solutions in CFD is explained. The study is complemented with examples of spurious behavior observed in CFD computations.

Finite element analysis of low speed viscous and inviscid aerodynamic flows

NASA Technical Reports Server (NTRS)

Baker, A. J.; Manhardt, P. D.

1977-01-01

A weak interaction solution algorithm was established for aerodynamic flow about an isolated airfoil. Finite element numerical methodology was applied to solution of each of differential equations governing potential flow, and viscous and turbulent boundary layer and wake flow downstream of the sharp trailing edge. The algorithm accounts for computed viscous displacement effects on the potential flow. Closure for turbulence was accomplished using both first and second order models. The COMOC finite element fluid mechanics computer program was modified to solve the identified equation systems for two dimensional flows. A numerical program was completed to determine factors affecting solution accuracy, convergence and stability for the combined potential, boundary layer, and parabolic Navier-Stokes equation systems. Good accuracy and convergence are demonstrated. Each solution is obtained within the identical finite element framework of COMOC.

A computing method for sound propagation through a nonuniform jet stream

NASA Technical Reports Server (NTRS)

Padula, S. L.; Liu, C. H.

1974-01-01

The classical formulation of sound propagation through a jet flow was found to be inadequate for computer solutions. Previous investigations selected the phase and amplitude of the acoustic pressure as dependent variables requiring the solution of a system of nonlinear algebraic equations. The nonlinearities complicated both the analysis and the computation. A reformulation of the convective wave equation in terms of a new set of dependent variables is developed with a special emphasis on its suitability for numerical solutions on fast computers. The technique is very attractive because the resulting equations are linear in nonwaving variables. The computer solution to such a linear system of algebraic equations may be obtained by well-defined and direct means which are conservative of computer time and storage space. Typical examples are illustrated and computational results are compared with available numerical and experimental data.

Numerical solution of supersonic three-dimensional free-mixing flows using the parabolic-elliptic Navier-Stokes equations

NASA Technical Reports Server (NTRS)

Hirsh, R. S.

1976-01-01

A numerical method is presented for solving the parabolic-elliptic Navier-Stokes equations. The solution procedure is applied to three-dimensional supersonic laminar jet flow issuing parallel with a supersonic free stream. A coordinate transformation is introduced which maps the boundaries at infinity into a finite computational domain in order to eliminate difficulties associated with the imposition of free-stream boundary conditions. Results are presented for an approximate circular jet, a square jet, varying aspect ratio rectangular jets, and interacting square jets. The solution behavior varies from axisymmetric to nearly two-dimensional in character. For cases where comparisons of the present results with those obtained from shear layer calculations could be made, agreement was good.

Dynamical analysis of the avian-human influenza epidemic model using the semi-analytical method

NASA Astrophysics Data System (ADS)

Jabbari, Azizeh; Kheiri, Hossein; Bekir, Ahmet

2015-03-01

In this work, we present a dynamic behavior of the avian-human influenza epidemic model by using efficient computational algorithm, namely the multistage differential transform method(MsDTM). The MsDTM is used here as an algorithm for approximating the solutions of the avian-human influenza epidemic model in a sequence of time intervals. In order to show the efficiency of the method, the obtained numerical results are compared with the fourth-order Runge-Kutta method (RK4M) and differential transform method(DTM) solutions. It is shown that the MsDTM has the advantage of giving an analytical form of the solution within each time interval which is not possible in purely numerical techniques like RK4M.

Numerical solution of Space Shuttle Orbiter flow field including real gas effects

NASA Technical Reports Server (NTRS)

Prabhu, D. K.; Tannehill, J. C.

1984-01-01

The hypersonic, laminar flow around the Space Shuttle Orbiter has been computed for both an ideal gas (gamma = 1.2) and equilibrium air using a real-gas, parabolized Navier-Stokes code. This code employs a generalized coordinate transformation; hence, it places no restrictions on the orientation of the solution surfaces. The initial solution in the nose region was computed using a 3-D, real-gas, time-dependent Navier-Stokes code. The thermodynamic and transport properties of equilibrium air were obtained from either approximate curve fits or a table look-up procedure. Numerical results are presented for flight conditions corresponding to the STS-3 trajectory. The computed surface pressures and convective heating rates are compared with data from the STS-3 flight.

Methodology of Numerical Optimization for Orbital Parameters of Binary Systems

NASA Astrophysics Data System (ADS)

Araya, I.; Curé, M.

2010-02-01

The use of a numerical method of maximization (or minimization) in optimization processes allows us to obtain a great amount of solutions. Therefore, we can find a global maximum or minimum of the problem, but this is only possible if we used a suitable methodology. To obtain the global optimum values, we use the genetic algorithm called PIKAIA (P. Charbonneau) and other four algorithms implemented in Mathematica. We demonstrate that derived orbital parameters of binary systems published in some papers, based on radial velocity measurements, are local minimum instead of global ones.

Interior point techniques for LP and NLP

DOE Office of Scientific and Technical Information (OSTI.GOV)

Evtushenko, Y.

By using surjective mapping the initial constrained optimization problem is transformed to a problem in a new space with only equality constraints. For the numerical solution of the latter problem we use the generalized gradient-projection method and Newton`s method. After inverse transformation to the initial space we obtain the family of numerical methods for solving optimization problems with equality and inequality constraints. In the linear programming case after some simplification we obtain Dikin`s algorithm, affine scaling algorithm and generalized primal dual interior point linear programming algorithm.

Asymptotic traveling wave solution for a credit rating migration problem

NASA Astrophysics Data System (ADS)

Liang, Jin; Wu, Yuan; Hu, Bei

2016-07-01

In this paper, an asymptotic traveling wave solution of a free boundary model for pricing a corporate bond with credit rating migration risk is studied. This is the first study to associate the asymptotic traveling wave solution to the credit rating migration problem. The pricing problem with credit rating migration risk is modeled by a free boundary problem. The existence, uniqueness and regularity of the solution are obtained. Under some condition, we proved that the solution of our credit rating problem is convergent to a traveling wave solution, which has an explicit form. Furthermore, numerical examples are presented.

Random element method for numerical modeling of diffusional processes

NASA Technical Reports Server (NTRS)

Ghoniem, A. F.; Oppenheim, A. K.

1982-01-01

The random element method is a generalization of the random vortex method that was developed for the numerical modeling of momentum transport processes as expressed in terms of the Navier-Stokes equations. The method is based on the concept that random walk, as exemplified by Brownian motion, is the stochastic manifestation of diffusional processes. The algorithm based on this method is grid-free and does not require the diffusion equation to be discritized over a mesh, it is thus devoid of numerical diffusion associated with finite difference methods. Moreover, the algorithm is self-adaptive in space and explicit in time, resulting in an improved numerical resolution of gradients as well as a simple and efficient computational procedure. The method is applied here to an assortment of problems of diffusion of momentum and energy in one-dimension as well as heat conduction in two-dimensions in order to assess its validity and accuracy. The numerical solutions obtained are found to be in good agreement with exact solution except for a statistical error introduced by using a finite number of elements, the error can be reduced by increasing the number of elements or by using ensemble averaging over a number of solutions.

Transient natural convection with density inversion from a horizontal cylinder

NASA Astrophysics Data System (ADS)

Wang, P.; Kahawita, R.; Nguyen, D. L.

1992-01-01

This paper is devoted to a numerical investigation of the free convection flow about a horizontal cylinder maintained at 0 °C in a water ambient close to the point of maximum density. Complete numerical solutions covering both the transient as well as steady state have been obtained. Principal results indicate that the proximity of the ambient temperature to the point of maximum density plays an important role in the type of convection pattern that may be obtained. When the ambient temperature is within 4.7 °C

Numerical modelling of physical processes in a ballistic laboratory setup with a tapered adapter and plastic piston used for obtaining high muzzle velocities

NASA Astrophysics Data System (ADS)

Bykov, N. V.

2014-12-01

Numerical modelling of a ballistic setup with a tapered adapter and plastic piston is considered. The processes in the firing chamber are described within the framework of quasi- one-dimensional gas dynamics and a geometrical law of propellant burn by means of Lagrangian mass coordinates. The deformable piston is considered to be an ideal liquid with specific equations of state. The numerical solution is obtained by means of a modified explicit von Neumann scheme. The calculation results given show that the ballistic setup with a tapered adapter and plastic piston produces increased shell muzzle velocities by a factor of more than 1.5-2.

Assessment Processes to Increase the Burden of Existing Buildings Using BIM

NASA Astrophysics Data System (ADS)

Szeląg, Romuald

2017-10-01

The process of implementation of the reconstruction of buildings is often associated with the need to adapt them to increased loads. In the restricted access to the archive project documentation it is necessary to use technical solutions to obtain a fairly short period of time necessary to implement the technical parameters of such processes. Dissemination of BIM in the design process can also be used effectively in the processes of identification of existing facilities for the implementation of the work of strengthening or adapting objects to the increased load requirements. Obtained in the process of research and macroscopic data is then used in the processes of numerical processing aimed at developing a numerical model reflects the actual parameters of the structure of the existing structure and, therefore, allows a better look at the object and the execution of the process to strengthen future. This article will identify possibilities for the use of BIM in processes of identification technology buildings and structures and indicated the necessary data to be obtained during the preliminary work. Introduced in model solutions enable the use of multi-criteria analysis of the choice of the most optimal solutions in terms of costs or expenditures of time during the process of construction. Taking the above work by building a numerical model of the object allows every step of verification by authorized person inventoried solutions and enables tracking and changes in the situation of those found derogations in relation to the parameters established at the primary stage. In the event of significant deviations, there is the possibility of rapid changes to the completed process of calculation and presentation of alternative solutions. Availability software using BIM technology is increasingly common here knowledge of the implementation of such solutions will become in a short time, the standard for most objects or engineering structures. The use of modern solutions using the described processes will be discussed on the example of an industrial facility where there was a need for installation of new equipment and adapt it to the technical parameters.

On the Solution of the Three-Dimensional Flowfield About a Flow-Through Nacelle. Ph.D. Thesis

NASA Technical Reports Server (NTRS)

Compton, William Bernard

1985-01-01

The solution of the three dimensional flow field for a flow through nacelle was studied. Both inviscid and viscous inviscid interacting solutions were examined. Inviscid solutions were obtained with two different computational procedures for solving the three dimensional Euler equations. The first procedure employs an alternating direction implicit numerical algorithm, and required the development of a complete computational model for the nacelle problem. The second computational technique employs a fourth order Runge-Kutta numerical algorithm which was modified to fit the nacelle problem. Viscous effects on the flow field were evaluated with a viscous inviscid interacting computational model. This model was constructed by coupling the explicit Euler solution procedure with a flag entrainment boundary layer solution procedure in a global iteration scheme. The computational techniques were used to compute the flow field for a long duct turbofan engine nacelle at free stream Mach numbers of 0.80 and 0.94 and angles of attack of 0 and 4 deg.

An enriched finite element method to fractional advection-diffusion equation

NASA Astrophysics Data System (ADS)

Luan, Shengzhi; Lian, Yanping; Ying, Yuping; Tang, Shaoqiang; Wagner, Gregory J.; Liu, Wing Kam

2017-08-01

In this paper, an enriched finite element method with fractional basis [ 1,x^{α }] for spatial fractional partial differential equations is proposed to obtain more stable and accurate numerical solutions. For pure fractional diffusion equation without advection, the enriched Galerkin finite element method formulation is demonstrated to simulate the exact solution successfully without any numerical oscillation, which is advantageous compared to the traditional Galerkin finite element method with integer basis [ 1,x] . For fractional advection-diffusion equation, the oscillatory behavior becomes complex due to the introduction of the advection term which can be characterized by a fractional element Peclet number. For the purpose of addressing the more complex numerical oscillation, an enriched Petrov-Galerkin finite element method is developed by using a dimensionless fractional stabilization parameter, which is formulated through a minimization of the residual of the nodal solution. The effectiveness and accuracy of the enriched finite element method are demonstrated by a series of numerical examples of fractional diffusion equation and fractional advection-diffusion equation, including both one-dimensional and two-dimensional, steady-state and time-dependent cases.

Numerical analysis of a main crack interactions with micro-defects/inhomogeneities using two-scale generalized/extended finite element method

NASA Astrophysics Data System (ADS)

Malekan, Mohammad; Barros, Felício B.

2017-12-01

Generalized or extended finite element method (G/XFEM) models the crack by enriching functions of partition of unity type with discontinuous functions that represent well the physical behavior of the problem. However, this enrichment functions are not available for all problem types. Thus, one can use numerically-built (global-local) enrichment functions to have a better approximate procedure. This paper investigates the effects of micro-defects/inhomogeneities on a main crack behavior by modeling the micro-defects/inhomogeneities in the local problem using a two-scale G/XFEM. The global-local enrichment functions are influenced by the micro-defects/inhomogeneities from the local problem and thus change the approximate solution of the global problem with the main crack. This approach is presented in detail by solving three different linear elastic fracture mechanics problems for different cases: two plane stress and a Reissner-Mindlin plate problems. The numerical results obtained with the two-scale G/XFEM are compared with the reference solutions from the analytical, numerical solution using standard G/XFEM method and ABAQUS as well, and from the literature.

Numerical solutions of the semiclassical Boltzmann ellipsoidal-statistical kinetic model equation

PubMed Central

Yang, Jaw-Yen; Yan, Chin-Yuan; Huang, Juan-Chen; Li, Zhihui

2014-01-01

Computations of rarefied gas dynamical flows governed by the semiclassical Boltzmann ellipsoidal-statistical (ES) kinetic model equation using an accurate numerical method are presented. The semiclassical ES model was derived through the maximum entropy principle and conserves not only the mass, momentum and energy, but also contains additional higher order moments that differ from the standard quantum distributions. A different decoding procedure to obtain the necessary parameters for determining the ES distribution is also devised. The numerical method in phase space combines the discrete-ordinate method in momentum space and the high-resolution shock capturing method in physical space. Numerical solutions of two-dimensional Riemann problems for two configurations covering various degrees of rarefaction are presented and various contours of the quantities unique to this new model are illustrated. When the relaxation time becomes very small, the main flow features a display similar to that of ideal quantum gas dynamics, and the present solutions are found to be consistent with existing calculations for classical gas. The effect of a parameter that permits an adjustable Prandtl number in the flow is also studied. PMID:25104904

Propagation of Finite Amplitude Sound in Multiple Waveguide Modes.

NASA Astrophysics Data System (ADS)

van Doren, Thomas Walter

1993-01-01

This dissertation describes a theoretical and experimental investigation of the propagation of finite amplitude sound in multiple waveguide modes. Quasilinear analytical solutions of the full second order nonlinear wave equation, the Westervelt equation, and the KZK parabolic wave equation are obtained for the fundamental and second harmonic sound fields in a rectangular rigid-wall waveguide. It is shown that the Westervelt equation is an acceptable approximation of the full nonlinear wave equation for describing guided sound waves of finite amplitude. A system of first order equations based on both a modal and harmonic expansion of the Westervelt equation is developed for waveguides with locally reactive wall impedances. Fully nonlinear numerical solutions of the system of coupled equations are presented for waveguides formed by two parallel planes which are either both rigid, or one rigid and one pressure release. These numerical solutions are compared to finite -difference solutions of the KZK equation, and it is shown that solutions of the KZK equation are valid only at frequencies which are high compared to the cutoff frequencies of the most important modes of propagation (i.e., for which sound propagates at small grazing angles). Numerical solutions of both the Westervelt and KZK equations are compared to experiments performed in an air-filled, rigid-wall, rectangular waveguide. Solutions of the Westervelt equation are in good agreement with experiment for low source frequencies, at which sound propagates at large grazing angles, whereas solutions of the KZK equation are not valid for these cases. At higher frequencies, at which sound propagates at small grazing angles, agreement between numerical solutions of the Westervelt and KZK equations and experiment is only fair, because of problems in specifying the experimental source condition with sufficient accuracy.

Novel asymmetric representation method for solving the higher-order Ginzburg-Landau equation

PubMed Central

Wong, Pring; Pang, Lihui; Wu, Ye; Lei, Ming; Liu, Wenjun

2016-01-01

In ultrafast optics, optical pulses are generated to be of shorter pulse duration, which has enormous significance to industrial applications and scientific research. The ultrashort pulse evolution in fiber lasers can be described by the higher-order Ginzburg-Landau (GL) equation. However, analytic soliton solutions for this equation have not been obtained by use of existing methods. In this paper, a novel method is proposed to deal with this equation. The analytic soliton solution is obtained for the first time, and is proved to be stable against amplitude perturbations. Through the split-step Fourier method, the bright soliton solution is studied numerically. The analytic results here may extend the integrable methods, and could be used to study soliton dynamics for some equations in other disciplines. It may also provide the other way to obtain two-soliton solutions for higher-order GL equations. PMID:27086841

Aeroacoustic Simulations of a Nose Landing Gear Using FUN3D on Pointwise Unstructured Grids

NASA Technical Reports Server (NTRS)

Vatsa, Veer N.; Khorrami, Mehdi R.; Rhoads, John; Lockard, David P.

2015-01-01

Numerical simulations have been performed for a partially-dressed, cavity-closed (PDCC) nose landing gear configuration that was tested in the University of Florida's open-jet acoustic facility known as the UFAFF. The unstructured-grid flow solver FUN3D is used to compute the unsteady flow field for this configuration. Mixed-element grids generated using the Pointwise(TradeMark) grid generation software are used for these simulations. Particular care is taken to ensure quality cells and proper resolution in critical areas of interest in an effort to minimize errors introduced by numerical artifacts. A hybrid Reynolds-averaged Navier-Stokes/large eddy simulation (RANS/LES) turbulence model is used for these simulations. Solutions are also presented for a wall function model coupled to the standard turbulence model. Time-averaged and instantaneous solutions obtained on these Pointwise grids are compared with the measured data and previous numerical solutions. The resulting CFD solutions are used as input to a Ffowcs Williams-Hawkings noise propagation code to compute the farfield noise levels in the flyover and sideline directions. The computed noise levels compare well with previous CFD solutions and experimental data.

Analytical solutions for coagulation and condensation kinetics of composite particles

NASA Astrophysics Data System (ADS)

Piskunov, Vladimir N.

2013-04-01

The processes of composite particles formation consisting of a mixture of different materials are essential for many practical problems: for analysis of the consequences of accidental releases in atmosphere; for simulation of precipitation formation in clouds; for description of multi-phase processes in chemical reactors and industrial facilities. Computer codes developed for numerical simulation of these processes require optimization of computational methods and verification of numerical programs. Kinetic equations of composite particle formation are given in this work in a concise form (impurity integrated). Coagulation, condensation and external sources associated with nucleation are taken into account. Analytical solutions were obtained in a number of model cases. The general laws for fraction redistribution of impurities were defined. The results can be applied to develop numerical algorithms considerably reducing the simulation effort, as well as to verify the numerical programs for calculation of the formation kinetics of composite particles in the problems of practical importance.

Numerical Solution of the Kzk Equation for Pulsed Finite Amplitude Sound Beams in Thermoviscous Fluids

NASA Astrophysics Data System (ADS)

Lee, Yang-Sub

A time-domain numerical algorithm for solving the KZK (Khokhlov-Zabolotskaya-Kuznetsov) nonlinear parabolic wave equation is developed for pulsed, axisymmetric, finite amplitude sound beams in thermoviscous fluids. The KZK equation accounts for the combined effects of diffraction, absorption, and nonlinearity at the same order of approximation. The accuracy of the algorithm is established via comparison with analytical solutions for several limiting cases, and with numerical results obtained from a widely used algorithm for solving the KZK equation in the frequency domain. The time domain algorithm is used to investigate waveform distortion and shock formation in directive sound beams radiated by pulsed circular piston sources. New results include predictions for the entire process of self-demodulation, and for the effect of frequency modulation on pulse envelope distortion. Numerical results are compared with measurements, and focused sources are investigated briefly.

Modeling of Electromagnetic Scattering by Discrete and Discretely Heterogeneous Random Media by Using Numerically Exact Solutions of the Maxwell Equations

NASA Technical Reports Server (NTRS)

Dlugach, Janna M.; Mishchenko, Michael I.

2017-01-01

In this paper, we discuss some aspects of numerical modeling of electromagnetic scattering by discrete random medium by using numerically exact solutions of the macroscopic Maxwell equations. Typical examples of such media are clouds of interstellar dust, clouds of interplanetary dust in the Solar system, dusty atmospheres of comets, particulate planetary rings, clouds in planetary atmospheres, aerosol particles with numerous inclusions and so on. Our study is based on the results of extensive computations of different characteristics of electromagnetic scattering obtained by using the superposition T-matrix method which represents a direct computer solver of the macroscopic Maxwell equations for an arbitrary multisphere configuration. As a result, in particular, we clarify the range of applicability of the low-density theories of radiative transfer and coherent backscattering as well as of widely used effective-medium approximations.

Numerical simulation of two-phase filtration in the near well bore zone

NASA Astrophysics Data System (ADS)

Maksat, Kalimoldayev; Kalipa, Kuspanova; Kulyash, Baisalbayeva; Orken, Mamyrbayev; Assel, Abdildayeva

2018-04-01

On the basis of the fundamental laws of energy conservation, nonstationary processes of filtration of two-phase liquids in multilayered reservoirs in the near well bore zone are considered. Number of reservoirs, fluid pressure in the given reservoirs, reservoir permeability, oil viscosity, etc. are taken into account upon that. Plane-parallel flow and axisymmetric cases have been studied. In the numerical solution, non-structured meshes are used. Closer to the well, the meshes thicken. The integration step over time is defined by the generalized Courant inequality. As a result, there are no large oscillations in the numerical solutions obtained. Oil production rates, Poisson's ratios, D-diameters of the well, filter height, filter permeability, and cumulative thickness of the filter cake and the area have been taken as the main inputs in numerical simulation of non-stationary processes of two-phase filtration.

A new numerical treatment based on Lucas polynomials for 1D and 2D sinh-Gordon equation

NASA Astrophysics Data System (ADS)

Oruç, Ömer

2018-04-01

In this paper, a new mixed method based on Lucas and Fibonacci polynomials is developed for numerical solutions of 1D and 2D sinh-Gordon equations. Firstly time variable discretized by central finite difference and then unknown function and its derivatives are expanded to Lucas series. With the help of these series expansion and Fibonacci polynomials, matrices for differentiation are derived. With this approach, finding the solution of sinh-Gordon equation transformed to finding the solution of an algebraic system of equations. Lucas series coefficients are acquired by solving this system of algebraic equations. Then by plugginging these coefficients into Lucas series expansion numerical solutions can be obtained consecutively. The main objective of this paper is to demonstrate that Lucas polynomial based method is convenient for 1D and 2D nonlinear problems. By calculating L2 and L∞ error norms of some 1D and 2D test problems efficiency and performance of the proposed method is monitored. Acquired accurate results confirm the applicability of the method.

A pertinent approach to solve nonlinear fuzzy integro-differential equations.

PubMed

Narayanamoorthy, S; Sathiyapriya, S P

2016-01-01

Fuzzy integro-differential equations is one of the important parts of fuzzy analysis theory that holds theoretical as well as applicable values in analytical dynamics and so an appropriate computational algorithm to solve them is in essence. In this article, we use parametric forms of fuzzy numbers and suggest an applicable approach for solving nonlinear fuzzy integro-differential equations using homotopy perturbation method. A clear and detailed description of the proposed method is provided. Our main objective is to illustrate that the construction of appropriate convex homotopy in a proper way leads to highly accurate solutions with less computational work. The efficiency of the approximation technique is expressed via stability and convergence analysis so as to guarantee the efficiency and performance of the methodology. Numerical examples are demonstrated to verify the convergence and it reveals the validity of the presented numerical technique. Numerical results are tabulated and examined by comparing the obtained approximate solutions with the known exact solutions. Graphical representations of the exact and acquired approximate fuzzy solutions clarify the accuracy of the approach.

Cubic spline numerical solution of an ablation problem with convective backface cooling

NASA Astrophysics Data System (ADS)

Lin, S.; Wang, P.; Kahawita, R.

1984-08-01

An implicit numerical technique using cubic splines is presented for solving an ablation problem on a thin wall with convective cooling. A non-uniform computational mesh with 6 grid points has been used for the numerical integration. The method has been found to be computationally efficient, providing for the care under consideration of an overall error of about 1 percent. The results obtained indicate that the convective cooling is an important factor in reducing the ablation thickness.

A meshless method for solving two-dimensional variable-order time fractional advection-diffusion equation

NASA Astrophysics Data System (ADS)

Tayebi, A.; Shekari, Y.; Heydari, M. H.

2017-07-01

Several physical phenomena such as transformation of pollutants, energy, particles and many others can be described by the well-known convection-diffusion equation which is a combination of the diffusion and advection equations. In this paper, this equation is generalized with the concept of variable-order fractional derivatives. The generalized equation is called variable-order time fractional advection-diffusion equation (V-OTFA-DE). An accurate and robust meshless method based on the moving least squares (MLS) approximation and the finite difference scheme is proposed for its numerical solution on two-dimensional (2-D) arbitrary domains. In the time domain, the finite difference technique with a θ-weighted scheme and in the space domain, the MLS approximation are employed to obtain appropriate semi-discrete solutions. Since the newly developed method is a meshless approach, it does not require any background mesh structure to obtain semi-discrete solutions of the problem under consideration, and the numerical solutions are constructed entirely based on a set of scattered nodes. The proposed method is validated in solving three different examples including two benchmark problems and an applied problem of pollutant distribution in the atmosphere. In all such cases, the obtained results show that the proposed method is very accurate and robust. Moreover, a remarkable property so-called positive scheme for the proposed method is observed in solving concentration transport phenomena.

A Zonal Approach for the Solution of Coupled Euler and Potential Solutions of Flows with Complex Geometries.

DTIC Science & Technology

1987-06-01

obtained from: A simple numerical intergration scheme is employed to perform the integral in Equations (B2) and (86) along the dividing streamline. A 11 4...angle of attack was small, the dividing streamline remained almost horizontal in this case. Results of a higher angle of attack case, in which the mesh

Weakly nonlinear behavior of a plate thickness-mode piezoelectric transformer.

PubMed

Yang, Jiashi; Chen, Ziguang; Hu, Yuantai; Jiang, Shunong; Guo, Shaohua

2007-04-01

We analyzed the weakly nonlinear behavior of a plate thickness-shear mode piezoelectric transformer near resonance. An approximate analytical solution was obtained. Numerical results based on the analytical solution are presented. It is shown that on one side of the resonant frequency the input-output relation becomes nonlinear, and on the other side the output voltage experiences jumps.

Group invariant solution for a pre-existing fracture driven by a power-law fluid in permeable rock

NASA Astrophysics Data System (ADS)

Fareo, A. G.; Mason, D. P.

2016-06-01

Group invariant analytical and numerical solutions for the evolution of a two-dimensional fracture with nonzero initial length in permeable rock and driven by an incompressible non-Newtonian fluid of power-law rheology are obtained. The effect of fluid leak-off on the evolution of the power-law fluid fracture is investigated.

Modified equations of finite-size layered plates made of orthotropic material. Comparison of the results of numerical calculations with analytical solutions

NASA Astrophysics Data System (ADS)

Volchkov, Yu. M.

2017-09-01

This paper describes the modified bending equations of layered orthotropic plates in the first approximation. The approximation of the solution of the equation of the three-dimensional theory of elasticity by the Legendre polynomial segments is used to obtain differential equations of the elastic layer. For the approximation of equilibrium equations and boundary conditions of three-dimensional theory of elasticity, several approximations of each desired function (stresses and displacements) are used. The stresses at the internal points of the plate are determined from the defining equations for the orthotropic material, averaged with respect to the plate thickness. The construction of the bending equations of layered plates for each layer is carried out with the help of the elastic layer equations and the conjugation conditions on the boundaries between layers, which are conditions for the continuity of normal stresses and displacements. The numerical solution of the problem of bending of the rectangular layered plate obtained with the help of modified equations is compared with an analytical solution. It is determined that the maximum error in determining the stresses does not exceed 3 %.

Singular perturbation analysis of the steady-state Poisson–Nernst–Planck system: Applications to ion channels

PubMed Central

SINGER, A.; GILLESPIE, D.; NORBURY, J.; EISENBERG, R. S.

2009-01-01

Ion channels are proteins with a narrow hole down their middle that control a wide range of biological function by controlling the flow of spherical ions from one macroscopic region to another. Ion channels do not change their conformation on the biological time scale once they are open, so they can be described by a combination of Poisson and drift-diffusion (Nernst–Planck) equations called PNP in biophysics. We use singular perturbation techniques to analyse the steady-state PNP system for a channel with a general geometry and a piecewise constant permanent charge profile. We construct an outer solution for the case of a constant permanent charge density in three dimensions that is also a valid solution of the one-dimensional system. The asymptotical current–voltage (I–V ) characteristic curve of the device (obtained by the singular perturbation analysis) is shown to be a very good approximation of the numerical I–V curve (obtained by solving the system numerically). The physical constraint of non-negative concentrations implies a unique solution, i.e., for each given applied potential there corresponds a unique electric current (relaxing this constraint yields non-physical multiple solutions for sufficiently large voltages). PMID:19809600

A stage structure pest management model with impulsive state feedback control

NASA Astrophysics Data System (ADS)

Pang, Guoping; Chen, Lansun; Xu, Weijian; Fu, Gang

2015-06-01

A stage structure pest management model with impulsive state feedback control is investigated. We get the sufficient condition for the existence of the order-1 periodic solution by differential equation geometry theory and successor function. Further, we obtain a new judgement method for the stability of the order-1 periodic solution of the semi-continuous systems by referencing the stability analysis for limit cycles of continuous systems, which is different from the previous method of analog of Poincarè criterion. Finally, we analyze numerically the theoretical results obtained.

Using exact solutions to develop an implicit scheme for the baroclinic primitive equations

NASA Technical Reports Server (NTRS)

Marchesin, D.

1984-01-01

The exact solutions presently obtained by means of a novel method for nonlinear initial value problems are used in the development of numerical schemes for the computer solution of these problems. The method is applied to a new, fully implicit scheme on a vertical slice of the isentropic baroclinic equations. It was not possible to find a global scale phenomenon that could be simulated by the baroclinic primitive equations on a vertical slice.

Analysis of periodically excited non-linear systems by a parametric continuation technique

NASA Astrophysics Data System (ADS)

Padmanabhan, C.; Singh, R.

1995-07-01

The dynamic behavior and frequency response of harmonically excited piecewise linear and/or non-linear systems has been the subject of several recent investigations. Most of the prior studies employed harmonic balance or Galerkin schemes, piecewise linear techniques, analog simulation and/or direct numerical integration (digital simulation). Such techniques are somewhat limited in their ability to predict all of the dynamic characteristics, including bifurcations leading to the occurrence of unstable, subharmonic, quasi-periodic and/or chaotic solutions. To overcome this problem, a parametric continuation scheme, based on the shooting method, is applied specifically to a periodically excited piecewise linear/non-linear system, in order to improve understanding as well as to obtain the complete dynamic response. Parameter regions exhibiting bifurcations to harmonic, subharmonic or quasi-periodic solutions are obtained quite efficiently and systematically. Unlike other techniques, the proposed scheme can follow period-doubling bifurcations, and with some modifications obtain stable quasi-periodic solutions and their bifurcations. This knowledge is essential in establishing conditions for the occurrence of chaotic oscillations in any non-linear system. The method is first validated through the Duffing oscillator example, the solutions to which are also obtained by conventional one-term harmonic balance and perturbation methods. The second example deals with a clearance non-linearity problem for both harmonic and periodic excitations. Predictions from the proposed scheme match well with available analog simulation data as well as with multi-term harmonic balance results. Potential savings in computational time over direct numerical integration is demonstrated for some of the example cases. Also, this work has filled in some of the solution regimes for an impact pair, which were missed previously in the literature. Finally, one main limitation associated with the proposed procedure is discussed.

Radiative flow of Carreau liquid in presence of Newtonian heating and chemical reaction

NASA Astrophysics Data System (ADS)

Hayat, T.; Ullah, Ikram; Ahmad, B.; Alsaedi, A.

Objective of this article is to investigate the magnetohydrodynamic (MHD) boundary layer stretched flow of Carreau fluid in the presence of Newtonian heating. Sheet is presumed permeable. Analysis is studied in the presence of chemical reaction and thermal radiation. Mathematical formulation is established by using the boundary layer approximations. The resultant nonlinear flow analysis is computed for the convergent solutions. Interval of convergence via numerical data and plots are obtained and verified. Impact of numerous pertinent variables on the velocity, temperature and concentration is outlined. Numerical data for surface drag coefficient, surface heat transfer (local Nusselt number) and mass transfer (local Sherwood number) is executed and inspected. Comparison of skin friction coefficient in limiting case is made for the verification of current derived solutions.

Exact solution for four-order acousto-optic Bragg diffraction with arbitrary initial conditions.

PubMed

Pieper, Ron; Koslover, Deborah; Poon, Ting-Chung

2009-03-01

An exact solution to the four-order acousto-optic (AO) Bragg diffraction problem with arbitrary initial conditions compatible with exact Bragg angle incident light is developed. The solution, obtained by solving a 4th-order differential equation, is formalized into a transition matrix operator predicting diffracted light orders at the exit of the AO cell in terms of the same diffracted light orders at the entrance. It is shown that the transition matrix is unitary and that this unitary matrix condition is sufficient to guarantee energy conservation. A comparison of analytical solutions with numerical predictions validates the formalism. Although not directly related to the approach used to obtain the solution, it was discovered that all four generated eigenvalues from the four-order AO differential matrix operator are expressed simply in terms of Euclid's Divine Proportion.

Evaluation of truncation error and adaptive grid generation for the transonic full potential flow calculations

NASA Technical Reports Server (NTRS)

Nakamura, S.

1983-01-01

The effects of truncation error on the numerical solution of transonic flows using the full potential equation are studied. The effects of adapting grid point distributions to various solution aspects including shock waves is also discussed. A conclusion is that a rapid change of grid spacing is damaging to the accuracy of the flow solution. Therefore, in a solution adaptive grid application an optimal grid is obtained as a tradeoff between the amount of grid refinement and the rate of grid stretching.

Error analysis of analytic solutions for self-excited near-symmetric rigid bodies - A numerical study

NASA Technical Reports Server (NTRS)

Kia, T.; Longuski, J. M.

1984-01-01

Analytic error bounds are presented for the solutions of approximate models for self-excited near-symmetric rigid bodies. The error bounds are developed for analytic solutions to Euler's equations of motion. The results are applied to obtain a simplified analytic solution for Eulerian rates and angles. The results of a sample application of the range and error bound expressions for the case of the Galileo spacecraft experiencing transverse torques demonstrate the use of the bounds in analyses of rigid body spin change maneuvers.

Finite Element Analysis of Magnetic Damping Effects on G-Jitter Induced Fluid Flow

NASA Technical Reports Server (NTRS)

Pan, Bo; Li, Ben Q.; deGroh, Henry C., III

1997-01-01

This paper reports some interim results on numerical modeling and analyses of magnetic damping of g-jitter driven fluid flow in microgravity. A finite element model is developed to represent the fluid flow, thermal and solute transport phenomena in a 2-D cavity under g-jitter conditions with and without an applied magnetic field. The numerical model is checked by comparing with analytical solutions obtained for a simple parallel plate channel flow driven by g-jitter in a transverse magnetic field. The model is then applied to study the effect of steady state g-jitter induced oscillation and on the solute redistribution in the liquid that bears direct relevance to the Bridgman-Stockbarger single crystal growth processes. A selection of computed results is presented and the results indicate that an applied magnetic field can effectively damp the velocity caused by g-jitter and help to reduce the time variation of solute redistribution.

Natural convection heat transfer in an oscillating vertical cylinder

PubMed Central

Ali Shah, Nehad; Tassaddiq, Asifa; Mustapha, Norzieha; Kechil, Seripah Awang

2018-01-01

This paper studies the heat transfer analysis caused due to free convection in a vertically oscillating cylinder. Exact solutions are determined by applying the Laplace and finite Hankel transforms. Expressions for temperature distribution and velocity field corresponding to cosine and sine oscillations are obtained. The solutions that have been obtained for velocity are presented in the forms of transient and post-transient solutions. Moreover, these solutions satisfy both the governing differential equation and all imposed initial and boundary conditions. Numerical computations and graphical illustrations are used in order to study the effects of Prandtl and Grashof numbers on velocity and temperature for various times. The transient solutions for both cosine and sine oscillations are also computed in tables. It is found that, the transient solutions are of considerable interest up to the times t = 15 for cosine oscillations and t = 1.75 for sine oscillations. After these moments, the transient solutions can be neglected and, the fluid moves according with the post-transient solutions. PMID:29304161

Natural convection heat transfer in an oscillating vertical cylinder.

PubMed

Khan, Ilyas; Ali Shah, Nehad; Tassaddiq, Asifa; Mustapha, Norzieha; Kechil, Seripah Awang

2018-01-01

This paper studies the heat transfer analysis caused due to free convection in a vertically oscillating cylinder. Exact solutions are determined by applying the Laplace and finite Hankel transforms. Expressions for temperature distribution and velocity field corresponding to cosine and sine oscillations are obtained. The solutions that have been obtained for velocity are presented in the forms of transient and post-transient solutions. Moreover, these solutions satisfy both the governing differential equation and all imposed initial and boundary conditions. Numerical computations and graphical illustrations are used in order to study the effects of Prandtl and Grashof numbers on velocity and temperature for various times. The transient solutions for both cosine and sine oscillations are also computed in tables. It is found that, the transient solutions are of considerable interest up to the times t = 15 for cosine oscillations and t = 1.75 for sine oscillations. After these moments, the transient solutions can be neglected and, the fluid moves according with the post-transient solutions.

In Search of Grid Converged Solutions

NASA Technical Reports Server (NTRS)

Lockard, David P.

2010-01-01

Assessing solution error continues to be a formidable task when numerically solving practical flow problems. Currently, grid refinement is the primary method used for error assessment. The minimum grid spacing requirements to achieve design order accuracy for a structured-grid scheme are determined for several simple examples using truncation error evaluations on a sequence of meshes. For certain methods and classes of problems, obtaining design order may not be sufficient to guarantee low error. Furthermore, some schemes can require much finer meshes to obtain design order than would be needed to reduce the error to acceptable levels. Results are then presented from realistic problems that further demonstrate the challenges associated with using grid refinement studies to assess solution accuracy.

Solution of second order quasi-linear boundary value problems by a wavelet method

DOE Office of Scientific and Technical Information (OSTI.GOV)

Zhang, Lei; Zhou, Youhe; Wang, Jizeng, E-mail: jzwang@lzu.edu.cn

2015-03-10

A wavelet Galerkin method based on expansions of Coiflet-like scaling function bases is applied to solve second order quasi-linear boundary value problems which represent a class of typical nonlinear differential equations. Two types of typical engineering problems are selected as test examples: one is about nonlinear heat conduction and the other is on bending of elastic beams. Numerical results are obtained by the proposed wavelet method. Through comparing to relevant analytical solutions as well as solutions obtained by other methods, we find that the method shows better efficiency and accuracy than several others, and the rate of convergence can evenmore » reach orders of 5.8.« less

Numerical study of the Kadomtsev-Petviashvili equation and dispersive shock waves

NASA Astrophysics Data System (ADS)

Grava, T.; Klein, C.; Pitton, G.

2018-02-01

A detailed numerical study of the long time behaviour of dispersive shock waves in solutions to the Kadomtsev-Petviashvili (KP) I equation is presented. It is shown that modulated lump solutions emerge from the dispersive shock waves. For the description of dispersive shock waves, Whitham modulation equations for KP are obtained. It is shown that the modulation equations near the soliton line are hyperbolic for the KPII equation while they are elliptic for the KPI equation leading to a focusing effect and the formation of lumps. Such a behaviour is similar to the appearance of breathers for the focusing nonlinear Schrödinger equation in the semiclassical limit.

Incremental analysis of large elastic deformation of a rotating cylinder

NASA Technical Reports Server (NTRS)

Buchanan, G. R.

1976-01-01

The effect of finite deformation upon a rotating, orthotropic cylinder was investigated using a general incremental theory. The incremental equations of motion are developed using the variational principle. The governing equations are derived using the principle of virtual work for a body with initial stress. The governing equations are reduced to those for the title problem and a numerical solution is obtained using finite difference approximations. Since the problem is defined in terms of one independent space coordinate, the finite difference grid can be modified as the incremental deformation occurs without serious numerical difficulties. The nonlinear problem is solved incrementally by totaling a series of linear solutions.

Flow of nanofluid past a Riga plate

NASA Astrophysics Data System (ADS)

Ahmad, Adeel; Asghar, Saleem; Afzal, Sumaira

2016-03-01

This paper studies the mixed convection boundary layer flow of a nanofluid past a vertical Riga plate in the presence of strong suction. The mathematical model incorporates the Brownian motion and thermophoresis effects due to nanofluid and the Grinberg-term for the wall parallel Lorentz force due to Riga plate. The analytical solution of the problem is presented using the perturbation method for small Brownian and thermophoresis diffusion parameters. The numerical solution is also presented to ensure the reliability of the asymptotic method. The comparison of the two solutions shows an excellent agreement. The correlation expressions for skin friction, Nusselt number and Sherwood number are developed by performing linear regression on the obtained numerical data. The effects of nanofluid and the Lorentz force due to Riga plate, on the skin friction are discussed.

Numerical solution of quadratic matrix equations for free vibration analysis of structures

NASA Technical Reports Server (NTRS)

Gupta, K. K.

1975-01-01

This paper is concerned with the efficient and accurate solution of the eigenvalue problem represented by quadratic matrix equations. Such matrix forms are obtained in connection with the free vibration analysis of structures, discretized by finite 'dynamic' elements, resulting in frequency-dependent stiffness and inertia matrices. The paper presents a new numerical solution procedure of the quadratic matrix equations, based on a combined Sturm sequence and inverse iteration technique enabling economical and accurate determination of a few required eigenvalues and associated vectors. An alternative procedure based on a simultaneous iteration procedure is also described when only the first few modes are the usual requirement. The employment of finite dynamic elements in conjunction with the presently developed eigenvalue routines results in a most significant economy in the dynamic analysis of structures.

Boundary Layer Model for Air Pollutant Concentrations Due to Highway Traffic

ERIC Educational Resources Information Center

Ragland, Kenneth W.; Peirce, J. Jeffrey

1975-01-01

A numerical solution of the three-dimensional steady-state diffusion equation for a finite width line source is presented. The wind speed and eddy diffusivity as a function of height above the roadway are obtained. Normalized ground level and elevated concentrations near a highway are obtained for winds perpendicular, parallel, and at 45 degrees.…

Vibrations and stresses in layered anisotropic cylinders

NASA Technical Reports Server (NTRS)

Mulholland, G. P.; Gupta, B. P.

1976-01-01

An equation describing the radial displacement in a k layered anisotropic cylinder was obtained. The cylinders are initially unstressed but are subjected to either a time dependent normal stress or a displacement at the external boundaries of the laminate. The solution is obtained by utilizing the Vodicka orthogonalization technique. Numerical examples are given to illustrate the procedure.

Lax-Friedrichs sweeping scheme for static Hamilton-Jacobi equations

NASA Astrophysics Data System (ADS)

Kao, Chiu Yen; Osher, Stanley; Qian, Jianliang

2004-05-01

We propose a simple, fast sweeping method based on the Lax-Friedrichs monotone numerical Hamiltonian to approximate viscosity solutions of arbitrary static Hamilton-Jacobi equations in any number of spatial dimensions. By using the Lax-Friedrichs numerical Hamiltonian, we can easily obtain the solution at a specific grid point in terms of its neighbors, so that a Gauss-Seidel type nonlinear iterative method can be utilized. Furthermore, by incorporating a group-wise causality principle into the Gauss-Seidel iteration by following a finite group of characteristics, we have an easy-to-implement, sweeping-type, and fast convergent numerical method. However, unlike other methods based on the Godunov numerical Hamiltonian, some computational boundary conditions are needed in the implementation. We give a simple recipe which enforces a version of discrete min-max principle. Some convergence analysis is done for the one-dimensional eikonal equation. Extensive 2-D and 3-D numerical examples illustrate the efficiency and accuracy of the new approach. To our knowledge, this is the first fast numerical method based on discretizing the Hamilton-Jacobi equation directly without assuming convexity and/or homogeneity of the Hamiltonian.

Minimum Altitude-Loss Soaring in a Specified Vertical Wind Distribution

NASA Technical Reports Server (NTRS)

Pierson, B. L.; Chen, I.

1979-01-01

Minimum altitude-loss flight of a sailplane through a given vertical wind distribution is discussed. The problem is posed as an optimal control problem, and several numerical solutions are obtained for a sinusoidal wind distribution.

A method of measuring the effective thermal conductivity of thermoplastic foams

NASA Astrophysics Data System (ADS)

Asséko, André Chateau Akué; Cosson, Benoit; Chaki, Salim; Duborper, Clément; Lacrampe, Marie-France; Krawczak, Patricia

2017-10-01

An inverse method for determining the in-plane effective thermal conductivity of porous thermoplastics was implemented by coupling infrared thermography experiments and numerical solution of heat transfer in straight fins having temperature-dependent convective heat transfer coefficient. The obtained effective thermal conductivity values were compared with previous results obtained using a numerical solution based on periodic homogenization techniques (NSHT) in which the microstructure heterogeneity of extruded polymeric polyethylene (PE) foam in which pores are filled with air with different levels of open and closed porosity was taken into account and Transient Plane Source Technique (TPS) in order to verify the accuracy of the proposed method. The new method proposed in the present study is in good agreement with both NSHT and TPS. It is also applicable to structural materials such as composites, e.g. unidirectional fiber-reinforced plastics, where heat transfer is very different according to the fiber direction (parallel or transverse to the fibers).

A new anisotropic mesh adaptation method based upon hierarchical a posteriori error estimates

NASA Astrophysics Data System (ADS)

Huang, Weizhang; Kamenski, Lennard; Lang, Jens

2010-03-01

A new anisotropic mesh adaptation strategy for finite element solution of elliptic differential equations is presented. It generates anisotropic adaptive meshes as quasi-uniform ones in some metric space, with the metric tensor being computed based on hierarchical a posteriori error estimates. A global hierarchical error estimate is employed in this study to obtain reliable directional information of the solution. Instead of solving the global error problem exactly, which is costly in general, we solve it iteratively using the symmetric Gauß-Seidel method. Numerical results show that a few GS iterations are sufficient for obtaining a reasonably good approximation to the error for use in anisotropic mesh adaptation. The new method is compared with several strategies using local error estimators or recovered Hessians. Numerical results are presented for a selection of test examples and a mathematical model for heat conduction in a thermal battery with large orthotropic jumps in the material coefficients.

Elasticity Solution of an Adhesively Bonded Cover Plate of Various Geometries

NASA Technical Reports Server (NTRS)

Aksel, G. N.; Erdogan, F.

1985-01-01

The plane strain of adhesively bonded structures consisting of two different isotropic adherends is considered. By expressing the x-y components of the displacements in terms of Fourier integrals and using the corresponding boundary and continuity conditions, the integral equations for the general problem are obtained and solved numerically by applying Gauss-Chebyshev integration scheme. The shear and the normal stresses in the adhesive are calculated for various geometries and material properties for a stiffened plate under uniaxial tension. Numerical results involving the stress intensity factors and the strain energy release rate are presented. The closed-form expressions for the Fredholm kernels are provided to obtain the solution for an arbitrary geometry and material properties. For the general geometry, the contribution of the normal stress is quite significant, while for symmetric geometries, the shear stress is dominant, the normal stress vanishes if the adherends are of the same material and the same thickness.

Holographic models with anisotropic scaling

NASA Astrophysics Data System (ADS)

Brynjolfsson, E. J.; Danielsson, U. H.; Thorlacius, L.; Zingg, T.

2013-12-01

We consider gravity duals to d+1 dimensional quantum critical points with anisotropic scaling. The primary motivation comes from strongly correlated electron systems in condensed matter theory but the main focus of the present paper is on the gravity models in their own right. Physics at finite temperature and fixed charge density is described in terms of charged black branes. Some exact solutions are known and can be used to obtain a maximally extended spacetime geometry, which has a null curvature singularity inside a single non-degenerate horizon, but generic black brane solutions in the model can only be obtained numerically. Charged matter gives rise to black branes with hair that are dual to the superconducting phase of a holographic superconductor. Our numerical results indicate that holographic superconductors with anisotropic scaling have vanishing zero temperature entropy when the back reaction of the hair on the brane geometry is taken into account.

Direct discontinuous Galerkin method and its variations for second order elliptic equations

DOE PAGES

Huang, Hongying; Chen, Zheng; Li, Jin; ...

2016-08-23

In this study, we study direct discontinuous Galerkin method (Liu and Yan in SIAM J Numer Anal 47(1):475–698, 2009) and its variations (Liu and Yan in Commun Comput Phys 8(3):541–564, 2010; Vidden and Yan in J Comput Math 31(6):638–662, 2013; Yan in J Sci Comput 54(2–3):663–683, 2013) for 2nd order elliptic problems. A priori error estimate under energy norm is established for all four methods. Optimal error estimate under L 2 norm is obtained for DDG method with interface correction (Liu and Yan in Commun Comput Phys 8(3):541–564, 2010) and symmetric DDG method (Vidden and Yan in J Comput Mathmore » 31(6):638–662, 2013). A series of numerical examples are carried out to illustrate the accuracy and capability of the schemes. Numerically we obtain optimal (k+1)th order convergence for DDG method with interface correction and symmetric DDG method on nonuniform and unstructured triangular meshes. An interface problem with discontinuous diffusion coefficients is investigated and optimal (k+1)th order accuracy is obtained. Peak solutions with sharp transitions are captured well. Highly oscillatory wave solutions of Helmholz equation are well resolved.« less

Direct discontinuous Galerkin method and its variations for second order elliptic equations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Huang, Hongying; Chen, Zheng; Li, Jin

In this study, we study direct discontinuous Galerkin method (Liu and Yan in SIAM J Numer Anal 47(1):475–698, 2009) and its variations (Liu and Yan in Commun Comput Phys 8(3):541–564, 2010; Vidden and Yan in J Comput Math 31(6):638–662, 2013; Yan in J Sci Comput 54(2–3):663–683, 2013) for 2nd order elliptic problems. A priori error estimate under energy norm is established for all four methods. Optimal error estimate under L 2 norm is obtained for DDG method with interface correction (Liu and Yan in Commun Comput Phys 8(3):541–564, 2010) and symmetric DDG method (Vidden and Yan in J Comput Mathmore » 31(6):638–662, 2013). A series of numerical examples are carried out to illustrate the accuracy and capability of the schemes. Numerically we obtain optimal (k+1)th order convergence for DDG method with interface correction and symmetric DDG method on nonuniform and unstructured triangular meshes. An interface problem with discontinuous diffusion coefficients is investigated and optimal (k+1)th order accuracy is obtained. Peak solutions with sharp transitions are captured well. Highly oscillatory wave solutions of Helmholz equation are well resolved.« less

Correlation between discrete probability and reaction front propagation rate in heterogeneous mixtures

NASA Astrophysics Data System (ADS)

Naine, Tarun Bharath; Gundawar, Manoj Kumar

2017-09-01

We demonstrate a very powerful correlation between the discrete probability of distances of neighboring cells and thermal wave propagation rate, for a system of cells spread on a one-dimensional chain. A gamma distribution is employed to model the distances of neighboring cells. In the absence of an analytical solution and the differences in ignition times of adjacent reaction cells following non-Markovian statistics, invariably the solution for thermal wave propagation rate for a one-dimensional system with randomly distributed cells is obtained by numerical simulations. However, such simulations which are based on Monte-Carlo methods require several iterations of calculations for different realizations of distribution of adjacent cells. For several one-dimensional systems, differing in the value of shaping parameter of the gamma distribution, we show that the average reaction front propagation rates obtained by a discrete probability between two limits, shows excellent agreement with those obtained numerically. With the upper limit at 1.3, the lower limit depends on the non-dimensional ignition temperature. Additionally, this approach also facilitates the prediction of burning limits of heterogeneous thermal mixtures. The proposed method completely eliminates the need for laborious, time intensive numerical calculations where the thermal wave propagation rates can now be calculated based only on macroscopic entity of discrete probability.

Techniques of orbital decay and long-term ephemeris prediction for satellites in earth orbit

NASA Technical Reports Server (NTRS)

Barry, B. F.; Pimm, R. S.; Rowe, C. K.

1971-01-01

In the special perturbation method, Cowell and variation-of-parameters formulations of the motion equations are implemented and numerically integrated. Variations in the orbital elements due to drag are computed using the 1970 Jacchia atmospheric density model, which includes the effects of semiannual variations, diurnal bulge, solar activity, and geomagnetic activity. In the general perturbation method, two-variable asymptotic series and automated manipulation capabilities are used to obtain analytical solutions to the variation-of-parameters equations. Solutions are obtained considering the effect of oblateness only and the combined effects of oblateness and drag. These solutions are then numerically evaluated by means of a FORTRAN program in which an updating scheme is used to maintain accurate epoch values of the elements. The atmospheric density function is approximated by a Fourier series in true anomaly, and the 1970 Jacchia model is used to periodically update the Fourier coefficients. The accuracy of both methods is demonstrated by comparing computed orbital elements to actual elements over time spans of up to 8 days for the special perturbation method and up to 356 days for the general perturbation method.

Time-dependent flow model of a generalized Burgers' fluid with fractional derivatives through a cylindrical domain: An exact and numerical approach

NASA Astrophysics Data System (ADS)

Safdar, Rabia; Imran, M.; Khalique, Chaudry Masood

2018-06-01

Exact solutions for velocity field and tangential stress for rotational flow of a generalized Burgers' fluid within an infinite circular pipe are derived by using the methods of Laplace and finite Hankel transformations. Firstly we take the position of fluid at rest and then the fluid flow due to the rotation of the pipe around the axis of flow having time dependant angular velocity. The exact solutions are presented in terms of the generalized Ga,b,c (., t) -functions. The corresponding results can be freely specified for the same results of Burgers', Oldroyd B, Maxwell, second grade and Newtonian fluids (performing the same motion) as particular cases of the results obtained earlier. The impact of the different parameters, individually and in comparison, are represented by graphical demonstrations. Secondly the numerical solutions for velocity and stress are also obtained with the help of Laplace transformation, Gaver Stehfest's algorithm and MATHCAD. Finally a comparison of both methods for the same problem is done and shows the consistency of results.

Concentration history during pumping from a leaky aquifer with stratified initial concentration

USGS Publications Warehouse

Goode, Daniel J.; Hsieh, Paul A.; Shapiro, Allen M.; Wood, Warren W.; Kraemer, Thomas F.

1993-01-01

Analytical and numerical solutions are employed to examine the concentration history of a dissolved substance in water pumped from a leaky aquifer. Many aquifer systems are characterized by stratification, for example, a sandy layer overlain by a clay layer. To obtain information about separate hydrogeologic units, aquifer pumping tests are often conducted with a well penetrating only one of the layers. When the initial concentration distribution is also stratified (the concentration varies with elevation only), the concentration breakthrough in the pumped well may be interpreted to provide information on aquifer hydraulic and transport properties. To facilitate this interpretation, we present some simple analytical and numerical solutions for limiting cases and illustrate their application to a fractured bedrock/glacial drift aquifer system where the solute of interest is dissolved radon gas. In addition to qualitative information on water source, this method may yield estimates of effective porosity and saturated thickness (or fracture transport aperture) from a single-hole test. Little information about dispersivity is obtained because the measured concentration is not significantly affected by dispersion in the aquifer.

Space-time domain solutions of the wave equation by a non-singular boundary integral method and Fourier transform.

PubMed

Klaseboer, Evert; Sepehrirahnama, Shahrokh; Chan, Derek Y C

2017-08-01

The general space-time evolution of the scattering of an incident acoustic plane wave pulse by an arbitrary configuration of targets is treated by employing a recently developed non-singular boundary integral method to solve the Helmholtz equation in the frequency domain from which the space-time solution of the wave equation is obtained using the fast Fourier transform. The non-singular boundary integral solution can enforce the radiation boundary condition at infinity exactly and can account for multiple scattering effects at all spacings between scatterers without adverse effects on the numerical precision. More generally, the absence of singular kernels in the non-singular integral equation confers high numerical stability and precision for smaller numbers of degrees of freedom. The use of fast Fourier transform to obtain the time dependence is not constrained to discrete time steps and is particularly efficient for studying the response to different incident pulses by the same configuration of scatterers. The precision that can be attained using a smaller number of Fourier components is also quantified.

Comptonization in Ultra-Strong Magnetic Fields: Numerical Solution to the Radiative Transfer Problem

NASA Technical Reports Server (NTRS)

Ceccobello, C.; Farinelli, R.; Titarchuk, L.

2014-01-01

We consider the radiative transfer problem in a plane-parallel slab of thermal electrons in the presence of an ultra-strong magnetic field (B approximately greater than B(sub c) approx. = 4.4 x 10(exp 13) G). Under these conditions, the magnetic field behaves like a birefringent medium for the propagating photons, and the electromagnetic radiation is split into two polarization modes, ordinary and extraordinary, that have different cross-sections. When the optical depth of the slab is large, the ordinary-mode photons are strongly Comptonized and the photon field is dominated by an isotropic component. Aims. The radiative transfer problem in strong magnetic fields presents many mathematical issues and analytical or numerical solutions can be obtained only under some given approximations. We investigate this problem both from the analytical and numerical point of view, provide a test of the previous analytical estimates, and extend these results with numerical techniques. Methods. We consider here the case of low temperature black-body photons propagating in a sub-relativistic temperature plasma, which allows us to deal with a semi-Fokker-Planck approximation of the radiative transfer equation. The problem can then be treated with the variable separation method, and we use a numerical technique to find solutions to the eigenvalue problem in the case of a singular kernel of the space operator. The singularity of the space kernel is the result of the strong angular dependence of the electron cross-section in the presence of a strong magnetic field. Results. We provide the numerical solution obtained for eigenvalues and eigenfunctions of the space operator, and the emerging Comptonization spectrum of the ordinary-mode photons for any eigenvalue of the space equation and for energies significantly lesser than the cyclotron energy, which is on the order of MeV for the intensity of the magnetic field here considered. Conclusions. We derived the specific intensity of the ordinary photons, under the approximation of large angle and large optical depth. These assumptions allow the equation to be treated using a diffusion-like approximation.

Numerical solution to generalized Burgers'-Fisher equation using Exp-function method hybridized with heuristic computation.

PubMed

Malik, Suheel Abdullah; Qureshi, Ijaz Mansoor; Amir, Muhammad; Malik, Aqdas Naveed; Haq, Ihsanul

2015-01-01

In this paper, a new heuristic scheme for the approximate solution of the generalized Burgers'-Fisher equation is proposed. The scheme is based on the hybridization of Exp-function method with nature inspired algorithm. The given nonlinear partial differential equation (NPDE) through substitution is converted into a nonlinear ordinary differential equation (NODE). The travelling wave solution is approximated by the Exp-function method with unknown parameters. The unknown parameters are estimated by transforming the NODE into an equivalent global error minimization problem by using a fitness function. The popular genetic algorithm (GA) is used to solve the minimization problem, and to achieve the unknown parameters. The proposed scheme is successfully implemented to solve the generalized Burgers'-Fisher equation. The comparison of numerical results with the exact solutions, and the solutions obtained using some traditional methods, including adomian decomposition method (ADM), homotopy perturbation method (HPM), and optimal homotopy asymptotic method (OHAM), show that the suggested scheme is fairly accurate and viable for solving such problems.

Numerical Solution to Generalized Burgers'-Fisher Equation Using Exp-Function Method Hybridized with Heuristic Computation

PubMed Central

Malik, Suheel Abdullah; Qureshi, Ijaz Mansoor; Amir, Muhammad; Malik, Aqdas Naveed; Haq, Ihsanul

2015-01-01

In this paper, a new heuristic scheme for the approximate solution of the generalized Burgers'-Fisher equation is proposed. The scheme is based on the hybridization of Exp-function method with nature inspired algorithm. The given nonlinear partial differential equation (NPDE) through substitution is converted into a nonlinear ordinary differential equation (NODE). The travelling wave solution is approximated by the Exp-function method with unknown parameters. The unknown parameters are estimated by transforming the NODE into an equivalent global error minimization problem by using a fitness function. The popular genetic algorithm (GA) is used to solve the minimization problem, and to achieve the unknown parameters. The proposed scheme is successfully implemented to solve the generalized Burgers'-Fisher equation. The comparison of numerical results with the exact solutions, and the solutions obtained using some traditional methods, including adomian decomposition method (ADM), homotopy perturbation method (HPM), and optimal homotopy asymptotic method (OHAM), show that the suggested scheme is fairly accurate and viable for solving such problems. PMID:25811858

Numerical computation of diffusion on a surface.

PubMed

Schwartz, Peter; Adalsteinsson, David; Colella, Phillip; Arkin, Adam Paul; Onsum, Matthew

2005-08-09

We present a numerical method for computing diffusive transport on a surface derived from image data. Our underlying discretization method uses a Cartesian grid embedded boundary method for computing the volume transport in a region consisting of all points a small distance from the surface. We obtain a representation of this region from image data by using a front propagation computation based on level set methods for solving the Hamilton-Jacobi and eikonal equations. We demonstrate that the method is second-order accurate in space and time and is capable of computing solutions on complex surface geometries obtained from image data of cells.

S4 solution of the transport equation for eigenvalues using Legendre polynomials

NASA Astrophysics Data System (ADS)

Öztürk, Hakan; Bülbül, Ahmet

2017-09-01

Numerical solution of the transport equation for monoenergetic neutrons scattered isotropically through the medium of a finite homogeneous slab is studied for the determination of the eigenvalues. After obtaining the discrete ordinates form of the transport equation, separated homogeneous and particular solutions are formed and then the eigenvalues are calculated using the Gauss-Legendre quadrature set. Then, the calculated eigenvalues for various values of the c0, the mean number of secondary neutrons per collision, are given in the tables.

Development of a time-dependent incompressible Navier-Stokes solver based on a fractional-step method

NASA Technical Reports Server (NTRS)

Rosenfeld, Moshe

1990-01-01

The main goals are the development, validation, and application of a fractional step solution method of the time-dependent incompressible Navier-Stokes equations in generalized coordinate systems. A solution method that combines a finite volume discretization with a novel choice of the dependent variables and a fractional step splitting to obtain accurate solutions in arbitrary geometries is extended to include more general situations, including cases with moving grids. The numerical techniques are enhanced to gain efficiency and generality.

On Solutions for the Transient Response of Beams

NASA Technical Reports Server (NTRS)

Leonard, Robert W.

1959-01-01

Williams type modal solutions of the elementary and Timoshenko beam equations are presented for the response of several uniform beams to a general applied load. Example computations are shown for a free-free beam subject to various concentrated loads at its center. Discussion includes factors influencing the convergence of modal solutions and factors to be considered in a choice of beam theory. Results obtained by two numerical procedures, the traveling-wave method and Houbolt's method, are also presented and discussed.

Impacts of the horizontal and vertical grids on the numerical solutions of the dynamical equations - Part 1: Nonhydrostatic inertia-gravity modes

NASA Astrophysics Data System (ADS)

Konor, Celal S.; Randall, David A.

2018-05-01

We have used a normal-mode analysis to investigate the impacts of the horizontal and vertical discretizations on the numerical solutions of the nonhydrostatic anelastic inertia-gravity modes on a midlatitude f plane. The dispersion equations are derived from the linearized anelastic equations that are discretized on the Z, C, D, CD, (DC), A, E and B horizontal grids, and on the L and CP vertical grids. The effects of both horizontal grid spacing and vertical wavenumber are analyzed, and the role of nonhydrostatic effects is discussed. We also compare the results of the normal-mode analyses with numerical solutions obtained by running linearized numerical models based on the various horizontal grids. The sources and behaviors of the computational modes in the numerical simulations are also examined.Our normal-mode analyses with the Z, C, D, A, E and B grids generally confirm the conclusions of previous shallow-water studies for the cyclone-resolving scales (with low horizontal wavenumbers). We conclude that, aided by nonhydrostatic effects, the Z and C grids become overall more accurate for cloud-resolving resolutions (with high horizontal wavenumbers) than for the cyclone-resolving scales.A companion paper, Part 2, discusses the impacts of the discretization on the Rossby modes on a midlatitude β plane.

Analytical theory of two-dimensional ring dark soliton in nonlocal nonlinear media

NASA Astrophysics Data System (ADS)

Chen, Wei; Wang, Qi; Shi, Jielong; Shen, Ming

2017-11-01

Completely stable two-dimensional ring dark soliton in nonlocal media with an arbitrary degree of nonlocality are investigated. The exact solution of the ring dark solitons is obtained with the variational method and a cylindrical nonlocal response function. The analytical results are confirmed by directly numerical simulations. We also analytically and numerically study the expansion dynamics of the gray ring dark solitons in detail.

Approximate method for calculating a thickwalled cylinder with rigidly clamped ends

NASA Astrophysics Data System (ADS)

Andreev, Vladimir

2018-03-01

Numerous papers dealing with the calculations of cylindrical bodies [1 -8 and others] have shown that analytic and numerical-analytical solutions in both homogeneous and inhomogeneous thick-walled shells can be obtained quite simply, using expansions in Fourier series on trigonometric functions, if the ends are hinged movable (sliding support). It is much more difficult to solve the problem of calculating shells with builtin ends.

A direct Primitive Variable Recovery Scheme for hyperbolic conservative equations: The case of relativistic hydrodynamics.

PubMed

Aguayo-Ortiz, A; Mendoza, S; Olvera, D

2018-01-01

In this article we develop a Primitive Variable Recovery Scheme (PVRS) to solve any system of coupled differential conservative equations. This method obtains directly the primitive variables applying the chain rule to the time term of the conservative equations. With this, a traditional finite volume method for the flux is applied in order avoid violation of both, the entropy and "Rankine-Hugoniot" jump conditions. The time evolution is then computed using a forward finite difference scheme. This numerical technique evades the recovery of the primitive vector by solving an algebraic system of equations as it is often used and so, it generalises standard techniques to solve these kind of coupled systems. The article is presented bearing in mind special relativistic hydrodynamic numerical schemes with an added pedagogical view in the appendix section in order to easily comprehend the PVRS. We present the convergence of the method for standard shock-tube problems of special relativistic hydrodynamics and a graphical visualisation of the errors using the fluctuations of the numerical values with respect to exact analytic solutions. The PVRS circumvents the sometimes arduous computation that arises from standard numerical methods techniques, which obtain the desired primitive vector solution through an algebraic polynomial of the charges.

A direct Primitive Variable Recovery Scheme for hyperbolic conservative equations: The case of relativistic hydrodynamics

PubMed Central

Mendoza, S.; Olvera, D.

2018-01-01

In this article we develop a Primitive Variable Recovery Scheme (PVRS) to solve any system of coupled differential conservative equations. This method obtains directly the primitive variables applying the chain rule to the time term of the conservative equations. With this, a traditional finite volume method for the flux is applied in order avoid violation of both, the entropy and “Rankine-Hugoniot” jump conditions. The time evolution is then computed using a forward finite difference scheme. This numerical technique evades the recovery of the primitive vector by solving an algebraic system of equations as it is often used and so, it generalises standard techniques to solve these kind of coupled systems. The article is presented bearing in mind special relativistic hydrodynamic numerical schemes with an added pedagogical view in the appendix section in order to easily comprehend the PVRS. We present the convergence of the method for standard shock-tube problems of special relativistic hydrodynamics and a graphical visualisation of the errors using the fluctuations of the numerical values with respect to exact analytic solutions. The PVRS circumvents the sometimes arduous computation that arises from standard numerical methods techniques, which obtain the desired primitive vector solution through an algebraic polynomial of the charges. PMID:29659602

Analysis of groundwater flow and stream depletion in L-shaped fluvial aquifers

NASA Astrophysics Data System (ADS)

Lin, Chao-Chih; Chang, Ya-Chi; Yeh, Hund-Der

2018-04-01

Understanding the head distribution in aquifers is crucial for the evaluation of groundwater resources. This article develops a model for describing flow induced by pumping in an L-shaped fluvial aquifer bounded by impermeable bedrocks and two nearly fully penetrating streams. A similar scenario for numerical studies was reported in Kihm et al. (2007). The water level of the streams is assumed to be linearly varying with distance. The aquifer is divided into two subregions and the continuity conditions of the hydraulic head and flux are imposed at the interface of the subregions. The steady-state solution describing the head distribution for the model without pumping is first developed by the method of separation of variables. The transient solution for the head distribution induced by pumping is then derived based on the steady-state solution as initial condition and the methods of finite Fourier transform and Laplace transform. Moreover, the solution for stream depletion rate (SDR) from each of the two streams is also developed based on the head solution and Darcy's law. Both head and SDR solutions in the real time domain are obtained by a numerical inversion scheme called the Stehfest algorithm. The software MODFLOW is chosen to compare with the proposed head solution for the L-shaped aquifer. The steady-state and transient head distributions within the L-shaped aquifer predicted by the present solution are compared with the numerical simulations and measurement data presented in Kihm et al. (2007).

Lattice Boltzmann Equation On a 2D Rectangular Grid

NASA Technical Reports Server (NTRS)

Bouzidi, MHamed; DHumieres, Dominique; Lallemand, Pierre; Luo, Li-Shi; Bushnell, Dennis M. (Technical Monitor)

2002-01-01

We construct a multi-relaxation lattice Boltzmann model on a two-dimensional rectangular grid. The model is partly inspired by a previous work of Koelman to construct a lattice BGK model on a two-dimensional rectangular grid. The linearized dispersion equation is analyzed to obtain the constraints on the isotropy of the transport coefficients and Galilean invariance for various wave propagations in the model. The linear stability of the model is also studied. The model is numerically tested for three cases: (a) a vortex moving with a constant velocity on a mesh periodic boundary conditions; (b) Poiseuille flow with an arbitrasy inclined angle with respect to the lattice orientation: and (c) a cylinder &symmetrically placed in a channel. The numerical results of these tests are compared with either analytic solutions or the results obtained by other methods. Satisfactory results are obtained for the numerical simulations.

Finite length Taylor Couette flow

NASA Technical Reports Server (NTRS)

Streett, C. L.; Hussaini, M. Y.

1987-01-01

Axisymmetric numerical solutions of the unsteady Navier-Stokes equations for flow between concentric rotating cylinders of finite length are obtained by a spectral collocation method. These representative results pertain to two-cell/one-cell exchange process, and are compared with recent experiments.

Prediction of discretization error using the error transport equation

NASA Astrophysics Data System (ADS)

Celik, Ismail B.; Parsons, Don Roscoe

2017-06-01

This study focuses on an approach to quantify the discretization error associated with numerical solutions of partial differential equations by solving an error transport equation (ETE). The goal is to develop a method that can be used to adequately predict the discretization error using the numerical solution on only one grid/mesh. The primary problem associated with solving the ETE is the formulation of the error source term which is required for accurately predicting the transport of the error. In this study, a novel approach is considered which involves fitting the numerical solution with a series of locally smooth curves and then blending them together with a weighted spline approach. The result is a continuously differentiable analytic expression that can be used to determine the error source term. Once the source term has been developed, the ETE can easily be solved using the same solver that is used to obtain the original numerical solution. The new methodology is applied to the two-dimensional Navier-Stokes equations in the laminar flow regime. A simple unsteady flow case is also considered. The discretization error predictions based on the methodology presented in this study are in good agreement with the 'true error'. While in most cases the error predictions are not quite as accurate as those from Richardson extrapolation, the results are reasonable and only require one numerical grid. The current results indicate that there is much promise going forward with the newly developed error source term evaluation technique and the ETE.

Black hole solutions in d = 5 Chern-Simons gravity

NASA Astrophysics Data System (ADS)

Brihaye, Yves; Radu, Eugen

2013-11-01

The five dimensional Einstein-Gauss-Bonnet gravity with a negative cosmological constant becomes, for a special value of the Gauss-Bonnet coupling constant, a Chern-Simons (CS) theory of gravity. In this work we discuss the properties of several different types of black object solutions of this model. Special attention is paid to the case of spinning black holes with equal-magnitude angular momenta which posses a regular horizon of spherical topology. Closed form solutions are obtained in the small angular momentum limit. Nonperturbative solutions are constructed by solving numerically the equations of the model. Apart from that, new exact solutions describing static squashed black holes and black strings are also discussed. The action and global charges of all configurations studied in this work are obtained by using the quasilocal formalism with boundary counterterms generalized for the case of a d = 5 CS theory.

Wave transience in a compressible atmosphere. I - Transient internal wave, mean-flow interaction. II - Transient equatorial waves in the quasi-biennial oscillation

NASA Technical Reports Server (NTRS)

Dunkerton, T. J.

1981-01-01

Analytical and numerical solutions are obtained in an approximate quasi-linear model, to describe the way in which vertically propagating waves give rise to mean flow accelerations in an atmosphere due to the effects of wave transience. These effects in turn result from compressibility and vertical group velocity feedback, and culminate in the spontaneous formation and descent of regions of strong mean wind shear. The numerical solutions display mean flow accelerations due to Kelvin waves in the equatorial stratosphere, with wave absorption altering the transience mechanism in such significant respects as causing the upper atmospheric mean flow acceleration to be very sensitive to the precise magnitude and distribution of the damping mechanisms. The numerical simulations of transient equatorial waves in the quasi-biennial oscillation are also considered.

Numerical modelling of bifurcation and localisation in cohesive-frictional materials

NASA Astrophysics Data System (ADS)

de Borst, René

1991-12-01

Methods are reviewed for analysing highly localised failure and bifurcation modes in discretised mechanical systems as typically arise in numerical simulations of failure in soils, rocks, metals and concrete. By the example of a plane-strain biaxial test it is shown that strain softening and lack of normality in elasto-plastic constitutive equations and the ensuing loss of ellipticity of the governing field equations cause a pathological mesh dependence of numerical solutions for such problems, thus rendering the results effectively meaningless. The need for introduction of higher-order continuum models is emphasised to remedy this shortcoming of the conventional approach. For one such a continuum model, namely the unconstrained Cosserat continuum, it is demonstrated that meaningful and convergent solutions (in the sense that a finite width of the localisation zone is computed upon mesh refinement) can be obtained.

Numerical method based on the lattice Boltzmann model for the Fisher equation.

PubMed

Yan, Guangwu; Zhang, Jianying; Dong, Yinfeng

2008-06-01

In this paper, a lattice Boltzmann model for the Fisher equation is proposed. First, the Chapman-Enskog expansion and the multiscale time expansion are used to describe higher-order moment of equilibrium distribution functions and a series of partial differential equations in different time scales. Second, the modified partial differential equation of the Fisher equation with the higher-order truncation error is obtained. Third, comparison between numerical results of the lattice Boltzmann models and exact solution is given. The numerical results agree well with the classical ones.

Approximating a retarded-advanced differential equation that models human phonation

NASA Astrophysics Data System (ADS)

Teodoro, M. Filomena

2017-11-01

In [1, 2, 3] we have got the numerical solution of a linear mixed type functional differential equation (MTFDE) introduced initially in [4], considering the autonomous and non-autonomous case by collocation, least squares and finite element methods considering B-splines basis set. The present work introduces a numerical scheme using least squares method (LSM) and Gaussian basis functions to solve numerically a nonlinear mixed type equation with symmetric delay and advance which models human phonation. The preliminary results are promising. We obtain an accuracy comparable with the previous results.

Numerical modeling of rapidly varying flows using HEC-RAS and WSPG models.

PubMed

Rao, Prasada; Hromadka, Theodore V

2016-01-01

The performance of two popular hydraulic models (HEC-RAS and WSPG) for modeling hydraulic jump in an open channel is investigated. The numerical solutions are compared with a new experimental data set obtained for varying channel bottom slopes and flow rates. Both the models satisfactorily predict the flow depths and location of the jump. The end results indicate that the numerical models output is sensitive to the value of chosen roughness coefficient. For this application, WSPG model is easier to implement with few input variables.

Equilibrium paths analysis of materials with rheological properties by using the chaos theory

NASA Astrophysics Data System (ADS)

Bednarek, Paweł; Rządkowski, Jan

2018-01-01

The numerical equilibrium path analysis of the material with random rheological properties by using standard procedures and specialist computer programs was not successful. The proper solution for the analysed heuristic model of the material was obtained on the base of chaos theory elements and neural networks. The paper deals with mathematical reasons of used computer programs and also are elaborated the properties of the attractor used in analysis. There are presented results of conducted numerical analysis both in a numerical and in graphical form for the used procedures.

Numerical investigation of stability of breather-type solutions of the nonlinear Schrödinger equation

NASA Astrophysics Data System (ADS)

Calini, A.; Schober, C. M.

2013-09-01

In this article we present the results of a broad numerical investigation on the stability of breather-type solutions of the nonlinear Schrödinger (NLS) equation, specifically the one- and two-mode breathers for an unstable plane wave, which are frequently used to model rogue waves. The numerical experiments involve large ensembles of perturbed initial data for six typical random perturbations. Ensemble estimates of the "closeness", A(t), of the perturbed solution to an element of the respective unperturbed family indicate that the only neutrally stable breathers are the ones of maximal dimension, that is: given an unstable background with N unstable modes, the only neutrally stable breathers are the N-dimensional ones (obtained as a superimposition of N simple breathers via iterated Backlund transformations). Conversely, breathers which are not fully saturated are sensitive to noisy environments and are unstable. Interestingly, A(t) is smallest for the coalesced two-mode breather indicating the coalesced case may be the most robust two-mode breather in a laboratory setting. The numerical simulations confirm and provide a realistic realization of the stability behavior established analytically by the authors.