A three-dimensional, time-dependent model of Mobile Bay

NASA Technical Reports Server (NTRS)

Pitts, F. H.; Farmer, R. C.

1976-01-01

A three-dimensional, time-variant mathematical model for momentum and mass transport in estuaries was developed and its solution implemented on a digital computer. The mathematical model is based on state and conservation equations applied to turbulent flow of a two-component, incompressible fluid having a free surface. Thus, bouyancy effects caused by density differences between the fresh and salt water, inertia from thare river and tidal currents, and differences in hydrostatic head are taken into account. The conservation equations, which are partial differential equations, are solved numerically by an explicit, one-step finite difference scheme and the solutions displayed numerically and graphically. To test the validity of the model, a specific estuary for which scaled model and experimental field data are available, Mobile Bay, was simulated. Comparisons of velocity, salinity and water level data show that the model is valid and a viable means of simulating the hydrodynamics and mass transport in non-idealized estuaries.

Numerical study of radiometric forces via the direct solution of the Boltzmann kinetic equation

NASA Astrophysics Data System (ADS)

Anikin, Yu. A.

2011-07-01

The two-dimensional rarefied gas motion in a Crookes radiometer and the resulting radiometric forces are studied by numerically solving the Boltzmann kinetic equation. The collision integral is directly evaluated using a projection method, and second-order accurate TVD schemes are used to solve the advection equation. The radiometric forces are found as functions of the Knudsen number and the temperatures, and their spatial distribution is analyzed.

A computer code for three-dimensional incompressible flows using nonorthogonal body-fitted coordinate systems

NASA Technical Reports Server (NTRS)

Chen, Y. S.

1986-01-01

In this report, a numerical method for solving the equations of motion of three-dimensional incompressible flows in nonorthogonal body-fitted coordinate (BFC) systems has been developed. The equations of motion are transformed to a generalized curvilinear coordinate system from which the transformed equations are discretized using finite difference approximations in the transformed domain. The hybrid scheme is used to approximate the convection terms in the governing equations. Solutions of the finite difference equations are obtained iteratively by using a pressure-velocity correction algorithm (SIMPLE-C). Numerical examples of two- and three-dimensional, laminar and turbulent flow problems are employed to evaluate the accuracy and efficiency of the present computer code. The user's guide and computer program listing of the present code are also included.

An Improved Treatment of External Boundary for Three-Dimensional Flow Computations

NASA Technical Reports Server (NTRS)

Tsynkov, Semyon V.; Vatsa, Veer N.

1997-01-01

We present an innovative numerical approach for setting highly accurate nonlocal boundary conditions at the external computational boundaries when calculating three-dimensional compressible viscous flows over finite bodies. The approach is based on application of the difference potentials method by V. S. Ryaben'kii and extends our previous technique developed for the two-dimensional case. The new boundary conditions methodology has been successfully combined with the NASA-developed code TLNS3D and used for the analysis of wing-shaped configurations in subsonic and transonic flow regimes. As demonstrated by the computational experiments, the improved external boundary conditions allow one to greatly reduce the size of the computational domain while still maintaining high accuracy of the numerical solution. Moreover, they may provide for a noticeable speedup of convergence of the multigrid iterations.

High-Order Semi-Discrete Central-Upwind Schemes for Multi-Dimensional Hamilton-Jacobi Equations

NASA Technical Reports Server (NTRS)

Bryson, Steve; Levy, Doron; Biegel, Bran R. (Technical Monitor)

2002-01-01

We present high-order semi-discrete central-upwind numerical schemes for approximating solutions of multi-dimensional Hamilton-Jacobi (HJ) equations. This scheme is based on the use of fifth-order central interpolants like those developed in [1], in fluxes presented in [3]. These interpolants use the weighted essentially nonoscillatory (WENO) approach to avoid spurious oscillations near singularities, and become "central-upwind" in the semi-discrete limit. This scheme provides numerical approximations whose error is as much as an order of magnitude smaller than those in previous WENO-based fifth-order methods [2, 1]. Thee results are discussed via examples in one, two and three dimensions. We also pregnant explicit N-dimensional formulas for the fluxes, discuss their monotonicity and tl!e connection between this method and that in [2].

Analytic solutions for Long's equation and its generalization

NASA Astrophysics Data System (ADS)

Humi, Mayer

2017-12-01

Two-dimensional, steady-state, stratified, isothermal atmospheric flow over topography is governed by Long's equation. Numerical solutions of this equation were derived and used by several authors. In particular, these solutions were applied extensively to analyze the experimental observations of gravity waves. In the first part of this paper we derive an extension of this equation to non-isothermal flows. Then we devise a transformation that simplifies this equation. We show that this simplified equation admits solitonic-type solutions in addition to regular gravity waves. These new analytical solutions provide new insights into the propagation and amplitude of gravity waves over topography.

A Model for Measured Traveling Waves at End-Diastole in Human Heart Wall by Ultrasonic Imaging Method

NASA Astrophysics Data System (ADS)

Bekki, Naoaki; Shintani, Seine A.; Ishiwata, Shin'ichi; Kanai, Hiroshi

2016-04-01

We observe traveling waves, measured by the ultrasonic noninvasive imaging method, in a longitudinal beam direction from the apex to the base side on the interventricular septum (IVS) during the period from the end-diastole to the beginning of systole for a healthy human heart wall. We present a possible phenomenological model to explain part of one-dimensional cardiac behaviors for the observed traveling waves around the time of R-wave of echocardiography (ECG) in the human heart. Although the observed two-dimensional patterns of traveling waves are extremely complex and no one knows yet the exact solutions for the traveling homoclinic plane wave in the one-dimensional complex Ginzburg-Landau equation (CGLE), we numerically find that part of the one-dimensional homoclinic dynamics of the phase and amplitude patterns in the observed traveling waves is similar to that of the numerical homoclinic plane-wave solutions in the CGLE with periodic boundary condition in a certain parameter space. It is suggested that part of the cardiac dynamics of the traveling waves on the IVS can be qualitatively described by the CGLE model as a paradigm for understanding biophysical nonlinear phenomena.

A two-dimensional analytical model for groundwater flow in a leaky aquifer extending finite distance under the estuary

NASA Astrophysics Data System (ADS)

Chuang, Mo-Hsiung; Hung, Chi-Tung; -Yen Lin, Wen; Ma, Kuo-chen

2017-04-01

In recent years, cities and industries in the vicinity of the estuarine region have developed rapidly, resulting in a sharp increase in the population concerned. The increasing demand for human activities, agriculture irrigation, and aquaculture relies on massive pumping of water in estuarine area. Since the 1950s, numerous studies have focused on the effects of tidal fluctuations on groundwater flow in the estuarine area. Tide-induced head fluctuation in a two-dimensional estuarine aquifer system is complicated and rather important in dealing with many groundwater management or remediation problems. The conceptual model of the aquifer system considered is multi-layered with estuarine bank and the leaky aquifer extend finite distance under the estuary. The solution of the model describing the groundwater head distribution in such an estuarine aquifer system and subject to the tidal fluctuation effects from estuarine river is developed based on the method of separation of variables along with river boundary. The solutions by Sun (Sun H. A two-dimensional analytical solution of groundwater response to tidal loading in an estuary, Water Resour. Res. 1997; 33:1429-35) as well as Tang and Jiao (Tang Z. and J. J. Jiao, A two-dimensional analytical solution for groundwater flow in a leaky confined aquifer system near open tidal water, Hydrological Processes, 2001; 15: 573-585) can be shown to be special cases of the present solution. On the basis of the analytical solution, the groundwater head distribution in response to estuarine boundary is examined and the influences of leakage, hydraulic parameters, and loading effect on the groundwater head fluctuation due to tide are investigated and discussed. KEYWORDS: analytical model, estuarine river, groundwater fluctuation, leaky aquifer.

Development of a three dimensional numerical water quality model for continental shelf applications

NASA Technical Reports Server (NTRS)

Spaulding, M.; Hunter, D.

1975-01-01

A model to predict the distribution of water quality parameters in three dimensions was developed. The mass transport equation was solved using a non-dimensional vertical axis and an alternating-direction-implicit finite difference technique. The reaction kinetics of the constituents were incorporated into a matrix method which permits computation of the interactions of multiple constituents. Methods for the computation of dispersion coefficients and coliform bacteria decay rates were determined. Numerical investigations of dispersive and dissipative effects showed that the three-dimensional model performs as predicted by analysis of simpler cases. The model was then applied to a two dimensional vertically averaged tidal dynamics model for the Providence River. It was also extended to a steady state application by replacing the time step with an iteration sequence. This modification was verified by comparison to analytical solutions and applied to a river confluence situation.

Modeling and numerical simulations of growth and morphologies of three dimensional aggregated silver films

NASA Astrophysics Data System (ADS)

Davis, L. J.; Boggess, M.; Kodpuak, E.; Deutsch, M.

2012-11-01

We report on a model for the deposition of three dimensional, aggregated nanocrystalline silver films, and an efficient numerical simulation method developed for visualizing such structures. We compare our results to a model system comprising chemically deposited silver films with morphologies ranging from dilute, uniform distributions of nanoparticles to highly porous aggregated networks. Disordered silver films grown in solution on silica substrates are characterized using digital image analysis of high resolution scanning electron micrographs. While the latter technique provides little volume information, plane-projected (two dimensional) island structure and surface coverage may be reliably determined. Three parameters governing film growth are evaluated using these data and used as inputs for the deposition model, greatly reducing computing requirements while still providing direct access to the complete (bulk) structure of the films throughout the growth process. We also show how valuable three dimensional characteristics of the deposited materials can be extracted using the simulated structures.

Calculation of flow about two-dimensional bodies by means of the velocity-vorticity formulation on a staggered grid

NASA Technical Reports Server (NTRS)

Stremel, Paul M.

1991-01-01

A method for calculating the incompressible viscous flow about two-dimensional bodies, utilizing the velocity-vorticity form of the Navier-Stokes equations using a staggered-grid formulation is presented. The solution is obtained by employing an alternative-direction implicit method for the solution of the block tridiagonal matrix resulting from the finite-difference representation of the governing equations. The boundary vorticity and the conservation of mass are calculated implicitly as a part of the solution. The mass conservation is calculated to machine zero for the duration of the computation. Calculations for the flow about a circular cylinder, a 2-pct thick flat plate at 90-deg incidence, an elliptic cylinder at 45-deg incidence, and a NACA 0012, with and without a deflected flap, at - 90-deg incidence are performed and compared with the results of other numerical investigations.

Nonequilibrium recombination after a curved shock wave

NASA Astrophysics Data System (ADS)

Wen, Chihyung; Hornung, Hans

2010-02-01

The effect of nonequilibrium recombination after a curved two-dimensional shock wave in a hypervelocity dissociating flow of an inviscid Lighthill-Freeman gas is considered. An analytical solution is obtained with the effective shock values derived by Hornung (1976) [5] and the assumption that the flow is ‘quasi-frozen’ after a thin dissociating layer near the shock. The solution gives the expression of dissociation fraction as a function of temperature on a streamline. A rule of thumb can then be provided to check the validity of binary scaling for experimental conditions and a tool to determine the limiting streamline that delineates the validity zone of binary scaling. The effects on the nonequilibrium chemical reaction of the large difference in free stream temperature between free-piston shock tunnel and equivalent flight conditions are discussed. Numerical examples are presented and the results are compared with solutions obtained with two-dimensional Euler equations using the code of Candler (1988) [10].

Conjugate-gradient optimization method for orbital-free density functional calculations.

PubMed

Jiang, Hong; Yang, Weitao

2004-08-01

Orbital-free density functional theory as an extension of traditional Thomas-Fermi theory has attracted a lot of interest in the past decade because of developments in both more accurate kinetic energy functionals and highly efficient numerical methodology. In this paper, we developed a conjugate-gradient method for the numerical solution of spin-dependent extended Thomas-Fermi equation by incorporating techniques previously used in Kohn-Sham calculations. The key ingredient of the method is an approximate line-search scheme and a collective treatment of two spin densities in the case of spin-dependent extended Thomas-Fermi problem. Test calculations for a quartic two-dimensional quantum dot system and a three-dimensional sodium cluster Na216 with a local pseudopotential demonstrate that the method is accurate and efficient. (c) 2004 American Institute of Physics.

Gravity-Driven Thin Film Flow of an Ellis Fluid.

PubMed

Kheyfets, Vitaly O; Kieweg, Sarah L

2013-12-01

The thin film lubrication approximation has been studied extensively for moving contact lines of Newtonian fluids. However, many industrial and biological applications of the thin film equation involve shear-thinning fluids, which often also exhibit a Newtonian plateau at low shear. This study presents new numerical simulations of the three-dimensional (i.e. two-dimensional spreading), constant-volume, gravity-driven, free surface flow of an Ellis fluid. The numerical solution was validated with a new similarity solution, compared to previous experiments, and then used in a parametric study. The parametric study centered around rheological data for an example biological application of thin film flow: topical drug delivery of anti-HIV microbicide formulations, e.g. hydroxyethylcellulose (HEC) polymer solutions. The parametric study evaluated how spreading length and front velocity saturation depend on Ellis parameters. A lower concentration polymer solution with smaller zero shear viscosity ( η 0 ), τ 1/2 , and λ values spread further. However, when comparing any two fluids with any possible combinations of Ellis parameters, the impact of changing one parameter on spreading length depends on the direction and magnitude of changes in the other two parameters. In addition, the isolated effect of the shear-thinning parameter, λ , on the front velocity saturation depended on τ 1/2 . This study highlighted the relative effects of the individual Ellis parameters, and showed that the shear rates in this flow were in both the shear-thinning and plateau regions of rheological behavior, emphasizing the importance of characterizing the full range of shear-rates in rheological measurements. The validated numerical model and parametric study provides a useful tool for future steps to optimize flow of a fluid with rheological behavior well-described by the Ellis constitutive model, in a range of industrial and biological applications.

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.

Perturbational blowup solutions to the compressible Euler equations with damping.

PubMed

Cheung, Ka Luen

2016-01-01

The N-dimensional isentropic compressible Euler system with a damping term is one of the most fundamental equations in fluid dynamics. Since it does not have a general solution in a closed form for arbitrary well-posed initial value problems. Constructing exact solutions to the system is a useful way to obtain important information on the properties of its solutions. In this article, we construct two families of exact solutions for the one-dimensional isentropic compressible Euler equations with damping by the perturbational method. The two families of exact solutions found include the cases [Formula: see text] and [Formula: see text], where [Formula: see text] is the adiabatic constant. With analysis of the key ordinary differential equation, we show that the classes of solutions include both blowup type and global existence type when the parameters are suitably chosen. Moreover, in the blowup cases, we show that the singularities are of essential type in the sense that they cannot be smoothed by redefining values at the odd points. The two families of exact solutions obtained in this paper can be useful to study of related numerical methods and algorithms such as the finite difference method, the finite element method and the finite volume method that are applied by scientists to simulate the fluids for applications.

Multigrid for hypersonic viscous two- and three-dimensional flows

NASA Technical Reports Server (NTRS)

Turkel, E.; Swanson, R. C.; Vatsa, V. N.; White, J. A.

1991-01-01

The use of a multigrid method with central differencing to solve the Navier-Stokes equations for hypersonic flows is considered. The time dependent form of the equations is integrated with an explicit Runge-Kutta scheme accelerated by local time stepping and implicit residual smoothing. Variable coefficients are developed for the implicit process that removes the diffusion limit on the time step, producing significant improvement in convergence. A numerical dissipation formulation that provides good shock capturing capability for hypersonic flows is presented. This formulation is shown to be a crucial aspect of the multigrid method. Solutions are given for two-dimensional viscous flow over a NACA 0012 airfoil and three-dimensional flow over a blunt biconic.

A MacCormack-TVD finite difference method to simulate the mass flow in mountainous terrain with variable computational domain

NASA Astrophysics Data System (ADS)

Ouyang, Chaojun; He, Siming; Xu, Qiang; Luo, Yu; Zhang, Wencheng

2013-03-01

A two-dimensional mountainous mass flow dynamic procedure solver (Massflow-2D) using the MacCormack-TVD finite difference scheme is proposed. The solver is implemented in Matlab on structured meshes with variable computational domain. To verify the model, a variety of numerical test scenarios, namely, the classical one-dimensional and two-dimensional dam break, the landslide in Hong Kong in 1993 and the Nora debris flow in the Italian Alps in 2000, are executed, and the model outputs are compared with published results. It is established that the model predictions agree well with both the analytical solution as well as the field observations.

Numerical method for solution of systems of non-stationary spatially one-dimensional nonlinear differential equations

NASA Technical Reports Server (NTRS)

Morozov, S. K.; Krasitskiy, O. P.

1978-01-01

A computational scheme and a standard program is proposed for solving systems of nonstationary spatially one-dimensional nonlinear differential equations using Newton's method. The proposed scheme is universal in its applicability and its reduces to a minimum the work of programming. The program is written in the FORTRAN language and can be used without change on electronic computers of type YeS and BESM-6. The standard program described permits the identification of nonstationary (or stationary) solutions to systems of spatially one-dimensional nonlinear (or linear) partial differential equations. The proposed method may be used to solve a series of geophysical problems which take chemical reactions, diffusion, and heat conductivity into account, to evaluate nonstationary thermal fields in two-dimensional structures when in one of the geometrical directions it can take a small number of discrete levels, and to solve problems in nonstationary gas dynamics.

Shock formation in the dispersionless Kadomtsev-Petviashvili equation

NASA Astrophysics Data System (ADS)

Grava, T.; Klein, C.; Eggers, J.

2016-04-01

The dispersionless Kadomtsev-Petviashvili (dKP) equation {{≤ft({{u}t}+u{{u}x}\\right)}x}={{u}yy} is one of the simplest nonlinear wave equations describing two-dimensional shocks. To solve the dKP equation numerically we use a coordinate transformation inspired by the method of characteristics for the one-dimensional Hopf equation {{u}t}+u{{u}x}=0 . We show numerically that the solutions to the transformed equation stays regular for longer times than the solution of the dKP equation. This permits us to extend the dKP solution as the graph of a multivalued function beyond the critical time when the gradients blow up. This overturned solution is multivalued in a lip shape region in the (x, y) plane, where the solution of the dKP equation exists in a weak sense only, and a shock front develops. A local expansion reveals the universal scaling structure of the shock, which after a suitable change of coordinates corresponds to a generic cusp catastrophe. We provide a heuristic derivation of the shock front position near the critical point for the solution of the dKP equation, and study the solution of the dKP equation when a small amount of dissipation is added. Using multiple-scale analysis, we show that in the limit of small dissipation and near the critical point of the dKP solution, the solution of the dissipative dKP equation converges to a Pearcey integral. We test and illustrate our results by detailed comparisons with numerical simulations of both the regularized equation, the dKP equation, and the asymptotic description given in terms of the Pearcey integral.

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.

Numerical Simulation of Interaction of Human Vocal Folds and Fluid Flow

NASA Astrophysics Data System (ADS)

Kosík, A.; Feistauer, M.; Horáček, J.; Sváček, P.

Our goal is to simulate airflow in human vocal folds and their flow-induced vibrations. We consider two-dimensional viscous incompressible flow in a time-dependent domain. The fluid flow is described by the Navier-Stokes equations in the arbitrary Lagrangian-Eulerian formulation. The flow problem is coupled with the elastic behaviour of the solid bodies. The developed solution of the coupled problem based on the finite element method is demonstrated by numerical experiments.

The Proteus Navier-Stokes code

NASA Technical Reports Server (NTRS)

Towne, Charles E.; Bui, Trong T.; Cavicchi, Richard H.; Conley, Julianne M.; Molls, Frank B.; Schwab, John R.

1992-01-01

An effort is currently underway at NASA Lewis to develop two- and three-dimensional Navier-Stokes codes, called Proteus, for aerospace propulsion applications. The emphasis in the development of Proteus is not algorithm development or research on numerical methods, but rather the development of the code itself. The objective is to develop codes that are user-oriented, easily-modified, and well-documented. Well-proven, state-of-the-art solution algorithms are being used. Code readability, documentation (both internal and external), and validation are being emphasized. This paper is a status report on the Proteus development effort. The analysis and solution procedure are described briefly, and the various features in the code are summarized. The results from some of the validation cases that have been run are presented for both the two- and three-dimensional codes.

Characterization of the heavily doped emitter and junction regions of silicon solar cells using an electron beam

NASA Technical Reports Server (NTRS)

Luke, K. L.; Cheng, L.-J.

1986-01-01

Heavily doped emitter and junction regions of silicon solar cells are investigated by means of the electron-beam-induced-current (EBIC) technique. Although the experimental EBIC data are collected under three-dimensional conditions, it is analytically demonstrated with two numerical examples that the solutions obtained with one-dimensional numerical modeling are adequate. EBIC data for bare and oxide-covered emitter surfaces are compared with theory. The improvement in collection efficiency when an emitter surface is covered with a 100-A SiO2 film varies with beam energy; for a cell with a junction depth of 0.35 microns, the improvement is about 54 percent at 2 keV.

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.

Numerical Solutions of the Complete Navier-Stokes Equations

NASA Technical Reports Server (NTRS)

Robinson, David F.; Hassan, H. A.

1997-01-01

This report details the development of a new two-equation turbulence closure model based on the exact turbulent kinetic energy k and the variance of vorticity, zeta. The model, which is applicable to three dimensional flowfields, employs one set of model constants and does not use damping or wall functions, or geometric factors.

Testing density-dependent groundwater models: Two-dimensional steady state unstable convection in infinite, finite and inclined porous layers

USGS Publications Warehouse

Weatherill, D.; Simmons, C.T.; Voss, C.I.; Robinson, N.I.

2004-01-01

This study proposes the use of several problems of unstable steady state convection with variable fluid density in a porous layer of infinite horizontal extent as two-dimensional (2-D) test cases for density-dependent groundwater flow and solute transport simulators. Unlike existing density-dependent model benchmarks, these problems have well-defined stability criteria that are determined analytically. These analytical stability indicators can be compared with numerical model results to test the ability of a code to accurately simulate buoyancy driven flow and diffusion. The basic analytical solution is for a horizontally infinite fluid-filled porous layer in which fluid density decreases with depth. The proposed test problems include unstable convection in an infinite horizontal box, in a finite horizontal box, and in an infinite inclined box. A dimensionless Rayleigh number incorporating properties of the fluid and the porous media determines the stability of the layer in each case. Testing the ability of numerical codes to match both the critical Rayleigh number at which convection occurs and the wavelength of convection cells is an addition to the benchmark problems currently in use. The proposed test problems are modelled in 2-D using the SUTRA [SUTRA-A model for saturated-unsaturated variable-density ground-water flow with solute or energy transport. US Geological Survey Water-Resources Investigations Report, 02-4231, 2002. 250 p] density-dependent groundwater flow and solute transport code. For the case of an infinite horizontal box, SUTRA results show a distinct change from stable to unstable behaviour around the theoretical critical Rayleigh number of 4??2 and the simulated wavelength of unstable convection agrees with that predicted by the analytical solution. The effects of finite layer aspect ratio and inclination on stability indicators are also tested and numerical results are in excellent agreement with theoretical stability criteria and with numerical results previously reported in traditional fluid mechanics literature. ?? 2004 Elsevier Ltd. All rights reserved.

An Assessment of Artificial Compressibility and Pressure Projection Methods for Incompressible Flow Simulations

NASA Technical Reports Server (NTRS)

Kwak, Dochan; Kiris, C.; Smith, Charles A. (Technical Monitor)

1998-01-01

Performance of the two commonly used numerical procedures, one based on artificial compressibility method and the other pressure projection method, are compared. These formulations are selected primarily because they are designed for three-dimensional applications. The computational procedures are compared by obtaining steady state solutions of a wake vortex and unsteady solutions of a curved duct flow. For steady computations, artificial compressibility was very efficient in terms of computing time and robustness. For an unsteady flow which requires small physical time step, pressure projection method was found to be computationally more efficient than an artificial compressibility method. This comparison is intended to give some basis for selecting a method or a flow solution code for large three-dimensional applications where computing resources become a critical issue.

Analytical and numerical solutions of the potential and electric field generated by different electrode arrays in a tumor tissue under electrotherapy

PubMed Central

2011-01-01

Background Electrotherapy is a relatively well established and efficient method of tumor treatment. In this paper we focus on analytical and numerical calculations of the potential and electric field distributions inside a tumor tissue in a two-dimensional model (2D-model) generated by means of electrode arrays with shapes of different conic sections (ellipse, parabola and hyperbola). Methods Analytical calculations of the potential and electric field distributions based on 2D-models for different electrode arrays are performed by solving the Laplace equation, meanwhile the numerical solution is solved by means of finite element method in two dimensions. Results Both analytical and numerical solutions reveal significant differences between the electric field distributions generated by electrode arrays with shapes of circle and different conic sections (elliptic, parabolic and hyperbolic). Electrode arrays with circular, elliptical and hyperbolic shapes have the advantage of concentrating the electric field lines in the tumor. Conclusion The mathematical approach presented in this study provides a useful tool for the design of electrode arrays with different shapes of conic sections by means of the use of the unifying principle. At the same time, we verify the good correspondence between the analytical and numerical solutions for the potential and electric field distributions generated by the electrode array with different conic sections. PMID:21943385

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 two-dimensional, TVD numerical scheme for inviscid, high Mach number flows in chemical equilibrium

NASA Technical Reports Server (NTRS)

Eberhardt, S.; Palmer, G.

1986-01-01

A new algorithm has been developed for hypervelocity flows in chemical equilibrium. Solutions have been achieved for Mach numbers up to 15 with no adverse effect on convergence. Two methods of coupling an equilibrium chemistry package have been tested, with the simpler method proving to be more robust. Improvements in boundary conditions are still required for a production-quality code.

An analytic-geometric model of the effect of spherically distributed injection errors for Galileo and Ulysses spacecraft - The multi-stage problem

NASA Technical Reports Server (NTRS)

Longuski, James M.; Mcronald, Angus D.

1988-01-01

In previous work the problem of injecting the Galileo and Ulysses spacecraft from low earth orbit into their respective interplanetary trajectories has been discussed for the single stage (Centaur) vehicle. The central issue, in the event of spherically distributed injection errors, is what happens to the vehicle? The difficulties addressed in this paper involve the multi-stage problem since both Galileo and Ulysses will be utilizing the two-stage IUS system. Ulysses will also include a third stage: the PAM-S. The solution is expressed in terms of probabilities for total percentage of escape, orbit decay and reentry trajectories. Analytic solutions are found for Hill's Equations of Relative Motion (more recently called Clohessy-Wiltshire Equations) for multi-stage injections. These solutions are interpreted geometrically on the injection sphere. The analytic-geometric models compare well with numerical solutions, provide insight into the behavior of trajectories mapped on the injection sphere and simplify the numerical two-dimensional search for trajectory families.

Swimming in a two-dimensional Brinkman fluid: Computational modeling and regularized solutions

NASA Astrophysics Data System (ADS)

Leiderman, Karin; Olson, Sarah D.

2016-02-01

The incompressible Brinkman equation represents the homogenized fluid flow past obstacles that comprise a small volume fraction. In nondimensional form, the Brinkman equation can be characterized by a single parameter that represents the friction or resistance due to the obstacles. In this work, we derive an exact fundamental solution for 2D Brinkman flow driven by a regularized point force and describe the numerical method to use it in practice. To test our solution and method, we compare numerical results with an analytic solution of a stationary cylinder in a uniform Brinkman flow. Our method is also compared to asymptotic theory; for an infinite-length, undulating sheet of small amplitude, we recover an increasing swimming speed as the resistance is increased. With this computational framework, we study a model swimmer of finite length and observe an enhancement in propulsion and efficiency for small to moderate resistance. Finally, we study the interaction of two swimmers where attraction does not occur when the initial separation distance is larger than the screening length.

A sharp interface method for compressible liquid–vapor flow with phase transition and surface tension

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

Fechter, Stefan, E-mail: stefan.fechter@iag.uni-stuttgart.de; Munz, Claus-Dieter, E-mail: munz@iag.uni-stuttgart.de; Rohde, Christian, E-mail: Christian.Rohde@mathematik.uni-stuttgart.de

The numerical approximation of non-isothermal liquid–vapor flow within the compressible regime is a difficult task because complex physical effects at the phase interfaces can govern the global flow behavior. We present a sharp interface approach which treats the interface as a shock-wave like discontinuity. Any mixing of fluid phases is avoided by using the flow solver in the bulk regions only, and a ghost-fluid approach close to the interface. The coupling states for the numerical solution in the bulk regions are determined by the solution of local two-phase Riemann problems across the interface. The Riemann solution accounts for the relevantmore » physics by enforcing appropriate jump conditions at the phase boundary. A wide variety of interface effects can be handled in a thermodynamically consistent way. This includes surface tension or mass/energy transfer by phase transition. Moreover, the local normal speed of the interface, which is needed to calculate the time evolution of the interface, is given by the Riemann solution. The interface tracking itself is based on a level-set method. The focus in this paper is the description of the two-phase Riemann solver and its usage within the sharp interface approach. One-dimensional problems are selected to validate the approach. Finally, the three-dimensional simulation of a wobbling droplet and a shock droplet interaction in two dimensions are shown. In both problems phase transition and surface tension determine the global bulk behavior.« less

Asymptotic structure and similarity solutions for three-dimensional turbulent boundary layers

NASA Technical Reports Server (NTRS)

Degani, A. T.; Walker, J. D. A.

1989-01-01

The asymptotic structure of the three-dimensional turbulent boundary layer is investigated in the limit of large Reynolds numbers. A self-consistent, but relatively complex, two-layer structure exists and the simplest situation, corresponding to a plane of symmetry, is considered in this paper as a first step. The adjustment of the streamwise velocity to relative rest, through an outer defect layer and then an inner wall layer, is similar to that in two-dimensional flow. The adjustment of the cross-streamwise velocity is more complicated and it is shown that two terms in the expansion are required to obtain useful results, and in particular to obtain the velocity skew angle at the wall near the symmetry plane. The conditions under which self-similarity is achieved near a plane of symmetry are investigated. A set of ordinary differential equations is developed which describe the streamwise and cross-streamwise velocities near a plane of symmetry in a self-similar flow through two orders of magnitude. Calculated numerical solutions of these equations yield trends which are consistent with experimental observations.

Fractional calculus phenomenology in two-dimensional plasma models

NASA Astrophysics Data System (ADS)

Gustafson, Kyle; Del Castillo Negrete, Diego; Dorland, Bill

2006-10-01

Transport processes in confined plasmas for fusion experiments, such as ITER, are not well-understood at the basic level of fully nonlinear, three-dimensional kinetic physics. Turbulent transport is invoked to describe the observed levels in tokamaks, which are orders of magnitude greater than the theoretical predictions. Recent results show the ability of a non-diffusive transport model to describe numerical observations of turbulent transport. For example, resistive MHD modeling of tracer particle transport in pressure-gradient driven turbulence for a three-dimensional plasma reveals that the superdiffusive (2̂˜t^α where α> 1) radial transport in this system is described quantitatively by a fractional diffusion equation Fractional calculus is a generalization involving integro-differential operators, which naturally describe non-local behaviors. Our previous work showed the quantitative agreement of special fractional diffusion equation solutions with numerical tracer particle flows in time-dependent linearized dynamics of the Hasegawa-Mima equation (for poloidal transport in a two-dimensional cold-ion plasma). In pursuit of a fractional diffusion model for transport in a gyrokinetic plasma, we now present numerical results from tracer particle transport in the nonlinear Hasegawa-Mima equation and a planar gyrokinetic model. Finite Larmor radius effects will be discussed. D. del Castillo Negrete, et al, Phys. Rev. Lett. 94, 065003 (2005).

Dynamic colloidal assembly pathways via low dimensional models

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

Yang, Yuguang; Bevan, Michael A., E-mail: mabevan@jhu.edu; Thyagarajan, Raghuram

2016-05-28

Here we construct a low-dimensional Smoluchowski model for electric field mediated colloidal crystallization using Brownian dynamic simulations, which were previously matched to experiments. Diffusion mapping is used to infer dimensionality and confirm the use of two order parameters, one for degree of condensation and one for global crystallinity. Free energy and diffusivity landscapes are obtained as the coefficients of a low-dimensional Smoluchowski equation to capture the thermodynamics and kinetics of microstructure evolution. The resulting low-dimensional model quantitatively captures the dynamics of different assembly pathways between fluid, polycrystal, and single crystals states, in agreement with the full N-dimensional data as characterizedmore » by first passage time distributions. Numerical solution of the low-dimensional Smoluchowski equation reveals statistical properties of the dynamic evolution of states vs. applied field amplitude and system size. The low-dimensional Smoluchowski equation and associated landscapes calculated here can serve as models for predictive control of electric field mediated assembly of colloidal ensembles into two-dimensional crystalline objects.« less

Solving the Helmholtz equation in conformal mapped ARROW structures using homotopy perturbation method.

PubMed

Reck, Kasper; Thomsen, Erik V; Hansen, Ole

2011-01-31

The scalar wave equation, or Helmholtz equation, describes within a certain approximation the electromagnetic field distribution in a given system. In this paper we show how to solve the Helmholtz equation in complex geometries using conformal mapping and the homotopy perturbation method. The solution of the mapped Helmholtz equation is found by solving an infinite series of Poisson equations using two dimensional Fourier series. The solution is entirely based on analytical expressions and is not mesh dependent. The analytical results are compared to a numerical (finite element method) solution.

Fluid flow and fuel-air mixing in a motored two-dimensional Wankel rotary engine

NASA Technical Reports Server (NTRS)

Shih, T. I.-P.; Nguyen, H. L.; Stegeman, J.

1986-01-01

The implicit-factored method of Beam and Warming was employed to obtain numerical solutions to the conservation equations of mass, species, momentum, and energy to study the unsteady, multidimensional flow and mixing of fuel and air inside the combustion chambers of a two-dimensional Wankel rotary engine under motored conditions. The effects of the following engine design and operating parameters on fluid flow and fuel-air mixing during the intake and compression cycles were studied: engine speed, angle of gaseous fuel injection during compression cycle, and speed of the fuel leaving fuel injector.

Fluid flow and fuel-air mixing in a motored two-dimensional Wankel rotary engine

NASA Astrophysics Data System (ADS)

Shih, T. I.-P.; Nguyen, H. L.; Stegeman, J.

1986-06-01

The implicit-factored method of Beam and Warming was employed to obtain numerical solutions to the conservation equations of mass, species, momentum, and energy to study the unsteady, multidimensional flow and mixing of fuel and air inside the combustion chambers of a two-dimensional Wankel rotary engine under motored conditions. The effects of the following engine design and operating parameters on fluid flow and fuel-air mixing during the intake and compression cycles were studied: engine speed, angle of gaseous fuel injection during compression cycle, and speed of the fuel leaving fuel injector.

User's guide to the NOZL3D and NOZLIC computer programs

NASA Technical Reports Server (NTRS)

Thomas, P. D.

1980-01-01

Complete FORTRAN listings and running instructions are given for a set of computer programs that perform an implicit numerical solution to the unsteady Navier-Stokes equations to predict the flow characteristics and performance of nonaxisymmetric nozzles. The set includes the NOZL3D program, which performs the flow computations; the NOZLIC program, which sets up the flow field initial conditions for general nozzle configurations, and also generates the computational grid for simple two dimensional and axisymmetric configurations; and the RGRIDD program, which generates the computational grid for complicated three dimensional configurations. The programs are designed specifically for the NASA-Langley CYBER 175 computer, and employ auxiliary disk files for primary data storage. Input instructions and computed results are given for four test cases that include two dimensional, three dimensional, and axisymmetric configurations.

Investigation of the complex electroviscous effects on electrolyte (single and multiphase) flow in porous medi.

NASA Astrophysics Data System (ADS)

Bolet, A. J. S.; Linga, G.; Mathiesen, J.

2017-12-01

Surface charge is an important control parameter for wall-bounded flow of electrolyte solution. The electroviscous effect has been studied theoretically in model geometries such as infinite capillaries. However, in more complex geometries a quantification of the electroviscous effect is a non-trival task due to strong non-linarites of the underlying equations. In general, one has to rely on numerical methods. Here we present numerical studies of the full three-dimensional steady state Stokes-Poisson-Nernst-Planck problem in order to model electrolyte transport in artificial porous samples. The simulations are performed using the finite element method. From the simulation, we quantity how the electroviscous effect changes the general flow permeability in complex three-dimensional porous media. The porous media we consider are mostly generated artificially by connecting randomly dispersed cylindrical pores. Furthermore, we present results of electric driven two-phase immiscible flow in two dimensions. The simulations are performed by augmenting the above equations with a phase field model to handle and track the interaction between the two fluids (using parameters corresponding to oil-water interfaces, where oil non-polar). In particular, we consider the electro-osmotic effect on imbibition due to charged walls and electrolyte-solution.

Simulation of confined magnetohydrodynamic flows with Dirichlet boundary conditions using a pseudo-spectral method with volume penalization

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

Morales, Jorge A.; Leroy, Matthieu; Bos, Wouter J.T.

A volume penalization approach to simulate magnetohydrodynamic (MHD) flows in confined domains is presented. Here the incompressible visco-resistive MHD equations are solved using parallel pseudo-spectral solvers in Cartesian geometries. The volume penalization technique is an immersed boundary method which is characterized by a high flexibility for the geometry of the considered flow. In the present case, it allows to use other than periodic boundary conditions in a Fourier pseudo-spectral approach. The numerical method is validated and its convergence is assessed for two- and three-dimensional hydrodynamic (HD) and MHD flows, by comparing the numerical results with results from literature and analyticalmore » solutions. The test cases considered are two-dimensional Taylor–Couette flow, the z-pinch configuration, three dimensional Orszag–Tang flow, Ohmic-decay in a periodic cylinder, three-dimensional Taylor–Couette flow with and without axial magnetic field and three-dimensional Hartmann-instabilities in a cylinder with an imposed helical magnetic field. Finally, we present a magnetohydrodynamic flow simulation in toroidal geometry with non-symmetric cross section and imposing a helical magnetic field to illustrate the potential of the method.« less

Methods in the study of discrete upper hybrid waves

NASA Astrophysics Data System (ADS)

Yoon, P. H.; Ye, S.; Labelle, J.; Weatherwax, A. T.; Menietti, J. D.

2007-11-01

Naturally occurring plasma waves characterized by fine frequency structure or discrete spectrum, detected by satellite, rocket-borne instruments, or ground-based receivers, can be interpreted as eigenmodes excited and trapped in field-aligned density structures. This paper overviews various theoretical methods to study such phenomena for a one-dimensional (1-D) density structure. Among the various methods are parabolic approximation, eikonal matching, eigenfunction matching, and full numerical solution based upon shooting method. Various approaches are compared against the full numerical solution. Among the analytic methods it is found that the eigenfunction matching technique best approximates the actual numerical solution. The analysis is further extended to 2-D geometry. A detailed comparative analysis between the eigenfunction matching and fully numerical methods is carried out for the 2-D case. Although in general the two methods compare favorably, significant differences are also found such that for application to actual observations it is prudent to employ the fully numerical method. Application of the methods developed in the present paper to actual geophysical problems will be given in a companion paper.

The Numerical Analysis of a Turbulent Compressible Jet. Degree awarded by Ohio State Univ., 2000

NASA Technical Reports Server (NTRS)

DeBonis, James R.

2001-01-01

A numerical method to simulate high Reynolds number jet flows was formulated and applied to gain a better understanding of the flow physics. Large-eddy simulation was chosen as the most promising approach to model the turbulent structures due to its compromise between accuracy and computational expense. The filtered Navier-Stokes equations were developed including a total energy form of the energy equation. Subgrid scale models for the momentum and energy equations were adapted from compressible forms of Smagorinsky's original model. The effect of using disparate temporal and spatial accuracy in a numerical scheme was discovered through one-dimensional model problems and a new uniformly fourth-order accurate numerical method was developed. Results from two- and three-dimensional validation exercises show that the code accurately reproduces both viscous and inviscid flows. Numerous axisymmetric jet simulations were performed to investigate the effect of grid resolution, numerical scheme, exit boundary conditions and subgrid scale modeling on the solution and the results were used to guide the three-dimensional calculations. Three-dimensional calculations of a Mach 1.4 jet showed that this LES simulation accurately captures the physics of the turbulent flow. The agreement with experimental data was relatively good and is much better than results in the current literature. Turbulent intensities indicate that the turbulent structures at this level of modeling are not isotropic and this information could lend itself to the development of improved subgrid scale models for LES and turbulence models for RANS simulations. A two point correlation technique was used to quantify the turbulent structures. Two point space correlations were used to obtain a measure of the integral length scale, which proved to be approximately 1/2 D(sub j). Two point space-time correlations were used to obtain the convection velocity for the turbulent structures. This velocity ranged from 0.57 to 0.71 U(sub j).

On a fourth order accurate implicit finite difference scheme for hyperbolic conservation laws. I - Nonstiff strongly dynamic problems

NASA Technical Reports Server (NTRS)

Harten, A.; Tal-Ezer, H.

1981-01-01

An implicit finite difference method of fourth order accuracy in space and time is introduced for the numerical solution of one-dimensional systems of hyperbolic conservation laws. The basic form of the method is a two-level scheme which is unconditionally stable and nondissipative. The scheme uses only three mesh points at level t and three mesh points at level t + delta t. The dissipative version of the basic method given is conditionally stable under the CFL (Courant-Friedrichs-Lewy) condition. This version is particularly useful for the numerical solution of problems with strong but nonstiff dynamic features, where the CFL restriction is reasonable on accuracy grounds. Numerical results are provided to illustrate properties of the proposed method.

New analytical solutions to the two-phase water faucet problem

DOE PAGES

Zou, Ling; Zhao, Haihua; Zhang, Hongbin

2016-06-17

Here, the one-dimensional water faucet problem is one of the classical benchmark problems originally proposed by Ransom to study the two-fluid two-phase flow model. With certain simplifications, such as massless gas phase and no wall and interfacial frictions, analytical solutions had been previously obtained for the transient liquid velocity and void fraction distribution. The water faucet problem and its analytical solutions have been widely used for the purposes of code assessment, benchmark and numerical verifications. In our previous study, the Ransom’s solutions were used for the mesh convergence study of a high-resolution spatial discretization scheme. It was found that, atmore » the steady state, an anticipated second-order spatial accuracy could not be achieved, when compared to the existing Ransom’s analytical solutions. A further investigation showed that the existing analytical solutions do not actually satisfy the commonly used two-fluid single-pressure two-phase flow equations. In this work, we present a new set of analytical solutions of the water faucet problem at the steady state, considering the gas phase density’s effect on pressure distribution. This new set of analytical solutions are used for mesh convergence studies, from which anticipated second-order of accuracy is achieved for the 2nd order spatial discretization scheme. In addition, extended Ransom’s transient solutions for the gas phase velocity and pressure are derived, with the assumption of decoupled liquid and gas pressures. Numerical verifications on the extended Ransom’s solutions are also presented.« less

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.

A two-dimensional composite grid numerical model based on the reduced system for oceanography

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

Xie, Y.F.; Browning, G.L.; Chesshire, G.

The proper mathematical limit of a hyperbolic system with multiple time scales, the reduced system, is a system that contains no high-frequency motions and is well posed if suitable boundary conditions are chosen for the initial-boundary value problem. The composite grid method, a robust and efficient grid-generation technique that smoothly and accurately treats general irregular boundaries, is used to approximate the two-dimensional version of the reduced system for oceanography on irregular ocean basins. A change-of-variable technique that substantially increases the accuracy of the model and a method for efficiently solving the elliptic equation for the geopotential are discussed. Numerical resultsmore » are presented for circular and kidney-shaped basins by using a set of analytic solutions constructed in this paper.« less

Matter-wave two-dimensional solitons in crossed linear and nonlinear optical lattices

NASA Astrophysics Data System (ADS)

da Luz, H. L. F.; Abdullaev, F. Kh.; Gammal, A.; Salerno, M.; Tomio, Lauro

2010-10-01

The existence of multidimensional matter-wave solitons in a crossed optical lattice (OL) with a linear optical lattice (LOL) in the x direction and a nonlinear optical lattice (NOL) in the y direction, where the NOL can be generated by a periodic spatial modulation of the scattering length using an optically induced Feshbach resonance is demonstrated. In particular, we show that such crossed LOLs and NOLs allow for stabilizing two-dimensional solitons against decay or collapse for both attractive and repulsive interactions. The solutions for the soliton stability are investigated analytically, by using a multi-Gaussian variational approach, with the Vakhitov-Kolokolov necessary criterion for stability; and numerically, by using the relaxation method and direct numerical time integrations of the Gross-Pitaevskii equation. Very good agreement of the results corresponding to both treatments is observed.

Numerical solution of the Navier-Stokes equations about three-dimensional configurations: A survey

NASA Technical Reports Server (NTRS)

Holst, Terry L.

1987-01-01

The numerical solution of the Navier-Stokes equations about three-dimensional configurations is reviewed. Formulational and computational requirements for the various Navier-Stokes approaches are examined for typical problems including the viscous flow field solution about a complete aerospace vehicle. Recent computed results, with experimental comparisons when available, are presented to highlight the presentation. The future of Navier-Stokes applications in three-dimensions is seen to be rapidly expanding across a broad front including internal and external flows, and flows across the entire speed regime from incompressible to hypersonic applications. Prospects for the future are described and recommendations for areas of concentrated research are indicated.

New numerical solutions of three-dimensional compressible hydrodynamic convection. [in stars

NASA Technical Reports Server (NTRS)

Hossain, Murshed; Mullan, D. J.

1990-01-01

Numerical solutions of three-dimensional compressible hydrodynamics (including sound waves) in a stratified medium with open boundaries are presented. Convergent/divergent points play a controlling role in the flows, which are dominated by a single frequency related to the mean sound crossing time. Superposed on these rapid compressive flows, slower eddy-like flows eventually create convective transport. The solutions contain small structures stacked on top of larger ones, with vertical scales equal to the local pressure scale heights, H sub p. Although convective transport starts later in the evolution, vertical scales of H sub p are apparently selected at much earlier times by nonlinear compressive effects.

Parallel 3-D numerical simulation of dielectric barrier discharge plasma actuators

NASA Astrophysics Data System (ADS)

Houba, Tomas

Dielectric barrier discharge plasma actuators have shown promise in a range of applications including flow control, sterilization and ozone generation. Developing numerical models of plasma actuators is of great importance, because a high-fidelity parallel numerical model allows new design configurations to be tested rapidly. Additionally, it provides a better understanding of the plasma actuator physics which is useful for further innovation. The physics of plasma actuators is studied numerically. A loosely coupled approach is utilized for the coupling of the plasma to the neutral fluid. The state of the art in numerical plasma modeling is advanced by the development of a parallel, three-dimensional, first-principles model with detailed air chemistry. The model incorporates 7 charged species and 18 reactions, along with a solution of the electron energy equation. To the author's knowledge, a parallel three-dimensional model of a gas discharge with a detailed air chemistry model and the solution of electron energy is unique. Three representative geometries are studied using the gas discharge model. The discharge of gas between two parallel electrodes is used to validate the air chemistry model developed for the gas discharge code. The gas discharge model is then applied to the discharge produced by placing a dc powered wire and grounded plate electrodes in a channel. Finally, a three-dimensional simulation of gas discharge produced by electrodes placed inside a riblet is carried out. The body force calculated with the gas discharge model is loosely coupled with a fluid model to predict the induced flow inside the riblet.

Separated flow over bodies of revolution using an unsteady discrete-vorticity cross wake. Part 2: Computer program description

NASA Technical Reports Server (NTRS)

Marshall, F. J.; Deffenbaugh, F. D.

1974-01-01

A method is developed to determine the flow field of a body of revolution in separated flow. The computer was used to integrate various solutions and solution properties of the sub-flow fields which made up the entire flow field without resorting to a finite difference solution to the complete Navier-Stokes equations. The technique entails the use of the unsteady cross flow analogy and a new solution to the two-dimensional unsteady separated flow problem based upon an unsteady, discrete-vorticity wake. Data for the forces and moments on aerodynamic bodies at low speeds and high angle of attack (outside the range of linear inviscid theories) such that the flow is substantially separated are produced which compare well with experimental data. In addition, three dimensional steady separated regions and wake vortex patterns are determined. The computer program developed to perform the numerical calculations is described.

Transient three-dimensional thermal-hydraulic analysis of nuclear reactor fuel rod arrays: general equations and numerical scheme

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

Wnek, W.J.; Ramshaw, J.D.; Trapp, J.A.

1975-11-01

A mathematical model and a numerical solution scheme for thermal- hydraulic analysis of fuel rod arrays are given. The model alleviates the two major deficiencies associated with existing rod array analysis models, that of a correct transverse momentum equation and the capability of handling reversing and circulatory flows. Possible applications of the model include steady state and transient subchannel calculations as well as analysis of flows in heat exchangers, other engineering equipment, and porous media. (auth)

Numerical realization of the variational method for generating self-trapped beams

NASA Astrophysics Data System (ADS)

Duque, Erick I.; Lopez-Aguayo, Servando; Malomed, Boris A.

2018-03-01

We introduce a numerical variational method based on the Rayleigh-Ritz optimization principle for predicting two-dimensional self-trapped beams in nonlinear media. This technique overcomes the limitation of the traditional variational approximation in performing analytical Lagrangian integration and differentiation. Approximate soliton solutions of a generalized nonlinear Schr\\"odinger equation are obtained, demonstrating robustness of the beams of various types (fundamental, vortices, multipoles, azimuthons) in the course of their propagation. The algorithm offers possibilities to produce more sophisticated soliton profiles in general nonlinear models.

Numerical model of two-dimensional heterogeneous combustion in porous media under natural convection or forced filtration

NASA Astrophysics Data System (ADS)

Lutsenko, Nickolay A.

2018-03-01

A novel mathematical model and original numerical method for investigating the two-dimensional waves of heterogeneous combustion in porous media are proposed and described in detail. The mathematical model is constructed within the framework of the model of interacting interpenetrating continua and includes equations of state, continuity, momentum conservation and energy for solid and gas phases. Combustion, considered in the paper, is due to the exothermic reaction between fuel in the porous solid medium and oxidiser contained in the gas flowing through the porous object. The original numerical method is based on a combination of explicit and implicit finite-difference schemes. A distinctive feature of the proposed model is that the gas velocity at the open boundaries (inlet and outlet) of the porous object is unknown and has to be found from the solution of the problem, i.e. the flow rate of the gas regulates itself. This approach allows processes to be modelled not only under forced filtration, but also under free convection, when there is no forced gas input in porous objects, which is typical for many natural or anthropogenic disasters (burning of peatlands, coal dumps, landfills, grain elevators). Some two-dimensional time-dependent problems of heterogeneous combustion in porous objects have been solved using the proposed numerical method. It is shown that two-dimensional waves of heterogeneous combustion in porous media can propagate in two modes with different characteristics, as in the case of one-dimensional combustion, but the combustion front can move in a complex manner, and gas dynamics within the porous objects can be complicated. When natural convection takes place, self-sustaining combustion waves can go through the all parts of the object regardless of where an ignition zone was located, so the all combustible material in each part of the object is burned out, in contrast to forced filtration.

Problems in nonlinear acoustics: Scattering of sound by sound, parametric receiving arrays, nonlinear effects in asymmetric sound beams and pulsed finite amplitude sound beams

NASA Astrophysics Data System (ADS)

Hamilton, Mark F.

1989-08-01

Four projects are discussed in this annual summary report, all of which involve basic research in nonlinear acoustics: Scattering of Sound by Sound, a theoretical study of two nonconlinear Gaussian beams which interact to produce sum and difference frequency sound; Parametric Receiving Arrays, a theoretical study of parametric reception in a reverberant environment; Nonlinear Effects in Asymmetric Sound Beams, a numerical study of two dimensional finite amplitude sound fields; and Pulsed Finite Amplitude Sound Beams, a numerical time domain solution of the KZK equation.

A Comparison of Numerical and Analytical Radiative-Transfer Solutions for Plane Albedo in Natural Waters

EPA Science Inventory

Several numerical and analytical solutions of the radiative transfer equation (RTE) for plane albedo were compared for solar light reflection by sea water. The study incorporated the simplest case, that being a semi-infinite one-dimensional plane-parallel absorbing and scattering...

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.

Numerical studies of interacting vortices

NASA Technical Reports Server (NTRS)

Liu, G. C.; Hsu, C. H.

1985-01-01

To get a basic understanding of the physics of flowfields modeled by vortex filaments with finite vortical cores, systematic numerical studies of the interactions of two dimensional vortices and pairs of coaxial axisymmetric circular vortex rings were made. Finite difference solutions of the unsteady incompressible Navier-Stokes equations were carried out using vorticity and stream function as primary variables. Special emphasis was placed on the formulation of appropriate boundary conditions necessary for the calculations in a finite computational domain. Numerical results illustrate the interaction of vortex filaments, demonstrate when and how they merge with each other, and establish the region of validity for an asymptotic analysis.

A Comparison of Simplified Two-dimensional Flow Models Exemplified by Water Flow in a Cavern

NASA Astrophysics Data System (ADS)

Prybytak, Dzmitry; Zima, Piotr

2017-12-01

The paper shows the results of a comparison of simplified models describing a two-dimensional water flow in the example of a water flow through a straight channel sector with a cavern. The following models were tested: the two-dimensional potential flow model, the Stokes model and the Navier-Stokes model. In order to solve the first two, the boundary element method was employed, whereas to solve the Navier-Stokes equations, the open-source code library OpenFOAM was applied. The results of numerical solutions were compared with the results of measurements carried out on a test stand in a hydraulic laboratory. The measurements were taken with an ADV probe (Acoustic Doppler Velocimeter). Finally, differences between the results obtained from the mathematical models and the results of laboratory measurements were analysed.

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.

Streamline integration as a method for two-dimensional elliptic grid generation

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

Wiesenberger, M., E-mail: Matthias.Wiesenberger@uibk.ac.at; Held, M.; Einkemmer, L.

We propose a new numerical algorithm to construct a structured numerical elliptic grid of a doubly connected domain. Our method is applicable to domains with boundaries defined by two contour lines of a two-dimensional function. Furthermore, we can adapt any analytically given boundary aligned structured grid, which specifically includes polar and Cartesian grids. The resulting coordinate lines are orthogonal to the boundary. Grid points as well as the elements of the Jacobian matrix can be computed efficiently and up to machine precision. In the simplest case we construct conformal grids, yet with the help of weight functions and monitor metricsmore » we can control the distribution of cells across the domain. Our algorithm is parallelizable and easy to implement with elementary numerical methods. We assess the quality of grids by considering both the distribution of cell sizes and the accuracy of the solution to elliptic problems. Among the tested grids these key properties are best fulfilled by the grid constructed with the monitor metric approach. - Graphical abstract: - Highlights: • Construct structured, elliptic numerical grids with elementary numerical methods. • Align coordinate lines with or make them orthogonal to the domain boundary. • Compute grid points and metric elements up to machine precision. • Control cell distribution by adaption functions or monitor metrics.« less

Long-Time Numerical Integration of the Three-Dimensional Wave Equation in the Vicinity of a Moving Source

NASA Technical Reports Server (NTRS)

Ryabenkii, V. S.; Turchaninov, V. I.; Tsynkov, S. V.

1999-01-01

We propose a family of algorithms for solving numerically a Cauchy problem for the three-dimensional wave equation. The sources that drive the equation (i.e., the right-hand side) are compactly supported in space for any given time; they, however, may actually move in space with a subsonic speed. The solution is calculated inside a finite domain (e.g., sphere) that also moves with a subsonic speed and always contains the support of the right-hand side. The algorithms employ a standard consistent and stable explicit finite-difference scheme for the wave equation. They allow one to calculate tile solution for arbitrarily long time intervals without error accumulation and with the fixed non-growing amount of tile CPU time and memory required for advancing one time step. The algorithms are inherently three-dimensional; they rely on the presence of lacunae in the solutions of the wave equation in oddly dimensional spaces. The methodology presented in the paper is, in fact, a building block for constructing the nonlocal highly accurate unsteady artificial boundary conditions to be used for the numerical simulation of waves propagating with finite speed over unbounded domains.

A multiple-block multigrid method for the solution of the three-dimensional Euler and Navier-Stokes equations

NASA Technical Reports Server (NTRS)

Atkins, Harold

1991-01-01

A multiple block multigrid method for the solution of the three dimensional Euler and Navier-Stokes equations is presented. The basic flow solver is a cell vertex method which employs central difference spatial approximations and Runge-Kutta time stepping. The use of local time stepping, implicit residual smoothing, multigrid techniques and variable coefficient numerical dissipation results in an efficient and robust scheme is discussed. The multiblock strategy places the block loop within the Runge-Kutta Loop such that accuracy and convergence are not affected by block boundaries. This has been verified by comparing the results of one and two block calculations in which the two block grid is generated by splitting the one block grid. Results are presented for both Euler and Navier-Stokes computations of wing/fuselage combinations.

Asymptotic and spectral analysis of the gyrokinetic-waterbag integro-differential operator in toroidal geometry

NASA Astrophysics Data System (ADS)

Besse, Nicolas; Coulette, David

2016-08-01

Achieving plasmas with good stability and confinement properties is a key research goal for magnetic fusion devices. The underlying equations are the Vlasov-Poisson and Vlasov-Maxwell (VPM) equations in three space variables, three velocity variables, and one time variable. Even in those somewhat academic cases where global equilibrium solutions are known, studying their stability requires the analysis of the spectral properties of the linearized operator, a daunting task. We have identified a model, for which not only equilibrium solutions can be constructed, but many of their stability properties are amenable to rigorous analysis. It uses a class of solution to the VPM equations (or to their gyrokinetic approximations) known as waterbag solutions which, in particular, are piecewise constant in phase-space. It also uses, not only the gyrokinetic approximation of fast cyclotronic motion around magnetic field lines, but also an asymptotic approximation regarding the magnetic-field-induced anisotropy: the spatial variation along the field lines is taken much slower than across them. Together, these assumptions result in a drastic reduction in the dimensionality of the linearized problem, which becomes a set of two nested one-dimensional problems: an integral equation in the poloidal variable, followed by a one-dimensional complex Schrödinger equation in the radial variable. We show here that the operator associated to the poloidal variable is meromorphic in the eigenparameter, the pulsation frequency. We also prove that, for all but a countable set of real pulsation frequencies, the operator is compact and thus behaves mostly as a finite-dimensional one. The numerical algorithms based on such ideas have been implemented in a companion paper [D. Coulette and N. Besse, "Numerical resolution of the global eigenvalue problem for gyrokinetic-waterbag model in toroidal geometry" (submitted)] and were found to be surprisingly close to those for the original gyrokinetic-Vlasov equations. The purpose of the present paper is to make these new ideas accessible to two readerships: applied mathematicians and plasma physicists.

A comparison of solute-transport solution techniques based on inverse modelling results

USGS Publications Warehouse

Mehl, S.; Hill, M.C.

2000-01-01

Five common numerical techniques (finite difference, predictor-corrector, total-variation-diminishing, method-of-characteristics, and modified-method-of-characteristics) were tested using simulations of a controlled conservative tracer-test experiment through a heterogeneous, two-dimensional sand tank. The experimental facility was constructed using randomly distributed homogeneous blocks of five sand types. This experimental model provides an outstanding opportunity to compare the solution techniques because of the heterogeneous hydraulic conductivity distribution of known structure, and the availability of detailed measurements with which to compare simulated concentrations. The present work uses this opportunity to investigate how three common types of results-simulated breakthrough curves, sensitivity analysis, and calibrated parameter values-change in this heterogeneous situation, given the different methods of simulating solute transport. The results show that simulated peak concentrations, even at very fine grid spacings, varied because of different amounts of numerical dispersion. Sensitivity analysis results were robust in that they were independent of the solution technique. They revealed extreme correlation between hydraulic conductivity and porosity, and that the breakthrough curve data did not provide enough information about the dispersivities to estimate individual values for the five sands. However, estimated hydraulic conductivity values are significantly influenced by both the large possible variations in model dispersion and the amount of numerical dispersion present in the solution technique.Five common numerical techniques (finite difference, predictor-corrector, total-variation-diminishing, method-of-characteristics, and modified-method-of-characteristics) were tested using simulations of a controlled conservative tracer-test experiment through a heterogeneous, two-dimensional sand tank. The experimental facility was constructed using randomly distributed homogeneous blocks of five sand types. This experimental model provides an outstanding opportunity to compare the solution techniques because of the heterogeneous hydraulic conductivity distribution of known structure, and the availability of detailed measurements with which to compare simulated concentrations. The present work uses this opportunity to investigate how three common types of results - simulated breakthrough curves, sensitivity analysis, and calibrated parameter values - change in this heterogeneous situation, given the different methods of simulating solute transport. The results show that simulated peak concentrations, even at very fine grid spacings, varied because of different amounts of numerical dispersion. Sensitivity analysis results were robust in that they were independent of the solution technique. They revealed extreme correlation between hydraulic conductivity and porosity, and that the breakthrough curve data did not provide enough information about the dispersivities to estimate individual values for the five sands. However, estimated hydraulic conductivity values are significantly influenced by both the large possible variations in model dispersion and the amount of numerical dispersion present in the solution technique.

An exact solution for the solidification of a liquid slab of binary mixture

NASA Technical Reports Server (NTRS)

Antar, B. N.; Collins, F. G.; Aumalia, A. E.

1986-01-01

The time dependent temperature and concentration profiles of a one dimensional finite slab of a binary liquid alloy is investigated during solidification. The governing equations are reduced to a set of coupled, nonlinear initial value problems using the method outlined by Meyer. Two methods will be used to solve these equations. The first method uses a Runge-Kutta-Fehlberg integrator to solve the equations numerically. The second method comprises of finding closed form solutions of the equations.

Discontinuous finite element method for vector radiative transfer

NASA Astrophysics Data System (ADS)

Wang, Cun-Hai; Yi, Hong-Liang; Tan, He-Ping

2017-03-01

The discontinuous finite element method (DFEM) is applied to solve the vector radiative transfer in participating media. The derivation in a discrete form of the vector radiation governing equations is presented, in which the angular space is discretized by the discrete-ordinates approach with a local refined modification, and the spatial domain is discretized into finite non-overlapped discontinuous elements. The elements in the whole solution domain are connected by modelling the boundary numerical flux between adjacent elements, which makes the DFEM numerically stable for solving radiative transfer equations. Several various problems of vector radiative transfer are tested to verify the performance of the developed DFEM, including vector radiative transfer in a one-dimensional parallel slab containing a Mie/Rayleigh/strong forward scattering medium and a two-dimensional square medium. The fact that DFEM results agree very well with the benchmark solutions in published references shows that the developed DFEM in this paper is accurate and effective for solving vector radiative transfer problems.

Diffraction of a plane wave by a three-dimensional corner

NASA Technical Reports Server (NTRS)

Ting, L.; Kung, F.

1971-01-01

By the superposition of the conical solution for the diffraction of a plane pulse by a three dimensional corner, the solution for a general incident plane wave is constructed. A numerical program is presented for the computation of the pressure distribution on the surface due to an incident plane wave of any wave form and at any incident angle. Numerical examples are presented to show the pressure signature at several points on the surface due to incident wave with a front shock wave, two shock waves in succession, or a compression wave with same peak pressure. The examples show that when the distance of a point on the surface from the edges or the vertex is comparable to the distance for the front pressure raise to reach the maximum, the peak pressure at that point can be much less than that given by a regular reflection, because the diffracted wave front arrives at that point prior to the arrival of the peak incident wave.

The three-dimensional turbulent boundary layer near a plane of symmetry

NASA Technical Reports Server (NTRS)

Degani, A. T.; Smith, F. T.; Walker, J. D. A.

1992-01-01

The asymptotic structure of the three-dimensional turbulent boundary layer near a plane of symmetry is considered in the limit of large Reynolds number. A self-consistent two-layer structure is shown to exist wherein the streamwise velocity is brought to rest through an outer defect layer and an inner wall layer in a manner similar to that in two-dimensional boundary layers. The cross-stream velocity distribution is more complex and two terms in the asymptotic expansion are required to yield a complete profile which is shown to exhibit a logarithmic region. The flow in the inner wall layer is demonstrated to be collateral to leading order; pressure-gradient effects are formally of higher order but can cause the velocity profile to skew substantially near the wall at the large but finite Reynolds numbers encountered in practice. The governing set of ordinary differential equations describing a self-similar flow is derived. The calculated numerical solutions of these equations are matched asymptotically to an inner wall-layer solution and the results show trends that are consistent with experimental observations.

Crack growth in bonded elastic half planes

NASA Technical Reports Server (NTRS)

Goree, J. G.

1975-01-01

Two solutions were developed for the two dimensional problem of bonded linearly elastic half-planes. For each solution, numerical results are presented for the stress intensity factors, strain energy release rate, stresses, and displacements. The behavior predicted by the studies was investigated experimentally using polymers for the material pairs. Close agreement was found for the critical stress intensity factor at fracture for the perpendicular crack near the interface. Fracture along the interface proved to be inconclusive due to difficulties in obtaining a brittle bond. Some interesting and predictable behavior regarding the potential for the crack to cross the interface was observed and is discussed.

Reflection and interference of electromagnetic waves in inhomogeneous media

NASA Technical Reports Server (NTRS)

Geiger, F. E.; Kyle, H. L.

1973-01-01

Solutions were obtained of the wave equation for a plane horizontally polarized electro-magnetic wave incident on a semi infinite two dimensional inhomogeneous medium. Two problems were considered: An inhomogeneous half space, and an inhomogeneous layer of arbitrary thickness. Solutions of the wave equation were obtained in terms of Hankel functions with complex arguments. Numerical calculations were made of the reflection coefficient R at the interface of the homogeneous medium. The startling results show that the reflection coefficient for a complex dielectric constant with gradient, can be less than that of the same medium with zero gradient.

Grid-converged solution and analysis of the unsteady viscous flow in a two-dimensional shock tube

NASA Astrophysics Data System (ADS)

Zhou, Guangzhao; Xu, Kun; Liu, Feng

2018-01-01

The flow in a shock tube is extremely complex with dynamic multi-scale structures of sharp fronts, flow separation, and vortices due to the interaction of the shock wave, the contact surface, and the boundary layer over the side wall of the tube. Prediction and understanding of the complex fluid dynamics are of theoretical and practical importance. It is also an extremely challenging problem for numerical simulation, especially at relatively high Reynolds numbers. Daru and Tenaud ["Evaluation of TVD high resolution schemes for unsteady viscous shocked flows," Comput. Fluids 30, 89-113 (2001)] proposed a two-dimensional model problem as a numerical test case for high-resolution schemes to simulate the flow field in a square closed shock tube. Though many researchers attempted this problem using a variety of computational methods, there is not yet an agreed-upon grid-converged solution of the problem at the Reynolds number of 1000. This paper presents a rigorous grid-convergence study and the resulting grid-converged solutions for this problem by using a newly developed, efficient, and high-order gas-kinetic scheme. Critical data extracted from the converged solutions are documented as benchmark data. The complex fluid dynamics of the flow at Re = 1000 are discussed and analyzed in detail. Major phenomena revealed by the numerical computations include the downward concentration of the fluid through the curved shock, the formation of the vortices, the mechanism of the shock wave bifurcation, the structure of the jet along the bottom wall, and the Kelvin-Helmholtz instability near the contact surface. Presentation and analysis of those flow processes provide important physical insight into the complex flow physics occurring in a shock tube.

High-order scheme for the source-sink term in a one-dimensional water temperature model

PubMed Central

Jing, Zheng; Kang, Ling

2017-01-01

The source-sink term in water temperature models represents the net heat absorbed or released by a water system. This term is very important because it accounts for solar radiation that can significantly affect water temperature, especially in lakes. However, existing numerical methods for discretizing the source-sink term are very simplistic, causing significant deviations between simulation results and measured data. To address this problem, we present a numerical method specific to the source-sink term. A vertical one-dimensional heat conduction equation was chosen to describe water temperature changes. A two-step operator-splitting method was adopted as the numerical solution. In the first step, using the undetermined coefficient method, a high-order scheme was adopted for discretizing the source-sink term. In the second step, the diffusion term was discretized using the Crank-Nicolson scheme. The effectiveness and capability of the numerical method was assessed by performing numerical tests. Then, the proposed numerical method was applied to a simulation of Guozheng Lake (located in central China). The modeling results were in an excellent agreement with measured data. PMID:28264005

High-order scheme for the source-sink term in a one-dimensional water temperature model.

PubMed

Jing, Zheng; Kang, Ling

2017-01-01

The source-sink term in water temperature models represents the net heat absorbed or released by a water system. This term is very important because it accounts for solar radiation that can significantly affect water temperature, especially in lakes. However, existing numerical methods for discretizing the source-sink term are very simplistic, causing significant deviations between simulation results and measured data. To address this problem, we present a numerical method specific to the source-sink term. A vertical one-dimensional heat conduction equation was chosen to describe water temperature changes. A two-step operator-splitting method was adopted as the numerical solution. In the first step, using the undetermined coefficient method, a high-order scheme was adopted for discretizing the source-sink term. In the second step, the diffusion term was discretized using the Crank-Nicolson scheme. The effectiveness and capability of the numerical method was assessed by performing numerical tests. Then, the proposed numerical method was applied to a simulation of Guozheng Lake (located in central China). The modeling results were in an excellent agreement with measured data.

NUMERICAL ANALYSES FOR TREATING DIFFUSION IN SINGLE-, TWO-, AND THREE-PHASE BINARY ALLOY SYSTEMS

NASA Technical Reports Server (NTRS)

Tenney, D. R.

1994-01-01

This package consists of a series of three computer programs for treating one-dimensional transient diffusion problems in single and multiple phase binary alloy systems. An accurate understanding of the diffusion process is important in the development and production of binary alloys. Previous solutions of the diffusion equations were highly restricted in their scope and application. The finite-difference solutions developed for this package are applicable for planar, cylindrical, and spherical geometries with any diffusion-zone size and any continuous variation of the diffusion coefficient with concentration. Special techniques were included to account for differences in modal volumes, initiation and growth of an intermediate phase, disappearance of a phase, and the presence of an initial composition profile in the specimen. In each analysis, an effort was made to achieve good accuracy while minimizing computation time. The solutions to the diffusion equations for single-, two-, and threephase binary alloy systems are numerically calculated by the three programs NAD1, NAD2, and NAD3. NAD1 treats the diffusion between pure metals which belong to a single-phase system. Diffusion in this system is described by a one-dimensional Fick's second law and will result in a continuous composition variation. For computational purposes, Fick's second law is expressed as an explicit second-order finite difference equation. Finite difference calculations are made by choosing the grid spacing small enough to give convergent solutions of acceptable accuracy. NAD2 treats diffusion between pure metals which form a two-phase system. Diffusion in the twophase system is described by two partial differential equations (a Fick's second law for each phase) and an interface-flux-balance equation which describes the location of the interface. Actual interface motion is obtained by a mass conservation procedure. To account for changes in the thicknesses of the two phases as diffusion progresses, a variable grid technique developed by Murray and Landis is employed. These equations are expressed in finite difference form and solved numerically. Program NAD3 treats diffusion between pure metals which form a two-phase system with an intermediate third phase. Diffusion in the three-phase system is described by three partial differential expressions of Fick's second law and two interface-flux-balance equations. As with the two-phase case, a variable grid finite difference is used to numerically solve the diffusion equations. Computation time is minimized without sacrificing solution accuracy by treating the three-phase problem as a two-phase problem when the thickness of the intermediate phase is less than a preset value. Comparisons between these programs and other solutions have shown excellent agreement. The programs are written in FORTRAN IV for batch execution on the CDC 6600 with a central memory requirement of approximately 51K (octal) 60 bit words.

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

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.

Impedance Eduction in Sound Fields With Peripherally Varying Liners and Flow

NASA Technical Reports Server (NTRS)

Watson, W. R.; Jones, M. G.

2015-01-01

A two-dimensional impedance eduction theory is extended to three-dimensional sound fields and peripherally varying duct liners. The approach is to first measure the acoustic pressure field at a series of flush-mounted wall microphones located around the periphery of the flow duct. The numerical solution for the acoustic pressure field at these microphones is also obtained by solving the three-dimensional convected Helmholtz equation using the finite element method. A quadratic objective function based on the difference between the measured and finite element solution is constructed and the unknown impedance function is obtained by minimizing this objective function. Impedance spectra educed for two uniform-structure liners (a wire-mesh and a conventional liner) and a hard-soft-hard peripherally varying liner (for which the soft segment is that of the conventional liner) are presented. Results are presented at three mean flow Mach numbers and fourteen sound source frequencies. The impedance spectra of the uniform-structure liners are also computed using a two-dimensional impedance eduction theory. The primary conclusions of the study are: 1) when measured data is used with the uniform-structure liners, the three-dimensional theory reproduces the same impedance spectra as the two-dimensional theory except for frequencies corresponding to very low or very high liner attenuation; and 2) good agreement between the educed impedance spectra of the uniform structure conventional liner and the soft segment of the peripherally varying liner is obtained.

Analysis of sound propagation in ducts using the wave envelope concept

NASA Technical Reports Server (NTRS)

Baumeister, K. J.

1974-01-01

A finite difference formulation is presented for sound propagation in a rectangular two-dimensional duct without steady flow for plane wave input. Before the difference equations are formulated, the governing Helmholtz equation is first transformed to a form whose solution does not oscillate along the length of the duct. This transformation reduces the required number of grid points by an order of magnitude, and the number of grid points becomes independent of the sound frequency. Physically, the transformed pressure represents the amplitude of the conventional sound wave. Example solutions are presented for sound propagation in a one-dimensional straight hard-wall duct and in a two-dimensional straight soft-wall duct without steady flow. The numerical solutions show evidence of the existence along the duct wall of a developing acoustic pressure diffusion boundary layer which is similar in nature to the conventional viscous flow boundary layer. In order to better illustrate this concept, the wave equation and boundary conditions are written such that the frequency no longer appears explicitly in them. The frequency effects in duct propagation can be visualized solely as an expansion and stretching of the suppressor duct.

Iterative solution of the inverse Cauchy problem for an elliptic equation by the conjugate gradient method

NASA Astrophysics Data System (ADS)

Vasil'ev, V. I.; Kardashevsky, A. M.; Popov, V. V.; Prokopev, G. A.

2017-10-01

This article presents results of computational experiment carried out using a finite-difference method for solving the inverse Cauchy problem for a two-dimensional elliptic equation. The computational algorithm involves an iterative determination of the missing boundary condition from the override condition using the conjugate gradient method. The results of calculations are carried out on the examples with exact solutions as well as at specifying an additional condition with random errors are presented. Results showed a high efficiency of the iterative method of conjugate gradients for numerical solution

A general algorithm using finite element method for aerodynamic configurations at low speeds

NASA Technical Reports Server (NTRS)

Balasubramanian, R.

1975-01-01

A finite element algorithm for numerical simulation of two-dimensional, incompressible, viscous flows was developed. The Navier-Stokes equations are suitably modelled to facilitate direct solution for the essential flow parameters. A leap-frog time differencing and Galerkin minimization of these model equations yields the finite element algorithm. The finite elements are triangular with bicubic shape functions approximating the solution space. The finite element matrices are unsymmetrically banded to facilitate savings in storage. An unsymmetric L-U decomposition is performed on the finite element matrices to obtain the solution for the boundary value problem.

Analytical solutions to time-fractional partial differential equations in a two-dimensional multilayer annulus

NASA Astrophysics Data System (ADS)

Chen, Shanzhen; Jiang, Xiaoyun

2012-08-01

In this paper, analytical solutions to time-fractional partial differential equations in a multi-layer annulus are presented. The final solutions are obtained in terms of Mittag-Leffler function by using the finite integral transform technique and Laplace transform technique. In addition, the classical diffusion equation (α=1), the Helmholtz equation (α→0) and the wave equation (α=2) are discussed as special cases. Finally, an illustrative example problem for the three-layer semi-circular annular region is solved and numerical results are presented graphically for various kind of order of fractional derivative.

A numerical solution of the Navier-Stokes equations for chemically nonequilibrium, merged stagnation shock layers on spheres and two-dimensional cylinders in air

NASA Technical Reports Server (NTRS)

Johnston, K. D.; Hendricks, W. L.

1978-01-01

Results of solving the Navier-Stokes equations for chemically nonequilibrium, merged stagnation shock layers on spheres and two-dimensional cylinders are presented. The effects of wall catalysis and slip are also examined. The thin shock layer assumption is not made, and the thick viscous shock is allowed to develop within the computational domain. The results show good comparison with existing data. Due to the more pronounced merging of shock layer and boundary layer for the sphere, the heating rates for spheres become higher than those for cylinders as the altitude is increased.

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

Zou, Ling; Zhao, Haihua; Zhang, Hongbin

Here, the one-dimensional water faucet problem is one of the classical benchmark problems originally proposed by Ransom to study the two-fluid two-phase flow model. With certain simplifications, such as massless gas phase and no wall and interfacial frictions, analytical solutions had been previously obtained for the transient liquid velocity and void fraction distribution. The water faucet problem and its analytical solutions have been widely used for the purposes of code assessment, benchmark and numerical verifications. In our previous study, the Ransom’s solutions were used for the mesh convergence study of a high-resolution spatial discretization scheme. It was found that, atmore » the steady state, an anticipated second-order spatial accuracy could not be achieved, when compared to the existing Ransom’s analytical solutions. A further investigation showed that the existing analytical solutions do not actually satisfy the commonly used two-fluid single-pressure two-phase flow equations. In this work, we present a new set of analytical solutions of the water faucet problem at the steady state, considering the gas phase density’s effect on pressure distribution. This new set of analytical solutions are used for mesh convergence studies, from which anticipated second-order of accuracy is achieved for the 2nd order spatial discretization scheme. In addition, extended Ransom’s transient solutions for the gas phase velocity and pressure are derived, with the assumption of decoupled liquid and gas pressures. Numerical verifications on the extended Ransom’s solutions are also presented.« less

A Conserving Discretization for the Free Boundary in a Two-Dimensional Stefan Problem

NASA Astrophysics Data System (ADS)

Segal, Guus; Vuik, Kees; Vermolen, Fred

1998-03-01

The dissolution of a disk-likeAl2Cuparticle is considered. A characteristic property is that initially the particle has a nonsmooth boundary. The mathematical model of this dissolution process contains a description of the particle interface, of which the position varies in time. Such a model is called a Stefan problem. It is impossible to obtain an analytical solution for a general two-dimensional Stefan problem, so we use the finite element method to solve this problem numerically. First, we apply a classical moving mesh method. Computations show that after some time steps the predicted particle interface becomes very unrealistic. Therefore, we derive a new method for the displacement of the free boundary based on the balance of atoms. This method leads to good results, also, for nonsmooth boundaries. Some numerical experiments are given for the dissolution of anAl2Cuparticle in anAl-Cualloy.

A numerical formulation and algorithm for limit and shakedown analysis of large-scale elastoplastic structures

NASA Astrophysics Data System (ADS)

Peng, Heng; Liu, Yinghua; Chen, Haofeng

2018-05-01

In this paper, a novel direct method called the stress compensation method (SCM) is proposed for limit and shakedown analysis of large-scale elastoplastic structures. Without needing to solve the specific mathematical programming problem, the SCM is a two-level iterative procedure based on a sequence of linear elastic finite element solutions where the global stiffness matrix is decomposed only once. In the inner loop, the static admissible residual stress field for shakedown analysis is constructed. In the outer loop, a series of decreasing load multipliers are updated to approach to the shakedown limit multiplier by using an efficient and robust iteration control technique, where the static shakedown theorem is adopted. Three numerical examples up to about 140,000 finite element nodes confirm the applicability and efficiency of this method for two-dimensional and three-dimensional elastoplastic structures, with detailed discussions on the convergence and the accuracy of the proposed algorithm.

ISCFD Nagoya 1989 - International Symposium on Computational Fluid Dynamics, 3rd, Nagoya, Japan, Aug. 28-31, 1989, Technical Papers

NASA Astrophysics Data System (ADS)

Recent advances in computational fluid dynamics are discussed in reviews and reports. Topics addressed include large-scale LESs for turbulent pipe and channel flows, numerical solutions of the Euler and Navier-Stokes equations on parallel computers, multigrid methods for steady high-Reynolds-number flow past sudden expansions, finite-volume methods on unstructured grids, supersonic wake flow on a blunt body, a grid-characteristic method for multidimensional gas dynamics, and CIC numerical simulation of a wave boundary layer. Consideration is given to vortex simulations of confined two-dimensional jets, supersonic viscous shear layers, spectral methods for compressible flows, shock-wave refraction at air/water interfaces, oscillatory flow in a two-dimensional collapsible channel, the growth of randomness in a spatially developing wake, and an efficient simplex algorithm for the finite-difference and dynamic linear-programming method in optimal potential control.

Three-dimensional low Reynolds number flows with a free surface

NASA Technical Reports Server (NTRS)

Degani, D.; Gutfinger, C.

1977-01-01

The two-dimensional leveling problem (Degani, Gutfinger, 1976) is extended to three dimensions in the case where the flow Re number is very low and attention is paid to the free surface boundary condition with surface tension effects included. The no-slip boundary condition on the wall is observed. The numerical solution falls back on the Marker and Cell (MAC) method (Harlow and Welch, 1965) with the computation region divided into a finite number of stationary rectangular cells (or boxes in the 3-D case) and fluid flow traverses the cells (or boxes).

The method of space-time conservation element and solution element-applications to one-dimensional and two-dimensional time-marching flow problems

NASA Technical Reports Server (NTRS)

Chang, Sin-Chung; Wang, Xiao-Yen; Chow, Chuen-Yen

1995-01-01

A nontraditional numerical method for solving conservation laws is being developed. The new method is designed from a physicist's perspective, i.e., its development is based more on physics than numerics. Even though it uses only the simplest approximation techniques, a 2D time-marching Euler solver developed recently using the new method is capable of generating nearly perfect solutions for a 2D shock reflection problem used by Helen Yee and others. Moreover, a recent application of this solver to computational aeroacoustics (CAA) problems reveals that: (1) accuracy of its results is comparable to that of a 6th order compact difference scheme even though nominally the current solver is only of 2nd-order accuracy; (2) generally, the non-reflecting boundary condition can be implemented in a simple way without involving characteristic variables; and (3) most importantly, the current solver is capable of handling both continuous and discontinuous flows very well and thus provides a unique numerical tool for solving those flow problems where the interactions between sound waves and shocks are important, such as the noise field around a supersonic over- or under-expansion jet.

Steering particles by breaking symmetries

NASA Astrophysics Data System (ADS)

Bet, Bram; Samin, Sela; Georgiev, Rumen; Burak Eral, Huseyin; van Roij, René

2018-06-01

We derive general equations of motions for highly-confined particles that perform quasi-two-dimensional motion in Hele-Shaw channels, which we solve analytically, aiming to derive design principles for self-steering particles. Based on symmetry properties of a particle, its equations of motion can be simplified, where we retrieve an earlier-known equation of motion for the orientation of dimer particles consisting of disks (Uspal et al 2013 Nat. Commun. 4), but now in full generality. Subsequently, these solutions are compared with particle trajectories that are obtained numerically. For mirror-symmetric particles, excellent agreement between the analytical and numerical solutions is found. For particles lacking mirror symmetry, the analytic solutions provide means to classify the motion based on particle geometry, while we find that taking the side-wall interactions into account is important to accurately describe the trajectories.

Numerical studies of unsteady two dimensional subsonic flows using the ICE method. Ph.D. Thesis - Toledo Univ.

NASA Technical Reports Server (NTRS)

Wieber, P. R.

1973-01-01

A numerical program was developed to compute transient compressible and incompressible laminar flows in two dimensions with multicomponent mixing and chemical reaction. The algorithm used the Los Alamos Scientific Laboratory ICE (Implicit Continuous-Fluid Eulerian) method as its base. The program can compute both high and low speed compressible flows. The numerical program incorporating the stabilization techniques was quite successful in treating both old and new problems. Detailed calculations of coaxial flow very close to the entry plane were possible. The program treated complex flows such as the formation and downstream growth of a recirculation cell. An implicit solution of the species equation predicted mixing and reaction rates which compared favorably with the literature.

A theoretical study of mixing downstream of transverse injection into a supersonic boundary layer

NASA Technical Reports Server (NTRS)

Baker, A. J.; Zelazny, S. W.

1972-01-01

A theoretical and analytical study was made of mixing downstream of transverse hydrogen injection, from single and multiple orifices, into a Mach 4 air boundary layer over a flat plate. Numerical solutions to the governing three-dimensional, elliptic boundary layer equations were obtained using a general purpose computer program. Founded upon a finite element solution algorithm. A prototype three-dimensional turbulent transport model was developed using mixing length theory in the wall region and the mass defect concept in the outer region. Excellent agreement between the computed flow field and experimental data for a jet/freestream dynamic pressure ratio of unity was obtained in the centerplane region of the single-jet configuration. Poorer agreement off centerplane suggests an inadequacy of the extrapolated two-dimensional turbulence model. Considerable improvement in off-centerplane computational agreement occured for a multi-jet configuration, using the same turbulent transport model.

Locating CVBEM collocation points for steady state heat transfer problems

USGS Publications Warehouse

Hromadka, T.V.

1985-01-01

The Complex Variable Boundary Element Method or CVBEM provides a highly accurate means of developing numerical solutions to steady state two-dimensional heat transfer problems. The numerical approach exactly solves the Laplace equation and satisfies the boundary conditions at specified points on the boundary by means of collocation. The accuracy of the approximation depends upon the nodal point distribution specified by the numerical analyst. In order to develop subsequent, refined approximation functions, four techniques for selecting additional collocation points are presented. The techniques are compared as to the governing theory, representation of the error of approximation on the problem boundary, the computational costs, and the ease of use by the numerical analyst. ?? 1985.

Modeling Three-Dimensional Flow in Confined Aquifers by Superposition of Both Two- and Three-Dimensional Analytic Functions

NASA Astrophysics Data System (ADS)

Haitjema, Henk M.

1985-10-01

A technique is presented to incorporate three-dimensional flow in a Dupuit-Forchheimer model. The method is based on superposition of approximate analytic solutions to both two- and three-dimensional flow features in a confined aquifer of infinite extent. Three-dimensional solutions are used in the domain of interest, while farfield conditions are represented by two-dimensional solutions. Approximate three- dimensional solutions have been derived for a partially penetrating well and a shallow creek. Each of these solutions satisfies the condition that no flow occurs across the confining layers of the aquifer. Because of this condition, the flow at some distance of a three-dimensional feature becomes nearly horizontal. Consequently, remotely from a three-dimensional feature, its three-dimensional solution is replaced by a corresponding two-dimensional one. The latter solution is trivial as compared to its three-dimensional counterpart, and its use greatly enhances the computational efficiency of the model. As an example, the flow is modeled between a partially penetrating well and a shallow creek that occur in a regional aquifer system.

A new extrapolation cascadic multigrid method for three dimensional elliptic boundary value problems

NASA Astrophysics Data System (ADS)

Pan, Kejia; He, Dongdong; Hu, Hongling; Ren, Zhengyong

2017-09-01

In this paper, we develop a new extrapolation cascadic multigrid method, which makes it possible to solve three dimensional elliptic boundary value problems with over 100 million unknowns on a desktop computer in half a minute. First, by combining Richardson extrapolation and quadratic finite element (FE) interpolation for the numerical solutions on two-level of grids (current and previous grids), we provide a quite good initial guess for the iterative solution on the next finer grid, which is a third-order approximation to the FE solution. And the resulting large linear system from the FE discretization is then solved by the Jacobi-preconditioned conjugate gradient (JCG) method with the obtained initial guess. Additionally, instead of performing a fixed number of iterations as used in existing cascadic multigrid methods, a relative residual tolerance is introduced in the JCG solver, which enables us to obtain conveniently the numerical solution with the desired accuracy. Moreover, a simple method based on the midpoint extrapolation formula is proposed to achieve higher-order accuracy on the finest grid cheaply and directly. Test results from four examples including two smooth problems with both constant and variable coefficients, an H3-regular problem as well as an anisotropic problem are reported to show that the proposed method has much better efficiency compared to the classical V-cycle and W-cycle multigrid methods. Finally, we present the reason why our method is highly efficient for solving these elliptic problems.

Extension and application of the Preissmann slot model to 2D transient mixed flows

NASA Astrophysics Data System (ADS)

Maranzoni, Andrea; Dazzi, Susanna; Aureli, Francesca; Mignosa, Paolo

2015-08-01

This paper presents an extension of the Preissmann slot concept for the modeling of highly transient two-dimensional (2D) mixed flows. The classic conservative formulation of the 2D shallow water equations for free surface flows is adapted by assuming that two fictitious vertical slots, aligned along the two Cartesian plane directions and normally intersecting, are added on the ceiling of each integration element. Accordingly, transitions between free surface and pressurized flow can be handled in a natural and straightforward way by using the same set of governing equations. The opportunity of coupling free surface and pressurized flows is actually useful not only in one-dimensional (1D) problems concerning sewer systems but also for modeling 2D flooding phenomena in which the pressurization of bridges, culverts, or other crossing hydraulic structures can be expected. Numerical simulations are performed by using a shock-capturing MUSCL-Hancock finite volume scheme combined with the FORCE (First-Order Centred) solver for the evaluation of the numerical fluxes. The validation of the mathematical model is accomplished on the basis of both exact solutions of 1D discontinuous initial value problems and reference radial solutions of idealized test cases with cylindrical symmetry. Furthermore, the capability of the model to deal with practical field-scale applications is assessed by simulating the transit of a bore under an arch bridge. Numerical results show that the proposed model is suitable for the prediction of highly transient 2D mixed flows.

An application of a two-equation model of turbulence to three-dimensional chemically reacting flows

NASA Technical Reports Server (NTRS)

Lee, J.

1994-01-01

A numerical study of three dimensional chemically reacting and non-reacting flowfields is conducted using a two-equation model of turbulence. A generalized flow solver using an implicit Lower-Upper (LU) diagonal decomposition numerical technique and finite-rate chemistry has been coupled with a low-Reynolds number two-equation model of turbulence. This flow solver is then used to study chemically reacting turbulent supersonic flows inside combustors with synergetic fuel injectors. The reacting and non-reacting turbulent combustor solutions obtained are compared with zero-equation turbulence model solutions and with available experimental data. The hydrogen-air chemistry is modeled using a nine-species/eighteen reaction model. A low-Reynolds number k-epsilon model was used to model the effect of turbulence because, in general, the low-Reynolds number k-epsilon models are easier to implement numerically and are far more general than algebraic models. However, low-Reynolds number k-epsilon models require a much finer near-wall grid resolution than high-Reynolds number models to resolve accurately the near-wall physics. This is especially true in complex flowfields, where the stiff nature of the near-wall turbulence must be resolved. Therefore, the limitations imposed by the near-wall characteristics and compressible model corrections need to be evaluated further. The gradient-diffusion hypothesis is used to model the effects of turbulence on the mass diffusion process. The influence of this low-Reynolds number turbulence model on the reacting flowfield predictions was studied parametrically.

Vortex methods for separated flows

NASA Technical Reports Server (NTRS)

Spalart, Philippe R.

1988-01-01

The numerical solution of the Euler or Navier-Stokes equations by Lagrangian vortex methods is discussed. The mathematical background is presented in an elementary fashion and includes the relationship with traditional point-vortex studies, the convergence to smooth solutions of the Euler equations, and the essential differences between two- and three-dimensional cases. The difficulties in extending the method to viscous or compressible flows are explained. The overlap with the excellent review articles available is kept to a minimum and more emphasis is placed on the area of expertise, namely two-dimensional flows around bluff bodies. When solid walls are present, complete mathematical models are not available and a more heuristic attitude must be adopted. The imposition of inviscid and viscous boundary conditions without conformal mappings or image vortices and the creation of vorticity along solid walls are examined in detail. Methods for boundary-layer treatment and the question of the Kutta condition are discussed. Practical aspects and tips helpful in creating a method that really works are explained. The topics include the robustness of the method and the assessment of accuracy, vortex-core profiles, timemarching schemes, numerical dissipation, and efficient programming. Calculations of flows past streamlined or bluff bodies are used as examples when appropriate.

Two dimensional modelling of flood flows and suspended sediment transport: the case of Brenta River

NASA Astrophysics Data System (ADS)

D'Alpaos, L.; Martini, P.; Carniello, L.

2003-04-01

The paper deals with numerical modelling of flood waves and suspended sediment in plain river basins. The two dimensional depth integrated momentum and continuity equations, modified to take into account of the bottom irregularities that strongly affect the hydrodynamic and the continuity in partially dry areas (for example, during the first stages of a plain flooding and in tidal flows), are solved with a standard Galerkin finite element method using a semi-implicit numerical scheme and considering the role both of the small channel network and the regulation dispositive on the flooding wave propagation. Transport of suspended sediment and bed evolution are coupled with the flood propagation through the convection-dispersion equation and the Exner's equation. Results of a real case study are presented in which the effects of extreme flood of Brenta River (Italy) are examinated. The flooded areas (urban and rural areas) are identified and a mitigation solution based on a diversion channel flowing into Venice Lagoon is proposed. We show that this solution strongly reduces the flood risk in the downstream areas and can provide an important sediment source to the Venice Lagoon. Finally, preliminary results of the sediment dispersion in the Venice Lagoon are presented.

On the Origins of Vortex Shedding in Two-dimensional Incompressible Flows

PubMed Central

Boghosian, M. E.; Cassel, K. W.

2016-01-01

An exegesis of a novel mechanism leading to vortex splitting and subsequent shedding that is valid for two-dimensional incompressible, inviscid or viscous, and external or internal or wall-bounded flows, is detailed in this research. The mechanism, termed the Vortex-Shedding Mechanism (VSM), is simple and intuitive, requiring only two coincident conditions in the flow: (1) the existence of a location with zero momentum and (2) the presence of a net force having a positive divergence. Numerical solutions of several model problems illustrate causality of the VSM. Moreover, the VSM criteria is proved to be a necessary and sufficient condition for a vortex splitting event in any two-dimensional, incompressible flow. The VSM is shown to exist in several canonical problems including the external flow past a circular cylinder. Suppression of the von Kármán vortex street is demonstrated for Reynolds numbers of 100 and 400 by mitigating the VSM. PMID:27795617

On the Origins of Vortex Shedding in Two-dimensional Incompressible Flows.

PubMed

Boghosian, M E; Cassel, K W

2016-12-01

An exegesis of a novel mechanism leading to vortex splitting and subsequent shedding that is valid for two-dimensional incompressible, inviscid or viscous, and external or internal or wall-bounded flows, is detailed in this research. The mechanism, termed the Vortex-Shedding Mechanism (VSM), is simple and intuitive, requiring only two coincident conditions in the flow: (1) the existence of a location with zero momentum and (2) the presence of a net force having a positive divergence. Numerical solutions of several model problems illustrate causality of the VSM. Moreover, the VSM criteria is proved to be a necessary and sufficient condition for a vortex splitting event in any two-dimensional, incompressible flow. The VSM is shown to exist in several canonical problems including the external flow past a circular cylinder. Suppression of the von Kármán vortex street is demonstrated for Reynolds numbers of 100 and 400 by mitigating the VSM.

Bernstein-Greene-Kruskal Modes in a Three-Dimensional Plasma

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

Ng, C.S.; Bhattacharjee, A.

2005-12-09

Bernstein-Greene-Kruskal modes in a three-dimensional (3D) unmagnetized plasma are constructed. It is shown that 3D solutions that depend only on energy do not exist. However, 3D solutions that depend on energy and additional constants of motion (such as angular momentum) do exist. Exact analytical as well as numerical solutions are constructed assuming spherical symmetry, and their properties are contrasted with those of 1D solutions. Possible extensions to solutions with cylindrical symmetry with or without a finite magnetic guide field are discussed.

Panel summary report

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

Gutjahr, A.L.; Kincaid, C.T.; Mercer, J.W.

1987-04-01

The objective of this report is to summarize the various modeling approaches that were used to simulate solute transport in a variably saturated emission. In particular, the technical strengths and weaknesses of each approach are discussed, and conclusions and recommendations for future studies are made. Five models are considered: (1) one-dimensional analytical and semianalytical solutions of the classical deterministic convection-dispersion equation (van Genuchten, Parker, and Kool, this report ); (2) one-dimensional simulation using a continuous-time Markov process (Knighton and Wagenet, this report); (3) one-dimensional simulation using the time domain method and the frequency domain method (Duffy and Al-Hassan, this report);more » (4) one-dimensional numerical approach that combines a solution of the classical deterministic convection-dispersion equation with a chemical equilibrium speciation model (Cederberg, this report); and (5) three-dimensional numerical solution of the classical deterministic convection-dispersion equation (Huyakorn, Jones, Parker, Wadsworth, and White, this report). As part of the discussion, the input data and modeling results are summarized. The models were used in a data analysis mode, as opposed to a predictive mode. Thus, the following discussion will concentrate on the data analysis aspects of model use. Also, all the approaches were similar in that they were based on a convection-dispersion model of solute transport. Each discussion addresses the modeling approaches in the order listed above.« less

A Numerical Investigation of the Burnett Equations Based on the Second Law

NASA Technical Reports Server (NTRS)

Comeaux, Keith A.; Chapman, Dean R.; MacCormack, Robert W.; Edwards, Thomas A. (Technical Monitor)

1995-01-01

The Burnett equations have been shown to potentially violate the second law of thermodynamics. The objective of this investigation is to correlate the numerical problems experienced by the Burnett equations to the negative production of entropy. The equations have had a long history of numerical instability to small wavelength disturbances. Recently, Zhong corrected the instability problem and made solutions attainable for one dimensional shock waves and hypersonic blunt bodies. Difficulties still exist when attempting to solve hypersonic flat plate boundary layers and blunt body wake flows, however. Numerical experiments will include one-dimensional shock waves, quasi-one dimensional nozzles, and expanding Prandlt-Meyer flows and specifically examine the entropy production for these cases.

The PAC-MAN model: Benchmark case for linear acoustics in computational physics

NASA Astrophysics Data System (ADS)

Ziegelwanger, Harald; Reiter, Paul

2017-10-01

Benchmark cases in the field of computational physics, on the one hand, have to contain a certain complexity to test numerical edge cases and, on the other hand, require the existence of an analytical solution, because an analytical solution allows the exact quantification of the accuracy of a numerical simulation method. This dilemma causes a need for analytical sound field formulations of complex acoustic problems. A well known example for such a benchmark case for harmonic linear acoustics is the ;Cat's Eye model;, which describes the three-dimensional sound field radiated from a sphere with a missing octant analytically. In this paper, a benchmark case for two-dimensional (2D) harmonic linear acoustic problems, viz., the ;PAC-MAN model;, is proposed. The PAC-MAN model describes the radiated and scattered sound field around an infinitely long cylinder with a cut out sector of variable angular width. While the analytical calculation of the 2D sound field allows different angular cut-out widths and arbitrarily positioned line sources, the computational cost associated with the solution of this problem is similar to a 1D problem because of a modal formulation of the sound field in the PAC-MAN model.

Exact Solutions for Wind-Driven Coastal Upwelling and Downwelling over Sloping Topography

NASA Astrophysics Data System (ADS)

Choboter, P.; Duke, D.; Horton, J.; Sinz, P.

2009-12-01

The dynamics of wind-driven coastal upwelling and downwelling are studied using a simplified dynamical model. Exact solutions are examined as a function of time and over a family of sloping topographies. Assumptions in the two-dimensional model include a frictionless ocean interior below the surface Ekman layer, and no alongshore dependence of the variables; however, dependence in the cross-shore and vertical directions is retained. Additionally, density and alongshore momentum are advected by the cross-shore velocity in order to maintain thermal wind. The time-dependent initial-value problem is solved with constant initial stratification and no initial alongshore flow. An alongshore pressure gradient is added to allow the cross-shore flow to be geostrophically balanced far from shore. Previously, this model has been used to study upwelling over flat-bottom and sloping topographies, but the novel feature in this work is the discovery of exact solutions for downwelling. These exact solutions are compared to numerical solutions from a primitive-equation ocean model, based on the Princeton Ocean Model, configured in a similar two-dimensional geometry. Many typical features of the evolution of density and velocity during downwelling are displayed by the analytical model.

An investigation on a two-dimensional problem of Mode-I crack in a thermoelastic medium

NASA Astrophysics Data System (ADS)

Kant, Shashi; Gupta, Manushi; Shivay, Om Namha; Mukhopadhyay, Santwana

2018-04-01

In this work, we consider a two-dimensional dynamical problem of an infinite space with finite linear Mode-I crack and employ a recently proposed heat conduction model: an exact heat conduction with a single delay term. The thermoelastic medium is taken to be homogeneous and isotropic. However, the boundary of the crack is subjected to a prescribed temperature and stress distributions. The Fourier and Laplace transform techniques are used to solve the problem. Mathematical modeling of the present problem reduces the solution of the problem into the solution of a system of four dual integral equations. The solution of these equations is equivalent to the solution of the Fredholm's integral equation of the first kind which has been solved by using the regularization method. Inverse Laplace transform is carried out by using the Bellman method, and we obtain the numerical solution for all the physical field variables in the physical domain. Results are shown graphically, and we highlight the effects of the presence of crack in the behavior of thermoelastic interactions inside the medium in the present context, and its results are compared with the results of the thermoelasticity of type-III.

Modification of a method-of-characteristics solute-transport model to incorporate decay and equilibrium-controlled sorption or ion exchange

USGS Publications Warehouse

Goode, D.J.; Konikow, Leonard F.

1989-01-01

The U.S. Geological Survey computer model of two-dimensional solute transport and dispersion in ground water (Konikow and Bredehoeft, 1978) has been modified to incorporate the following types of chemical reactions: (1) first-order irreversible rate-reaction, such as radioactive decay; (2) reversible equilibrium-controlled sorption with linear, Freundlich, or Langmuir isotherms; and (3) reversible equilibrium-controlled ion exchange for monovalent or divalent ions. Numerical procedures are developed to incorporate these processes in the general solution scheme that uses method-of- characteristics with particle tracking for advection and finite-difference methods for dispersion. The first type of reaction is accounted for by an exponential decay term applied directly to the particle concentration. The second and third types of reactions are incorporated through a retardation factor, which is a function of concentration for nonlinear cases. The model is evaluated and verified by comparison with analytical solutions for linear sorption and decay, and by comparison with other numerical solutions for nonlinear sorption and ion exchange.

Influence of heterogeneity on second-kind self-similar solutions for viscous gravity currents

DOE PAGES

Zheng, Zhong; Christov, Ivan  C.; Stone, Howard  A.

2014-05-01

We report experimental, theoretical and numerical results on the effects of horizontal heterogeneities on the propagation of viscous gravity currents. We use two geometries to highlight these effects: (a) a horizontal channel (or crack) whose gap thickness varies as a power-law function of the streamwise coordinate; (b) a heterogeneous porous medium whose permeability and porosity have power-law variations. We demonstrate that two types of self-similar behaviours emerge as a result of horizontal heterogeneity: (a) a first-kind self-similar solution is found using dimensional analysis (scaling) for viscous gravity currents that propagate away from the origin (a point of zero permeability); (b)more » a second-kind self-similar solution is found using a phase-plane analysis for viscous gravity currents that propagate toward the origin. These theoretical predictions, obtained using the ideas of self-similar intermediate asymptotics, are compared with experimental results and numerical solutions of the governing partial differential equation developed under the lubrication approximation. All three results are found to be in good agreement.« less

Documentation of the Surface-Water Routing (SWR1) Process for modeling surface-water flow with the U.S. Geological Survey Modular Ground-Water Model (MODFLOW-2005)

USGS Publications Warehouse

Hughes, Joseph D.; Langevin, Christian D.; Chartier, Kevin L.; White, Jeremy T.

2012-01-01

A flexible Surface-Water Routing (SWR1) Process that solves the continuity equation for one-dimensional and two-dimensional surface-water flow routing has been developed for the U.S. Geological Survey three-dimensional groundwater model, MODFLOW-2005. Simple level- and tilted-pool reservoir routing and a diffusive-wave approximation of the Saint-Venant equations have been implemented. Both methods can be implemented in the same model and the solution method can be simplified to represent constant-stage elements that are functionally equivalent to the standard MODFLOW River or Drain Package boundary conditions. A generic approach has been used to represent surface-water features (reaches) and allows implementation of a variety of geometric forms. One-dimensional geometric forms include rectangular, trapezoidal, and irregular cross section reaches to simulate one-dimensional surface-water features, such as canals and streams. Two-dimensional geometric forms include reaches defined using specified stage-volume-area-perimeter (SVAP) tables and reaches covering entire finite-difference grid cells to simulate two-dimensional surface-water features, such as wetlands and lakes. Specified SVAP tables can be used to represent reaches that are smaller than the finite-difference grid cell (for example, isolated lakes), or reaches that cannot be represented accurately using the defined top of the model. Specified lateral flows (which can represent point and distributed flows) and stage-dependent rainfall and evaporation can be applied to each reach. The SWR1 Process can be used with the MODFLOW Unsaturated Zone Flow (UZF1) Package to permit dynamic simulation of runoff from the land surface to specified reaches. Surface-water/groundwater interactions in the SWR1 Process are mathematically defined to be a function of the difference between simulated stages and groundwater levels, and the specific form of the reach conductance equation used in each reach. Conductance can be specified directly or calculated as a function of the simulated wetted perimeter and defined reach bed hydraulic properties, or as a weighted combination of both reach bed hydraulic properties and horizontal hydraulic conductivity. Each reach can be explicitly coupled to a single specific groundwater-model layer or coupled to multiple groundwater-model layers based on the reach geometry and groundwater-model layer elevations in the row and column containing the reach. Surface-water flow between reservoirs is simulated using control structures. Surface-water flow between reaches, simulated by the diffusive-wave approximation, can also be simulated using control structures. A variety of control structures have been included in the SWR1 Process and include (1) excess-volume structures, (2) uncontrolled-discharge structures, (3) pumps, (4) defined stage-discharge relations, (5) culverts, (6) fixed- or movable-crest weirs, and (7) fixed or operable gated spillways. Multiple control structures can be implemented in individual reaches and are treated as composite flow structures. Solution of the continuity equation at the reach-group scale (a single reach or a user-defined collection of individual reaches) is achieved using exact Newton methods with direct solution methods or exact and inexact Newton methods with Krylov sub-space methods. Newton methods have been used in the SWR1 Process because of their ability to solve nonlinear problems. Multiple SWR1 time steps can be simulated for each MODFLOW time step, and a simple adaptive time-step algorithm, based on user-specified rainfall, stage, flow, or convergence constraints, has been implemented to better resolve surface-water response. A simple linear- or sigmoid-depth scaling approach also has been implemented to account for increased bed roughness at small surface-water depths and to increase numerical stability. A line-search algorithm also has been included to improve the quality of the Newton-step upgrade vector, if possible. The SWR1 Process has been benchmarked against one- and two-dimensional numerical solutions from existing one- and two-dimensional numerical codes that solve the dynamic-wave approximation of the Saint-Venant equations. Two-dimensional solutions test the ability of the SWR1 Process to simulate the response of a surface-water system to (1) steady flow conditions for an inclined surface (solution of Manning's equation), and (2) transient inflow and rainfall for an inclined surface. The one-dimensional solution tests the ability of the SWR1 Process to simulate a looped network with multiple upstream inflows and several control structures. The SWR1 Process also has been compared to a level-pool reservoir solution. A synthetic test problem was developed to evaluate a number of different SWR1 solution options and simulate surface-water/groundwater interaction. The solution approach used in the SWR1 Process may not be applicable for all surface-water/groundwater problems. The SWR1 Process is best suited for modeling long-term changes (days to years) in surface-water and groundwater flow. Use of the SWR1 Process is not recommended for modeling the transient exchange of water between streams and aquifers when local and convective acceleration and other secondary effects (for example, wind and Coriolis forces) are substantial. Dam break evaluations and two-dimensional evaluations of spatially extensive domains are examples where acceleration terms and secondary effects would be significant, respectively.

An unconditionally stable method for numerically solving solar sail spacecraft equations of motion

NASA Astrophysics Data System (ADS)

Karwas, Alex

Solar sails use the endless supply of the Sun's radiation to propel spacecraft through space. The sails use the momentum transfer from the impinging solar radiation to provide thrust to the spacecraft while expending zero fuel. Recently, the first solar sail spacecraft, or sailcraft, named IKAROS completed a successful mission to Venus and proved the concept of solar sail propulsion. Sailcraft experimental data is difficult to gather due to the large expenses of space travel, therefore, a reliable and accurate computational method is needed to make the process more efficient. Presented in this document is a new approach to simulating solar sail spacecraft trajectories. The new method provides unconditionally stable numerical solutions for trajectory propagation and includes an improved physical description over other methods. The unconditional stability of the new method means that a unique numerical solution is always determined. The improved physical description of the trajectory provides a numerical solution and time derivatives that are continuous throughout the entire trajectory. The error of the continuous numerical solution is also known for the entire trajectory. Optimal control for maximizing thrust is also provided within the framework of the new method. Verification of the new approach is presented through a mathematical description and through numerical simulations. The mathematical description provides details of the sailcraft equations of motion, the numerical method used to solve the equations, and the formulation for implementing the equations of motion into the numerical solver. Previous work in the field is summarized to show that the new approach can act as a replacement to previous trajectory propagation methods. A code was developed to perform the simulations and it is also described in this document. Results of the simulations are compared to the flight data from the IKAROS mission. Comparison of the two sets of data show that the new approach is capable of accurately simulating sailcraft motion. Sailcraft and spacecraft simulations are compared to flight data and to other numerical solution techniques. The new formulation shows an increase in accuracy over a widely used trajectory propagation technique. Simulations for two-dimensional, three-dimensional, and variable attitude trajectories are presented to show the multiple capabilities of the new technique. An element of optimal control is also part of the new technique. An additional equation is added to the sailcraft equations of motion that maximizes thrust in a specific direction. A technical description and results of an example optimization problem are presented. The spacecraft attitude dynamics equations take the simulation a step further by providing control torques using the angular rate and acceleration outputs of the numerical formulation.

Numerical solutions of 2-D multi-stage rotor/stator unsteady flow interactions

NASA Astrophysics Data System (ADS)

Yang, R.-J.; Lin, S.-J.

1991-01-01

The Rai method of single-stage rotor/stator flow interaction is extended to handle multistage configurations. In this study, a two-dimensional Navier-Stokes multi-zone approach was used to investigate unsteady flow interactions within two multistage axial turbines. The governing equations are solved by an iterative, factored, implicit finite-difference, upwind algorithm. Numerical accuracy is checked by investigating the effect of time step size, the effect of subiteration in the Newton-Raphson technique, and the effect of full viscous versus thin-layer approximation. Computer results compared well with experimental data. Unsteady flow interactions, wake cutting, and the associated evolution of vortical entities are discussed.

Unsteady flow model for circulation-control airfoils

NASA Technical Reports Server (NTRS)

Rao, B. M.

1979-01-01

An analysis and a numerical lifting surface method are developed for predicting the unsteady airloads on two-dimensional circulation control airfoils in incompressible flow. The analysis and the computer program are validated by correlating the computed unsteady airloads with test data and also with other theoretical solutions. Additionally, a mathematical model for predicting the bending-torsion flutter of a two-dimensional airfoil (a reference section of a wing or rotor blade) and a computer program using an iterative scheme are developed. The flutter program has a provision for using the CC airfoil airloads program or the Theodorsen hard flap solution to compute the unsteady lift and moment used in the flutter equations. The adopted mathematical model and the iterative scheme are used to perform a flutter analysis of a typical CC rotor blade reference section. The program seems to work well within the basic assumption of the incompressible flow.

Finite volume model for two-dimensional shallow environmental flow

USGS Publications Warehouse

Simoes, F.J.M.

2011-01-01

This paper presents the development of a two-dimensional, depth integrated, unsteady, free-surface model based on the shallow water equations. The development was motivated by the desire of balancing computational efficiency and accuracy by selective and conjunctive use of different numerical techniques. The base framework of the discrete model uses Godunov methods on unstructured triangular grids, but the solution technique emphasizes the use of a high-resolution Riemann solver where needed, switching to a simpler and computationally more efficient upwind finite volume technique in the smooth regions of the flow. Explicit time marching is accomplished with strong stability preserving Runge-Kutta methods, with additional acceleration techniques for steady-state computations. A simplified mass-preserving algorithm is used to deal with wet/dry fronts. Application of the model is made to several benchmark cases that show the interplay of the diverse solution techniques.

Numerical Solution of Time-Dependent Problems with a Fractional-Power Elliptic Operator

NASA Astrophysics Data System (ADS)

Vabishchevich, P. N.

2018-03-01

A time-dependent problem in a bounded domain for a fractional diffusion equation is considered. The first-order evolution equation involves a fractional-power second-order elliptic operator with Robin boundary conditions. A finite-element spatial approximation with an additive approximation of the operator of the problem is used. The time approximation is based on a vector scheme. The transition to a new time level is ensured by solving a sequence of standard elliptic boundary value problems. Numerical results obtained for a two-dimensional model problem are presented.

Numerical solution for the temperature distribution in a cooled guide vane blade of a radial gas turbine

NASA Technical Reports Server (NTRS)

Hosny, W. M.; Tabakoff, W.

1977-01-01

A two dimensional finite difference numerical technique is presented to determine the temperature distribution of an internal cooled blade of radial turbine guide vanes. A simple convection cooling is assumed inside the guide vane blade. Such cooling has relatively small cooling effectiveness at the leading edge and at the trailing edge. Heat transfer augmentation in these critical areas may be achieved by using impingement jets and film cooling. A computer program is written in FORTRAN IV for IBM 370/165 computer.

Numerical realization of the variational method for generating self-trapped beams.

PubMed

Duque, Erick I; Lopez-Aguayo, Servando; Malomed, Boris A

2018-03-19

We introduce a numerical variational method based on the Rayleigh-Ritz optimization principle for predicting two-dimensional self-trapped beams in nonlinear media. This technique overcomes the limitation of the traditional variational approximation in performing analytical Lagrangian integration and differentiation. Approximate soliton solutions of a generalized nonlinear Schrödinger equation are obtained, demonstrating robustness of the beams of various types (fundamental, vortices, multipoles, azimuthons) in the course of their propagation. The algorithm offers possibilities to produce more sophisticated soliton profiles in general nonlinear models.

Numerical Schemes and Computational Studies for Dynamically Orthogonal Equations (Multidisciplinary Simulation, Estimation, and Assimilation Systems: Reports in Ocean Science and Engineering)

DTIC Science & Technology

2011-08-01

heat transfers [49, 52]. However, the DO method has not yet been applied to Boussinesq flows, and the numerical challenges of the DO decomposition for...used a PCE scheme to study mixing in a two-dimensional (2D) microchannel and improved the efficiency of their solution scheme by decoupling the...to several Navier-Stokes flows and their stochastic dynamics has been studied, including mean-mode and mode-mode energy transfers for 2D flows and

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.

Electro-magneto interaction in fractional Green-Naghdi thermoelastic solid with a cylindrical cavity

NASA Astrophysics Data System (ADS)

Ezzat, M. A.; El-Bary, A. A.

2018-01-01

A unified mathematical model of Green-Naghdi's thermoelasticty theories (GN), based on fractional time-derivative of heat transfer is constructed. The model is applied to solve a one-dimensional problem of a perfect conducting unbounded body with a cylindrical cavity subjected to sinusoidal pulse heating in the presence of an axial uniform magnetic field. Laplace transform techniques are used to get the general analytical solutions in Laplace domain, and the inverse Laplace transforms based on Fourier expansion techniques are numerically implemented to obtain the numerical solutions in time domain. Comparisons are made with the results predicted by the two theories. The effects of the fractional derivative parameter on thermoelastic fields for different theories are discussed.

A diffusion model of protected population on bilocal habitat with generalized resource

NASA Astrophysics Data System (ADS)

Vasilyev, Maxim D.; Trofimtsev, Yuri I.; Vasilyeva, Natalya V.

2017-11-01

A model of population distribution in a two-dimensional area divided by an ecological barrier, i.e. the boundaries of natural reserve, is considered. Distribution of the population is defined by diffusion, directed migrations and areal resource. The exchange of specimens occurs between two parts of the habitat. The mathematical model is presented in the form of a boundary value problem for a system of non-linear parabolic equations with variable parameters of diffusion and growth function. The splitting space variables, sweep method and simple iteration methods were used for the numerical solution of a system. A set of programs was coded in Python. Numerical simulation results for the two-dimensional unsteady non-linear problem are analyzed in detail. The influence of migration flow coefficients and functions of natural birth/death ratio on the distributions of population densities is investigated. The results of the research would allow to describe the conditions of the stable and sustainable existence of populations in bilocal habitat containing the protected and non-protected zones.

Spatial solitons and stability in the one-dimensional and the two-dimensional generalized nonlinear Schrödinger equation with fourth-order diffraction and parity-time-symmetric potentials

NASA Astrophysics Data System (ADS)

Tiofack, C. G. L.; Ndzana, F., II; Mohamadou, A.; Kofane, T. C.

2018-03-01

We investigate the existence and stability of solitons in parity-time (PT )-symmetric optical media characterized by a generic complex hyperbolic refractive index distribution and fourth-order diffraction (FOD). For the linear case, we demonstrate numerically that the FOD parameter can alter the PT -breaking points. For nonlinear cases, the exact analytical expressions of the localized modes are obtained both in one- and two-dimensional nonlinear Schrödinger equations with self-focusing and self-defocusing Kerr nonlinearity. The effect of FOD on the stability structure of these localized modes is discussed with the help of linear stability analysis followed by the direct numerical simulation of the governing equation. Examples of stable and unstable solutions are given. The transverse power flow density associated with these localized modes is also discussed. It is found that the relative strength of the FOD coefficient can utterly change the direction of the power flow, which may be used to control the energy exchange among gain or loss regions.

Stabilized finite element methods to simulate the conductances of ion channels

NASA Astrophysics Data System (ADS)

Tu, Bin; Xie, Yan; Zhang, Linbo; Lu, Benzhuo

2015-03-01

We have previously developed a finite element simulator, ichannel, to simulate ion transport through three-dimensional ion channel systems via solving the Poisson-Nernst-Planck equations (PNP) and Size-modified Poisson-Nernst-Planck equations (SMPNP), and succeeded in simulating some ion channel systems. However, the iterative solution between the coupled Poisson equation and the Nernst-Planck equations has difficulty converging for some large systems. One reason we found is that the NP equations are advection-dominated diffusion equations, which causes troubles in the usual FE solution. The stabilized schemes have been applied to compute fluids flow in various research fields. However, they have not been studied in the simulation of ion transport through three-dimensional models based on experimentally determined ion channel structures. In this paper, two stabilized techniques, the SUPG and the Pseudo Residual-Free Bubble function (PRFB) are introduced to enhance the numerical robustness and convergence performance of the finite element algorithm in ichannel. The conductances of the voltage dependent anion channel (VDAC) and the anthrax toxin protective antigen pore (PA) are simulated to validate the stabilization techniques. Those two stabilized schemes give reasonable results for the two proteins, with decent agreement with both experimental data and Brownian dynamics (BD) simulations. For a variety of numerical tests, it is found that the simulator effectively avoids previous numerical instability after introducing the stabilization methods. Comparison based on our test data set between the two stabilized schemes indicates both SUPG and PRFB have similar performance (the latter is slightly more accurate and stable), while SUPG is relatively more convenient to implement.

Three-dimensional unsteady lifting surface theory in the subsonic range

NASA Technical Reports Server (NTRS)

Kuessner, H. G.

1985-01-01

The methods of the unsteady lifting surface theory are surveyed. Linearized Euler's equations are simplified by means of a Galileo-Lorentz transformation and a Laplace transformation so that the time and the compressibility of the fluid are limited to two constants. The solutions to this simplified problem are represented as integrals with a differential nucleus; these results in tolerance conditions, for which any exact solution must suffice. It is shown that none of the existing three-dimensional lifting surface theories in subsonic range satisfy these conditions. An oscillating elliptic lifting surface which satisfies the tolerance conditions is calculated through the use of Lame's functions. Numerical examples are calculated for the borderline cases of infinitely stretched elliptic lifting surfaces and of circular lifting surfaces. Out of the harmonic solutions any such temporal changes of the down current are calculated through the use of an inverse Laplace transformation.

Reaction mechanisms of methylene-blue degradation in three-dimensionally integrated micro-solution plasma

NASA Astrophysics Data System (ADS)

Shirafuji, Tatsuru; Ishida, Yodai; Nomura, Ayano; Hayashi, Yui; Goto, Motonobu

2017-06-01

We have performed matrix-assisted laser desorption ionization time-of-flight (MALDI-TOF) mass spectrometry (MS) on methylene-blue aqueous solutions treated with three-dimensionally integrated micro-solution plasma, in which we have acquired the time evolution of mass spectra as a function of treatment time. The time evolution of mass spectral peak intensities for major detected species has clearly indicated that the parent methylene-blue molecules are degraded through consecutive reactions. The primary reaction is the oxidation of the parent molecules. The oxidized species still have two benzene rings in the parent molecules. The secondary reactions are the separation of the oxidized species and the formation of compounds with one benzene ring. We have also performed the numerical fitting of the time evolution of the mass spectral peak intensities, the results of which have indicated that we must assume additional primary reactions before the primary oxidation for better agreement with experimental results.

Salty popcorn in a homogeneous low-dimensional toy model of holographic QCD

NASA Astrophysics Data System (ADS)

Elliot-Ripley, Matthew

2017-04-01

Recently, a homogeneous ansatz has been used to study cold dense nuclear matter in the Sakai-Sugimoto model of holographic QCD. To justify this homogeneous approximation we here investigate a homogeneous ansatz within a low-dimensional toy version of Sakai-Sugimoto to study finite baryon density configurations and compare it to full numerical solutions. We find the ansatz corresponds to enforcing a dyon salt arrangement in which the soliton solutions are split into half-soliton layers. Within this ansatz we find analogues of the proposed baryonic popcorn transitions, in which solutions split into multiple layers in the holographic direction. The homogeneous results are found to qualitatively match the full numerical solutions, lending confidence to the homogeneous approximations of the full Sakai-Sugimoto model. In addition, we find exact compact solutions in the high density, flat space limit which demonstrate the existence of further popcorn transitions to three layers and beyond.

A new flux conserving Newton's method scheme for the two-dimensional, steady Navier-Stokes equations

NASA Technical Reports Server (NTRS)

Scott, James R.; Chang, Sin-Chung

1993-01-01

A new numerical method is developed for the solution of the two-dimensional, steady Navier-Stokes equations. The method that is presented differs in significant ways from the established numerical methods for solving the Navier-Stokes equations. The major differences are described. First, the focus of the present method is on satisfying flux conservation in an integral formulation, rather than on simulating conservation laws in their differential form. Second, the present approach provides a unified treatment of the dependent variables and their unknown derivatives. All are treated as unknowns together to be solved for through simulating local and global flux conservation. Third, fluxes are balanced at cell interfaces without the use of interpolation or flux limiters. Fourth, flux conservation is achieved through the use of discrete regions known as conservation elements and solution elements. These elements are not the same as the standard control volumes used in the finite volume method. Fifth, the discrete approximation obtained on each solution element is a functional solution of both the integral and differential form of the Navier-Stokes equations. Finally, the method that is presented is a highly localized approach in which the coupling to nearby cells is only in one direction for each spatial coordinate, and involves only the immediately adjacent cells. A general third-order formulation for the steady, compressible Navier-Stokes equations is presented, and then a Newton's method scheme is developed for the solution of incompressible, low Reynolds number channel flow. It is shown that the Jacobian matrix is nearly block diagonal if the nonlinear system of discrete equations is arranged approximately and a proper pivoting strategy is used. Numerical results are presented for Reynolds numbers of 100, 1000, and 2000. Finally, it is shown that the present scheme can resolve the developing channel flow boundary layer using as few as six to ten cells per channel width, depending on the Reynolds number.

An Implicit Characteristic Based Method for Electromagnetics

NASA Technical Reports Server (NTRS)

Beggs, John H.; Briley, W. Roger

2001-01-01

An implicit characteristic-based approach for numerical solution of Maxwell's time-dependent curl equations in flux conservative form is introduced. This method combines a characteristic based finite difference spatial approximation with an implicit lower-upper approximate factorization (LU/AF) time integration scheme. This approach is advantageous for three-dimensional applications because the characteristic differencing enables a two-factor approximate factorization that retains its unconditional stability in three space dimensions, and it does not require solution of tridiagonal systems. Results are given both for a Fourier analysis of stability, damping and dispersion properties, and for one-dimensional model problems involving propagation and scattering for free space and dielectric materials using both uniform and nonuniform grids. The explicit Finite Difference Time Domain Method (FDTD) algorithm is used as a convenient reference algorithm for comparison. The one-dimensional results indicate that for low frequency problems on a highly resolved uniform or nonuniform grid, this LU/AF algorithm can produce accurate solutions at Courant numbers significantly greater than one, with a corresponding improvement in efficiency for simulating a given period of time. This approach appears promising for development of dispersion optimized LU/AF schemes for three dimensional applications.

Rapid and efficient mixing in a slip-driven three-dimensional flow in a rectangular channel

NASA Astrophysics Data System (ADS)

Pacheco, J. Rafael; Ping Chen, Kang; Hayes, Mark A.

2006-08-01

A method for generating mixing in an electroosmotic flow of an electrolytic solution in a three-dimensional channel is proposed. When the width-to-height aspect ratio of the channel cross-section is large, mixing of a blob of a solute in a slip-driven three-dimensional flow in a rectangular channel can be used to model and assess the effectiveness of this method. It is demonstrated through numerical simulations that under certain operating conditions, rapid and efficient mixing can be achieved. Future investigation will include the solution of the exact equations and experimentation.

Modeling radium and radon transport through soil and vegetation

USGS Publications Warehouse

Kozak, J.A.; Reeves, H.W.; Lewis, B.A.

2003-01-01

A one-dimensional flow and transport model was developed to describe the movement of two fluid phases, gas and water, within a porous medium and the transport of 226Ra and 222Rn within and between these two phases. Included in this model is the vegetative uptake of water and aqueous 226Ra and 222Rn that can be extracted from the soil via the transpiration stream. The mathematical model is formulated through a set of phase balance equations and a set of species balance equations. Mass exchange, sink terms and the dependence of physical properties upon phase composition couple the two sets of equations. Numerical solution of each set, with iteration between the sets, is carried out leading to a set-iterative compositional model. The Petrov-Galerkin finite element approach is used to allow for upstream weighting if required for a given simulation. Mass lumping improves solution convergence and stability behavior. The resulting numerical model was applied to four problems and was found to produce accurate, mass conservative solutions when compared to published experimental and numerical results and theoretical column experiments. Preliminary results suggest that the model can be used as an investigative tool to determine the feasibility of phytoremediating radium and radon-contaminated soil. ?? 2003 Elsevier Science B.V. All rights reserved.

A low-dimensional analogue of holographic baryons

NASA Astrophysics Data System (ADS)

Bolognesi, Stefano; Sutcliffe, Paul

2014-04-01

Baryons in holographic QCD correspond to topological solitons in the bulk. The most prominent example is the Sakai-Sugimoto model, where the bulk soliton in the five-dimensional spacetime of AdS-type can be approximated by the flat space self-dual Yang-Mills instanton with a small size. Recently, the validity of this approximation has been verified by comparison with the numerical field theory solution. However, multi-solitons and solitons with finite density are currently beyond numerical field theory computations. Various approximations have been applied to investigate these important issues and have led to proposals for finite density configurations that include dyonic salt and baryonic popcorn. Here we introduce and investigate a low-dimensional analogue of the Sakai-Sugimoto model, in which the bulk soliton can be approximated by a flat space sigma model instanton. The bulk theory is a baby Skyrme model in a three-dimensional spacetime with negative curvature. The advantage of the lower-dimensional theory is that numerical simulations of multi-solitons and finite density solutions can be performed and compared with flat space instanton approximations. In particular, analogues of dyonic salt and baryonic popcorn configurations are found and analysed.

Gravity and large black holes in Randall-Sundrum II braneworlds.

PubMed

Figueras, Pau; Wiseman, Toby

2011-08-19

We show how to construct low energy solutions to the Randall-Sundrum II (RSII) model by using an associated five-dimensional anti-de Sitter space (AdS(5)) and/or four-dimensional conformal field theory (CFT(4)) problem. The RSII solution is given as a perturbation of the AdS(5)-CFT(4) solution, with the perturbation parameter being the radius of curvature of the brane metric compared to the AdS length ℓ. The brane metric is then a specific perturbation of the AdS(5)-CFT(4) boundary metric. For low curvatures the RSII solution reproduces 4D general relativity on the brane. Recently, AdS(5)-CFT(4) solutions with a 4D Schwarzschild boundary metric were numerically constructed. We modify the boundary conditions to numerically construct large RSII static black holes with radius up to ~20ℓ. For a large radius, the RSII solutions are indeed close to the associated AdS(5)-CFT(4) solution. © 2011 American Physical Society

Numerical Analysis of an H 1-Galerkin Mixed Finite Element Method for Time Fractional Telegraph Equation

PubMed Central

Wang, Jinfeng; Zhao, Meng; Zhang, Min; Liu, Yang; Li, Hong

2014-01-01

We discuss and analyze an H 1-Galerkin mixed finite element (H 1-GMFE) method to look for the numerical solution of time fractional telegraph equation. We introduce an auxiliary variable to reduce the original equation into lower-order coupled equations and then formulate an H 1-GMFE scheme with two important variables. We discretize the Caputo time fractional derivatives using the finite difference methods and approximate the spatial direction by applying the H 1-GMFE method. Based on the discussion on the theoretical error analysis in L 2-norm for the scalar unknown and its gradient in one dimensional case, we obtain the optimal order of convergence in space-time direction. Further, we also derive the optimal error results for the scalar unknown in H 1-norm. Moreover, we derive and analyze the stability of H 1-GMFE scheme and give the results of a priori error estimates in two- or three-dimensional cases. In order to verify our theoretical analysis, we give some results of numerical calculation by using the Matlab procedure. PMID:25184148

Numerical solution of inverse scattering for near-field optics.

PubMed

Bao, Gang; Li, Peijun

2007-06-01

A novel regularized recursive linearization method is developed for a two-dimensional inverse medium scattering problem that arises in near-field optics, which reconstructs the scatterer of an inhomogeneous medium located on a substrate from data accessible through photon scanning tunneling microscopy experiments. Based on multiple frequency scattering data, the method starts from the Born approximation corresponding to weak scattering at a low frequency, and each update is obtained by continuation on the wavenumber from solutions of one forward problem and one adjoint problem of the Helmholtz equation.

A Runge-Kutta discontinuous finite element method for high speed flows

NASA Technical Reports Server (NTRS)

Bey, Kim S.; Oden, J. T.

1991-01-01

A Runge-Kutta discontinuous finite element method is developed for hyperbolic systems of conservation laws in two space variables. The discontinuous Galerkin spatial approximation to the conservation laws results in a system of ordinary differential equations which are marched in time using Runge-Kutta methods. Numerical results for the two-dimensional Burger's equation show that the method is (p+1)-order accurate in time and space, where p is the degree of the polynomial approximation of the solution within an element and is capable of capturing shocks over a single element without oscillations. Results for this problem also show that the accuracy of the solution in smooth regions is unaffected by the local projection and that the accuracy in smooth regions increases as p increases. Numerical results for the Euler equations show that the method captures shocks without oscillations and with higher resolution than a first-order scheme.

A Two-Dimensional Helmholtz Equation Solution for the Multiple Cavity Scattering Problem

DTIC Science & Technology

2013-02-01

obtained by using the block Gauss – Seidel iterative meth- od. To show the convergence of the iterative method, we define the error between two...models to the general multiple cavity setting. Numerical examples indicate that the convergence of the Gauss – Seidel iterative method depends on the...variational approach. A block Gauss – Seidel iterative method is introduced to solve the cou- pled system of the multiple cavity scattering problem, where

On the relationship between kinetic and fluid formalisms for convection in the inner magnetosphere

NASA Astrophysics Data System (ADS)

Song, Yang; Sazykin, Stanislav; Wolf, Richard A.

2008-08-01

In the inner magnetosphere, the plasma flows are mostly slow compared to thermal or Alfvén speeds, but the convection is far away from the ideal magnetohydrodynamic regime since the gradient/curvature drifts become significant. Both kinetic (Wolf, 1983) and two-fluid (Peymirat and Fontaine, 1994; Heinemann, 1999) formalisms have been used to describe plasma dynamics, but it is not fully understood how they relate to each other. We explore the relations among kinetic, fluid, and recently developed "average" (Liu, 2006) models in an attempt to find the simplest yet realistic way to describe the convection. First, we prove analytically that the model of (Liu, 2006), when closed with the assumption of a Maxwellian distribution, is equivalent to the fluid model of (Heinemann, 1999). Second, we analyze the transport of both one-dimensional and two-dimensional Gaussian-shaped blob of hot plasma. For the kinetic case, it is known that the time evolution of such a blob is gradual spreading in time. For the fluid case, Heinemann and Wolf (2001a, 2001b) showed that in a one-dimensional idealized case, the blob separates into two drifting at different speeds. We present a fully nonlinear solution of this case, confirming this behavior but demonstrating what appears to be a shocklike steepening of the faster drifting secondary blob. A new, more realistic two-dimensional example using the dipole geometry with a uniform electric field confirms the one-dimensional solutions. Implications for the numerical simulations of magnetospheric dynamics are discussed.

Finite Deformation of Magnetoelastic Film

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

Barham, Matthew Ian

2011-05-31

A nonlinear two-dimensional theory is developed for thin magnetoelastic lms capable of large deformations. This is derived directly from three-dimensional theory. Signi cant simpli cations emerge in the descent from three dimensions to two, permitting the self eld generated by the body to be computed a posteriori. The model is specialized to isotropic elastomers with two material models. First weak magnetization is investigated leading to a free energy where magnetization and deformation are un-coupled. The second closely couples the magnetization and deformation. Numerical solutions are obtained to equilibrium boundary-value problems in which the membrane is subjected to lateral pressure andmore » an applied magnetic eld. An instability is inferred and investigated for the weak magnetization material model.« less

Two-dimensional solitons in conservative and parity-time-symmetric triple-core waveguides with cubic-quintic nonlinearity

NASA Astrophysics Data System (ADS)

Feijoo, David; Zezyulin, Dmitry A.; Konotop, Vladimir V.

2015-12-01

We analyze a system of three two-dimensional nonlinear Schrödinger equations coupled by linear terms and with the cubic-quintic (focusing-defocusing) nonlinearity. We consider two versions of the model: conservative and parity-time (PT ) symmetric. These models describe triple-core nonlinear optical waveguides, with balanced gain and losses in the PT -symmetric case. We obtain families of soliton solutions and discuss their stability. The latter study is performed using a linear stability analysis and checked with direct numerical simulations of the evolutional system of equations. Stable solitons are found in the conservative and PT -symmetric cases. Interactions and collisions between the conservative and PT -symmetric solitons are briefly investigated, as well.

Least-squares finite element solutions for three-dimensional backward-facing step flow

NASA Technical Reports Server (NTRS)

Jiang, Bo-Nan; Hou, Lin-Jun; Lin, Tsung-Liang

1993-01-01

Comprehensive numerical solutions of the steady state incompressible viscous flow over a three-dimensional backward-facing step up to Re equals 800 are presented. The results are obtained by the least-squares finite element method (LSFEM) which is based on the velocity-pressure-vorticity formulation. The computed model is of the same size as that of Armaly's experiment. Three-dimensional phenomena are observed even at low Reynolds number. The calculated values of the primary reattachment length are in good agreement with experimental results.

Modelling Detailed-Chemistry Effects on Turbulent Diffusion Flames using a Parallel Solution-Adaptive Scheme

NASA Astrophysics Data System (ADS)

Jha, Pradeep Kumar

Capturing the effects of detailed-chemistry on turbulent combustion processes is a central challenge faced by the numerical combustion community. However, the inherent complexity and non-linear nature of both turbulence and chemistry require that combustion models rely heavily on engineering approximations to remain computationally tractable. This thesis proposes a computationally efficient algorithm for modelling detailed-chemistry effects in turbulent diffusion flames and numerically predicting the associated flame properties. The cornerstone of this combustion modelling tool is the use of parallel Adaptive Mesh Refinement (AMR) scheme with the recently proposed Flame Prolongation of Intrinsic low-dimensional manifold (FPI) tabulated-chemistry approach for modelling complex chemistry. The effect of turbulence on the mean chemistry is incorporated using a Presumed Conditional Moment (PCM) approach based on a beta-probability density function (PDF). The two-equation k-w turbulence model is used for modelling the effects of the unresolved turbulence on the mean flow field. The finite-rate of methane-air combustion is represented here by using the GRI-Mech 3.0 scheme. This detailed mechanism is used to build the FPI tables. A state of the art numerical scheme based on a parallel block-based solution-adaptive algorithm has been developed to solve the Favre-averaged Navier-Stokes (FANS) and other governing partial-differential equations using a second-order accurate, fully-coupled finite-volume formulation on body-fitted, multi-block, quadrilateral/hexahedral mesh for two-dimensional and three-dimensional flow geometries, respectively. A standard fourth-order Runge-Kutta time-marching scheme is used for time-accurate temporal discretizations. Numerical predictions of three different diffusion flames configurations are considered in the present work: a laminar counter-flow flame; a laminar co-flow diffusion flame; and a Sydney bluff-body turbulent reacting flow. Comparisons are made between the predicted results of the present FPI scheme and Steady Laminar Flamelet Model (SLFM) approach for diffusion flames. The effects of grid resolution on the predicted overall flame solutions are also assessed. Other non-reacting flows have also been considered to further validate other aspects of the numerical scheme. The present schemes predict results which are in good agreement with published experimental results and reduces the computational cost involved in modelling turbulent diffusion flames significantly, both in terms of storage and processing time.

The modified semi-discrete two-dimensional Toda lattice with self-consistent sources

NASA Astrophysics Data System (ADS)

Gegenhasi

2017-07-01

In this paper, we derive the Grammian determinant solutions to the modified semi-discrete two-dimensional Toda lattice equation, and then construct the semi-discrete two-dimensional Toda lattice equation with self-consistent sources via source generation procedure. The algebraic structure of the resulting coupled modified differential-difference equation is clarified by presenting its Grammian determinant solutions and Casorati determinant solutions. As an application of the Grammian determinant and Casorati determinant solution, the explicit one-soliton and two-soliton solution of the modified semi-discrete two-dimensional Toda lattice equation with self-consistent sources are given. We also construct another form of the modified semi-discrete two-dimensional Toda lattice equation with self-consistent sources which is the Bäcklund transformation for the semi-discrete two-dimensional Toda lattice equation with self-consistent sources.

Numerical simulation of the hydrodynamical combustion to strange quark matter

NASA Astrophysics Data System (ADS)

Niebergal, Brian; Ouyed, Rachid; Jaikumar, Prashanth

2010-12-01

We present results from a numerical solution to the burning of neutron matter inside a cold neutron star into stable u,d,s quark matter. Our method solves hydrodynamical flow equations in one dimension with neutrino emission from weak equilibrating reactions, and strange quark diffusion across the burning front. We also include entropy change from heat released in forming the stable quark phase. Our numerical results suggest burning front laminar speeds of 0.002-0.04 times the speed of light, much faster than previous estimates derived using only a reactive-diffusive description. Analytic solutions to hydrodynamical jump conditions with a temperature-dependent equation of state agree very well with our numerical findings for fluid velocities. The most important effect of neutrino cooling is that the conversion front stalls at lower density (below ≈2 times saturation density). In a two-dimensional setting, such rapid speeds and neutrino cooling may allow for a flame wrinkle instability to develop, possibly leading to detonation.

Numerical modelling techniques of soft soil improvement via stone columns: A brief review

NASA Astrophysics Data System (ADS)

Zukri, Azhani; Nazir, Ramli

2018-04-01

There are a number of numerical studies on stone column systems in the literature. Most of the studies found were involved with two-dimensional analysis of the stone column behaviour, while only a few studies used three-dimensional analysis. The most popular software utilised in those studies was Plaxis 2D and 3D. Other types of software that used for numerical analysis are DIANA, EXAMINE, ZSoil, ABAQUS, ANSYS, NISA, GEOSTUDIO, CRISP, TOCHNOG, CESAR, GEOFEM (2D & 3D), FLAC, and FLAC 3. This paper will review the methodological approaches to model stone column numerically, both in two-dimensional and three-dimensional analyses. The numerical techniques and suitable constitutive model used in the studies will also be discussed. In addition, the validation methods conducted were to verify the numerical analysis conducted will be presented. This review paper also serves as a guide for junior engineers through the applicable procedures and considerations when constructing and running a two or three-dimensional numerical analysis while also citing numerous relevant references.

Navier-Stokes solution on the CYBER-203 by a pseudospectral technique

NASA Technical Reports Server (NTRS)

Lambiotte, J. J.; Hussaini, M. Y.; Bokhari, S.; Orszag, S. A.

1983-01-01

A three-level, time-split, mixed spectral/finite difference method for the numerical solution of the three-dimensional, compressible Navier-Stokes equations has been developed and implemented on the Control Data Corporation (CDC) CYBER-203. This method uses a spectral representation for the flow variables in the streamwise and spanwise coordinates, and central differences in the normal direction. The five dependent variables are interleaved one horizontal plane at a time and the array of their values at the grid points of each horizontal plane is a typical vector in the computation. The code is organized so as to require, per time step, a single forward-backward pass through the entire data base. The one-and two-dimensional Fast Fourier Transforms are performed using software especially developed for the CYBER-203.

Embedding methods for the steady Euler equations

NASA Technical Reports Server (NTRS)

Chang, S. H.; Johnson, G. M.

1983-01-01

An approach to the numerical solution of the steady Euler equations is to embed the first-order Euler system in a second-order system and then to recapture the original solution by imposing additional boundary conditions. Initial development of this approach and computational experimentation with it were previously based on heuristic physical reasoning. This has led to the construction of a relaxation procedure for the solution of two-dimensional steady flow problems. The theoretical justification for the embedding approach is addressed. It is proven that, with the appropriate choice of embedding operator and additional boundary conditions, the solution to the embedded system is exactly the one to the original Euler equations. Hence, solving the embedded version of the Euler equations will not produce extraneous solutions.

An approximate theoretical method for modeling the static thrust performance of non-axisymmetric two-dimensional convergent-divergent nozzles. M.S. Thesis - George Washington Univ.

NASA Technical Reports Server (NTRS)

Hunter, Craig A.

1995-01-01

An analytical/numerical method has been developed to predict the static thrust performance of non-axisymmetric, two-dimensional convergent-divergent exhaust nozzles. Thermodynamic nozzle performance effects due to over- and underexpansion are modeled using one-dimensional compressible flow theory. Boundary layer development and skin friction losses are calculated using an approximate integral momentum method based on the classic karman-Polhausen solution. Angularity effects are included with these two models in a computational Nozzle Performance Analysis Code, NPAC. In four different case studies, results from NPAC are compared to experimental data obtained from subscale nozzle testing to demonstrate the capabilities and limitations of the NPAC method. In several cases, the NPAC prediction matched experimental gross thrust efficiency data to within 0.1 percent at a design NPR, and to within 0.5 percent at off-design conditions.

ANALYZING NUMERICAL ERRORS IN DOMAIN HEAT TRANSPORT MODELS USING THE CVBEM.

USGS Publications Warehouse

Hromadka, T.V.; ,

1985-01-01

Besides providing an exact solution for steady-state heat conduction processes (Laplace Poisson equations), the CVBEM (complex variable boundary element method) can be used for the numerical error analysis of domain model solutions. For problems where soil water phase change latent heat effects dominate the thermal regime, heat transport can be approximately modeled as a time-stepped steady-state condition in the thawed and frozen regions, respectively. The CVBEM provides an exact solution of the two-dimensional steady-state heat transport problem, and also provides the error in matching the prescribed boundary conditions by the development of a modeling error distribution or an approximative boundary generation. This error evaluation can be used to develop highly accurate CVBEM models of the heat transport process, and the resulting model can be used as a test case for evaluating the precision of domain models based on finite elements or finite differences.

Dissipative discrete breathers: periodic, quasiperiodic, chaotic, and mobile.

PubMed

Martínez, P J; Meister, M; Floría, L M; Falo, F

2003-06-01

The properties of discrete breathers in dissipative one-dimensional lattices of nonlinear oscillators subject to periodic driving forces are reviewed. We focus on oscillobreathers in the Frenkel-Kontorova chain and rotobreathers in a ladder of Josephson junctions. Both types of exponentially localized solutions are easily obtained numerically using adiabatic continuation from the anticontinuous limit. Linear stability (Floquet) analysis allows the characterization of different types of bifurcations experienced by periodic discrete breathers. Some of these bifurcations produce nonperiodic localized solutions, namely, quasiperiodic and chaotic discrete breathers, which are generally impossible as exact solutions in Hamiltonian systems. Within a certain range of parameters, propagating breathers occur as attractors of the dissipative dynamics. General features of these excitations are discussed and the Peierls-Nabarro barrier is addressed. Numerical scattering experiments with mobile breathers reveal the existence of two-breather bound states and allow a first glimpse at the intricate phenomenology of these special multibreather configurations. (c) 2003 American Institute of Physics.

Huygens' inspired multi-pendulum setups: Experiments and stability analysis

NASA Astrophysics Data System (ADS)

Hoogeboom, F. N.; Pogromsky, A. Y.; Nijmeijer, H.

2016-11-01

This paper examines synchronization of a set of metronomes placed on a lightweight foam platform. Two configurations of the set of metronomes are considered: a row setup containing one-dimensional coupling and a cross setup containing two-dimensional coupling. Depending on the configuration and coupling between the metronomes, i.e., the platform parameters, in- and/or anti-phase synchronized behavior is observed in the experiments. To explain this behavior, mathematical models of a metronome and experimental setups have been derived and used in a local stability analysis. It is numerically and experimentally demonstrated that varying the coupling parameters for both configurations has a significant influence on the stability of the synchronized solutions.

Dynamics of bright-bright solitons in Bose-Einstein condensate with Raman-induced one-dimensional spin-orbit coupling

NASA Astrophysics Data System (ADS)

Wen, Lin; Zhang, Xiao-Fei; Hu, Ai-Yuan; Zhou, Jing; Yu, Peng; Xia, Lei; Sun, Qing; Ji, An-Chun

2018-03-01

We investigate the dynamics of bright-bright solitons in one-dimensional two-component Bose-Einstein condensates with Raman-induced spin-orbit coupling, via the variational approximation and the numerical simulation of Gross-Pitaevskii equations. For the uniform system without trapping potential, we obtain two population balanced stationary solitons. By performing the linear stability analysis, we find a Goldstone eigenmode and an oscillation eigenmode around these stationary solitons. Moreover, we derive a general dynamical solution to describe the center-of-mass motion and spin evolution of the solitons under the action of spin-orbit coupling. The effects of a harmonic trap have also been discussed.

Programming of the complex logarithm function in the solution of the cracked anisotropic plate loaded by a point force

NASA Astrophysics Data System (ADS)

Zaal, K. J. J. M.

1991-06-01

In programming solutions of complex function theory, the complex logarithm function is replaced by the complex logarithmic function, introducing a discontinuity along the branch cut into the programmed solution which was not present in the mathematical solution. Recently, Liaw and Kamel presented their solution of the infinite anisotropic centrally cracked plate loaded by an arbitrary point force, which they used as Green's function in a boundary element method intended to evaluate the stress intensity factor at the tip of a crack originating from an elliptical home. Their solution may be used as Green's function of many more numerical methods involving anisotropic elasticity. In programming applications of Liaw and Kamel's solution, the standard definition of the logarithmic function with the branch cut at the nonpositive real axis cannot provide a reliable computation of the displacement field for Liaw and Kamel's solution. Either the branch cut should be redefined outside the domain of the logarithmic function, after proving that the domain is limited to a part of the plane, or the logarithmic function should be defined on its Riemann surface. A two dimensional line fractal can provide the link between all mesh points on the plane essential to evaluate the logarithm function on its Riemann surface. As an example, a two dimensional line fractal is defined for a mesh once used by Erdogan and Arin.

On the effects of mass and momentum transfer from droplets impacting on steady two-dimensional rimming flow in a horizontal cylinder

NASA Astrophysics Data System (ADS)

Williams, J.; Hibberd, S.; Power, H.; Riley, D. S.

2012-05-01

Motivated by applications in aero-engines, steady two-dimensional thin-film flow on the inside of a circular cylinder is studied when the film surface is subject to mass and momentum transfer from impacting droplets. Asymptotic analysis is used systematically to identify distinguished limits that incorporate these transfer effects at leading order and to provide a new mathematical model. Applying both analytical and numerical approaches to the model, a set of stable steady, two-dimensional solutions that fit within the rational framework is determined. A number of these solutions feature steep fronts and associated recirculating pools, which are undesirable in an aeroengine since oil may be stripped away from the steep fronts when there is a core flow external to the film, and recirculation may lead to oil degradation. The model, however, provides a means of investigating whether the formation of the steep fronts on the film surface and of internal recirculation pools can be delayed, or inhibited altogether, by designing jets to deliver prescribed distributions of oil droplets or by the judicious siting of oil sinks. Moreover, by studying pathlines, oil-residence times can be predicted and systems optimized.

A dimensionally split Cartesian cut cell method for hyperbolic conservation laws

NASA Astrophysics Data System (ADS)

Gokhale, Nandan; Nikiforakis, Nikos; Klein, Rupert

2018-07-01

We present a dimensionally split method for solving hyperbolic conservation laws on Cartesian cut cell meshes. The approach combines local geometric and wave speed information to determine a novel stabilised cut cell flux, and we provide a full description of its three-dimensional implementation in the dimensionally split framework of Klein et al. [1]. The convergence and stability of the method are proved for the one-dimensional linear advection equation, while its multi-dimensional numerical performance is investigated through the computation of solutions to a number of test problems for the linear advection and Euler equations. When compared to the cut cell flux of Klein et al., it was found that the new flux alleviates the problem of oscillatory boundary solutions produced by the former at higher Courant numbers, and also enables the computation of more accurate solutions near stagnation points. Being dimensionally split, the method is simple to implement and extends readily to multiple dimensions.

A Semi-Implicit, Three-Dimensional Model for Estuarine Circulation

USGS Publications Warehouse

Smith, Peter E.

2006-01-01

A semi-implicit, finite-difference method for the numerical solution of the three-dimensional equations for circulation in estuaries is presented and tested. The method uses a three-time-level, leapfrog-trapezoidal scheme that is essentially second-order accurate in the spatial and temporal numerical approximations. The three-time-level scheme is shown to be preferred over a two-time-level scheme, especially for problems with strong nonlinearities. The stability of the semi-implicit scheme is free from any time-step limitation related to the terms describing vertical diffusion and the propagation of the surface gravity waves. The scheme does not rely on any form of vertical/horizontal mode-splitting to treat the vertical diffusion implicitly. At each time step, the numerical method uses a double-sweep method to transform a large number of small tridiagonal equation systems and then uses the preconditioned conjugate-gradient method to solve a single, large, five-diagonal equation system for the water surface elevation. The governing equations for the multi-level scheme are prepared in a conservative form by integrating them over the height of each horizontal layer. The layer-integrated volumetric transports replace velocities as the dependent variables so that the depth-integrated continuity equation that is used in the solution for the water surface elevation is linear. Volumetric transports are computed explicitly from the momentum equations. The resulting method is mass conservative, efficient, and numerically accurate.

An improved exceedance theory for combined random stresses

NASA Technical Reports Server (NTRS)

Lester, H. C.

1974-01-01

An extension is presented of Rice's classic solution for the exceedances of a constant level by a single random process to its counterpart for an n-dimensional vector process. An interaction boundary, analogous to the constant level considered by Rice for the one-dimensional case, is assumed in the form of a hypersurface. The theory for the numbers of boundary exceedances is developed by using a joint statistical approach which fully accounts for all cross-correlation effects. An exact expression is derived for the n-dimensional exceedance density function, which is valid for an arbitrary interaction boundary. For application to biaxial states of combined random stress, the general theory is reduced to the two-dimensional case. An elliptical stress interaction boundary is assumed and the exact expression for the density function is presented. The equations are expressed in a format which facilitates calculating the exceedances by numerically evaluating a line integral. The behavior of the density function for the two-dimensional case is briefly discussed.

Fundamental Solution For The Self-healing Fracture Pulse

NASA Astrophysics Data System (ADS)

Nielsen, S.; Madariaga, R.

We find the analytical solution for a fundamental fracture mode in the form of a self- similar, self-healing pulse. The existence of such a fracture mode was strongly sug- gested by recent numerical findings but, to our knwledge, no formal proof had been proposed up to date. We present a two dimensional, anti-plane solution for fixed rup- ture and healing velocities, that satisfies both wave equation and stress conditions; we argue that such a solution is plausible even in the absence of rate-weakening in the friction, as an alternative to the classic crack solution. In practice, the impulsive mode rather than the expanding crack mode is selected depending on details of fracture initiation, and is therafter self-maintained. We discuss stress concentration, fracture energy, rupture velocity and compare them to the case of a crack. The analytical study is complemented by various numerical examples and comparisons. On more general grounds, we argue that an infinity of marginally stable fracture modes may exist other than the crack solution or the impulseive fracture described here.

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

Giunta, G.; Belouettar, S.

In this paper, the static response of three-dimensional beams made of functionally graded materials is investigated through a family of hierarchical one-dimensional finite elements. A wide variety of elements is proposed differing by the kinematic formulation and the number of nodes per elements along the beam axis. Elements’ stiffness matrix and load vector are derived in a unified nuclear form that does not depend upon the a priori expansion order over the cross-section nor the finite element approximation along the beam axis. Results are validated towards three-dimensional finite element models as well as equivalent Navier-type analytical solutions. The numerical investigationsmore » show that accurate and efficient solutions (when compared with full three-dimensional FEM solutions) can be obtained by the proposed family of hierarchical one-dimensional elements’ family.« less

On the prediction of far field computational aeroacoustics of advanced propellers

NASA Technical Reports Server (NTRS)

Jaeger, Stephen M.; Korkan, Kenneth D.

1990-01-01

A numerical method for determining the acoustic far field generated by a high-speed subsonic aircraft propeller was developed. The approach used in this method was to generate the entire three-dimensional pressure field about the propeller (using an Euler flowfield solver) and then to apply a solution of the wave equation on a cylindrical surface enveloping the propeller. The method is applied to generate the three-dimensional flowfield between two blades of an advanced propeller. The results are compared with experimental data obtained in a wind-tunnel test at a Mach number of 0.6.

Aquifer response to stream-stage and recharge variations. I. Analytical step-response functions

USGS Publications Warehouse

Moench, A.F.; Barlow, P.M.

2000-01-01

Laplace transform step-response functions are presented for various homogeneous confined and leaky aquifer types and for anisotropic, homogeneous unconfined aquifers interacting with perennial streams. Flow is one-dimensional, perpendicular to the stream in the confined and leaky aquifers, and two-dimensional in a plane perpendicular to the stream in the water-table aquifers. The stream is assumed to penetrate the full thickness of the aquifer. The aquifers may be semi-infinite or finite in width and may or may not be bounded at the stream by a semipervious streambank. The solutions are presented in a unified manner so that mathematical relations among the various aquifer configurations are clearly demonstrated. The Laplace transform solutions are inverted numerically to obtain the real-time step-response functions for use in the convolution (or superposition) integral. To maintain linearity in the case of unconfined aquifers, fluctuations in the elevation of the water table are assumed to be small relative to the saturated thickness, and vertical flow into or out of the zone above the water table is assumed to occur instantaneously. Effects of hysteresis in the moisture distribution above the water table are therefore neglected. Graphical comparisons of the new solutions are made with known closed-form solutions.Laplace transform step-response functions are presented for various homogeneous confined and leaky aquifer types and for anisotropic, homogeneous unconfined aquifers interacting with perennial streams. Flow is one-dimensional, perpendicular to the stream in the confined and leaky aquifers, and two-dimensional in a plane perpendicular to the stream in the water-table aquifers. The stream is assumed to penetrate the full thickness of the aquifer. The aquifers may be semi-infinite or finite in width and may or may not be bounded at the stream by a semipervious streambank. The solutions are presented in a unified manner so that mathematical relations among the various aquifer configurations are clearly demonstrated. The Laplace transform solutions are inverted numerically to obtain the real-time step-response functions for use in the convolution (or superposition) integral. To maintain linearity in the case of unconfined aquifers, fluctuations in the elevation of the water table are assumed to be small relative to the saturated thickness, and vertical flow into or out of the zone above the water table is assumed to occur instantaneously. Effects of hysteresis in the moisture distribution above the water table are therefore neglected. Graphical comparisons of the new solutions are made with known closed-form solutions.

Solitons in Bose-Einstein Condensates

NASA Astrophysics Data System (ADS)

Carr, Lincoln D.

2003-05-01

The stationary form, dynamical properties, and experimental criteria for creation of matter-wave bright and dark solitons, both singly and in trains, are studied numerically and analytically in the context of Bose-Einstein condensates [1]. The full set of stationary solutions in closed analytic form to the mean field model in the quasi-one-dimensional regime, which is a nonlinear Schrodinger equation equally relevant in nonlinear optics, is developed under periodic and box boundary conditions [2]. These solutions are extended numerically into the two and three dimensional regimes, where it is shown that dark solitons can be used to create vortex-anti-vortex pairs under realistic conditions. Specific experimental prescriptions for creating viable dark and bright solitons in the quasi-one-dimensional regime are provided. These analytic methods are then extended to treat the nonlinear Schrodinger equation with a generalized lattice potential, which models a Bose-Einstein condensate trapped in the potential generated by a standing light wave. A novel solution family is developed and stability criterion are presented. Experiments which successfully carried out these ideas are briefly discussed [3]. [1] Dissertation research completed at the University of Washington Physics Department under the advisorship of Prof. William P. Reinhardt. [2] L. D. Carr, C. W. Clark, and W. P. Reinhardt, Phys. Rev. A v. 62 p. 063610-1--10 and Phys. Rev. A v.62, p.063611-1--10 (2000). [3] L. Khaykovich, F. Schreck, T. Bourdel, J. Cubizolles, G. Ferrari, L. D. Carr, Y. Castin, and C. Salomon, Science v. 296, p.1290--1293 (2002).

Numerical calculation of the internal flow field in a centrifugal compressor impeller

NASA Technical Reports Server (NTRS)

Walitt, L.; Harp, J. L., Jr.; Liu, C. Y.

1975-01-01

An iterative numerical method has been developed for the calculation of steady, three-dimensional, viscous, compressible flow fields in centrifugal compressor impellers. The computer code, which embodies the method, solves the steady three dimensional, compressible Navier-Stokes equations in rotating, curvilinear coordinates. The solution takes place on blade-to-blade surfaces of revolution which move from the hub to the shroud during each iteration.

A novel method for correcting scanline-observational bias of discontinuity orientation

PubMed Central

Huang, Lei; Tang, Huiming; Tan, Qinwen; Wang, Dingjian; Wang, Liangqing; Ez Eldin, Mutasim A. M.; Li, Changdong; Wu, Qiong

2016-01-01

Scanline observation is known to introduce an angular bias into the probability distribution of orientation in three-dimensional space. In this paper, numerical solutions expressing the functional relationship between the scanline-observational distribution (in one-dimensional space) and the inherent distribution (in three-dimensional space) are derived using probability theory and calculus under the independence hypothesis of dip direction and dip angle. Based on these solutions, a novel method for obtaining the inherent distribution (also for correcting the bias) is proposed, an approach which includes two procedures: 1) Correcting the cumulative probabilities of orientation according to the solutions, and 2) Determining the distribution of the corrected orientations using approximation methods such as the one-sample Kolmogorov-Smirnov test. The inherent distribution corrected by the proposed method can be used for discrete fracture network (DFN) modelling, which is applied to such areas as rockmass stability evaluation, rockmass permeability analysis, rockmass quality calculation and other related fields. To maximize the correction capacity of the proposed method, the observed sample size is suggested through effectiveness tests for different distribution types, dispersions and sample sizes. The performance of the proposed method and the comparison of its correction capacity with existing methods are illustrated with two case studies. PMID:26961249

Thermoelastic damping in thin microrings with two-dimensional heat conduction

NASA Astrophysics Data System (ADS)

Fang, Yuming; Li, Pu

2015-05-01

Accurate determination of thermoelastic damping (TED) is very challenging in the design of micro-resonators. Microrings are widely used in many micro-resonators. In the past, to model the TED effect on the microrings, some analytical models have been developed. However, in the previous works, the heat conduction within the microring is modeled by using the one-dimensional approach. The governing equation for heat conduction is solved only for the one-dimensional heat conduction along the radial thickness of the microring. This paper presents a simple analytical model for TED in microrings. The two-dimensional heat conduction over the thermoelastic temperature gradients along the radial thickness and the circumferential direction are considered in the present model. A two-dimensional heat conduction equation is developed. The solution of the equation is represented by the product of an assumed sine series along the radial thickness and an assumed trigonometric series along the circumferential direction. The analytical results obtained by the present 2-D model show a good agreement with the numerical (FEM) results. The limitations of the previous 1-D model are assessed.

Finite-volume application of high order ENO schemes to multi-dimensional boundary-value problems

NASA Technical Reports Server (NTRS)

Casper, Jay; Dorrepaal, J. Mark

1990-01-01

The finite volume approach in developing multi-dimensional, high-order accurate essentially non-oscillatory (ENO) schemes is considered. In particular, a two dimensional extension is proposed for the Euler equation of gas dynamics. This requires a spatial reconstruction operator that attains formal high order of accuracy in two dimensions by taking account of cross gradients. Given a set of cell averages in two spatial variables, polynomial interpolation of a two dimensional primitive function is employed in order to extract high-order pointwise values on cell interfaces. These points are appropriately chosen so that correspondingly high-order flux integrals are obtained through each interface by quadrature, at each point having calculated a flux contribution in an upwind fashion. The solution-in-the-small of Riemann's initial value problem (IVP) that is required for this pointwise flux computation is achieved using Roe's approximate Riemann solver. Issues to be considered in this two dimensional extension include the implementation of boundary conditions and application to general curvilinear coordinates. Results of numerical experiments are presented for qualitative and quantitative examination. These results contain the first successful application of ENO schemes to boundary value problems with solid walls.

A generalized volumetric dispersion model for a class of two-phase separation/reaction: finite difference solutions

NASA Astrophysics Data System (ADS)

Siripatana, Chairat; Thongpan, Hathaikarn; Promraksa, Arwut

2017-03-01

This article explores a volumetric approach in formulating differential equations for a class of engineering flow problems involving component transfer within or between two phases. In contrast to conventional formulation which is based on linear velocities, this work proposed a slightly different approach based on volumetric flow-rate which is essentially constant in many industrial processes. In effect, many multi-dimensional flow problems found industrially can be simplified into multi-component or multi-phase but one-dimensional flow problems. The formulation is largely generic, covering counter-current, concurrent or batch, fixed and fluidized bed arrangement. It was also intended to use for start-up, shut-down, control and steady state simulation. Since many realistic and industrial operation are dynamic with variable velocity and porosity in relation to position, analytical solutions are rare and limited to only very simple cases. Thus we also provide a numerical solution using Crank-Nicolson finite difference scheme. This solution is inherently stable as tested against a few cases published in the literature. However, it is anticipated that, for unconfined flow or non-constant flow-rate, traditional formulation should be applied.

Second-order small-disturbance solutions for hypersonic flow over power-law bodies

NASA Technical Reports Server (NTRS)

Townsend, J. C.

1975-01-01

Similarity solutions were found which give the adiabatic flow of an ideal gas about two-dimensional and axisymmetric power-law bodies at infinite Mach number to second order in the body slenderness parameter. The flow variables were expressed as a sum of zero-order and perturbation similarity functions for which the axial variations in the flow equations separated out. The resulting similarity equations were integrated numerically. The solutions, which are universal functions, are presented in graphic and tabular form. To avoid a singularity in the calculations, the results are limited to body power-law exponents greater than about 0.85 for the two-dimensional case and 0.75 for the axisymmetric case. Because of the entropy layer induced by the nose bluntness (for power-law bodies other than cones and wedges), only the pressure function is valid at the body surface. The similarity results give excellent agreement with the exact solutions for inviscid flow over wedges and cones having half-angles up to about 20 deg. They give good agreement with experimental shock-wave shapes and surface-pressure distributions for 3/4-power axisymmetric bodies, considering that Mach number and boundary-layer displacement effects are not included in the theory.

Ensemble solute transport in two-dimensional operator-scaling random fields

NASA Astrophysics Data System (ADS)

Monnig, Nathan D.; Benson, David A.; Meerschaert, Mark M.

2008-02-01

Motivated by field measurements of aquifer hydraulic conductivity (K), recent techniques were developed to construct anisotropic fractal random fields in which the scaling, or self-similarity parameter, varies with direction and is defined by a matrix. Ensemble numerical results are analyzed for solute transport through these two-dimensional "operator-scaling" fractional Brownian motion ln(K) fields. Both the longitudinal and transverse Hurst coefficients, as well as the "radius of isotropy" are important to both plume growth rates and the timing and duration of breakthrough. It is possible to create operator-scaling fractional Brownian motion fields that have more "continuity" or stratification in the direction of transport. The effects on a conservative solute plume are continually faster-than-Fickian growth rates, highly non-Gaussian shapes, and a heavier tail early in the breakthrough curve. Contrary to some analytic stochastic theories for monofractal K fields, the plume growth rates never exceed A. Mercado's (1967) purely stratified aquifer growth rate of plume apparent dispersivity proportional to mean distance. Apparent superstratified growth must be the result of other demonstrable factors, such as initial plume size.

Time-Accurate Local Time Stepping and High-Order Time CESE Methods for Multi-Dimensional Flows Using Unstructured Meshes

NASA Technical Reports Server (NTRS)

Chang, Chau-Lyan; Venkatachari, Balaji Shankar; Cheng, Gary

2013-01-01

With the wide availability of affordable multiple-core parallel supercomputers, next generation numerical simulations of flow physics are being focused on unsteady computations for problems involving multiple time scales and multiple physics. These simulations require higher solution accuracy than most algorithms and computational fluid dynamics codes currently available. This paper focuses on the developmental effort for high-fidelity multi-dimensional, unstructured-mesh flow solvers using the space-time conservation element, solution element (CESE) framework. Two approaches have been investigated in this research in order to provide high-accuracy, cross-cutting numerical simulations for a variety of flow regimes: 1) time-accurate local time stepping and 2) highorder CESE method. The first approach utilizes consistent numerical formulations in the space-time flux integration to preserve temporal conservation across the cells with different marching time steps. Such approach relieves the stringent time step constraint associated with the smallest time step in the computational domain while preserving temporal accuracy for all the cells. For flows involving multiple scales, both numerical accuracy and efficiency can be significantly enhanced. The second approach extends the current CESE solver to higher-order accuracy. Unlike other existing explicit high-order methods for unstructured meshes, the CESE framework maintains a CFL condition of one for arbitrarily high-order formulations while retaining the same compact stencil as its second-order counterpart. For large-scale unsteady computations, this feature substantially enhances numerical efficiency. Numerical formulations and validations using benchmark problems are discussed in this paper along with realistic examples.

Non-linear instability analysis of the two-dimensional Navier-Stokes equation: The Taylor-Green vortex problem

NASA Astrophysics Data System (ADS)

Sengupta, Tapan K.; Sharma, Nidhi; Sengupta, Aditi

2018-05-01

An enstrophy-based non-linear instability analysis of the Navier-Stokes equation for two-dimensional (2D) flows is presented here, using the Taylor-Green vortex (TGV) problem as an example. This problem admits a time-dependent analytical solution as the base flow, whose instability is traced here. The numerical study of the evolution of the Taylor-Green vortices shows that the flow becomes turbulent, but an explanation for this transition has not been advanced so far. The deviation of the numerical solution from the analytical solution is studied here using a high accuracy compact scheme on a non-uniform grid (NUC6), with the fourth-order Runge-Kutta method. The stream function-vorticity (ψ, ω) formulation of the governing equations is solved here in a periodic square domain with four vortices at t = 0. Simulations performed at different Reynolds numbers reveal that numerical errors in computations induce a breakdown of symmetry and simultaneous fragmentation of vortices. It is shown that the actual physical instability is triggered by the growth of disturbances and is explained by the evolution of disturbance mechanical energy and enstrophy. The disturbance evolution equations have been traced by looking at (a) disturbance mechanical energy of the Navier-Stokes equation, as described in the work of Sengupta et al., "Vortex-induced instability of an incompressible wall-bounded shear layer," J. Fluid Mech. 493, 277-286 (2003), and (b) the creation of rotationality via the enstrophy transport equation in the work of Sengupta et al., "Diffusion in inhomogeneous flows: Unique equilibrium state in an internal flow," Comput. Fluids 88, 440-451 (2013).

Asymptotic and spectral analysis of the gyrokinetic-waterbag integro-differential operator in toroidal geometry

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

Besse, Nicolas, E-mail: Nicolas.Besse@oca.eu; Institut Jean Lamour, UMR CNRS/UL 7198, Université de Lorraine, BP 70239 54506 Vandoeuvre-lès-Nancy Cedex; Coulette, David, E-mail: David.Coulette@ipcms.unistra.fr

2016-08-15

Achieving plasmas with good stability and confinement properties is a key research goal for magnetic fusion devices. The underlying equations are the Vlasov–Poisson and Vlasov–Maxwell (VPM) equations in three space variables, three velocity variables, and one time variable. Even in those somewhat academic cases where global equilibrium solutions are known, studying their stability requires the analysis of the spectral properties of the linearized operator, a daunting task. We have identified a model, for which not only equilibrium solutions can be constructed, but many of their stability properties are amenable to rigorous analysis. It uses a class of solution to themore » VPM equations (or to their gyrokinetic approximations) known as waterbag solutions which, in particular, are piecewise constant in phase-space. It also uses, not only the gyrokinetic approximation of fast cyclotronic motion around magnetic field lines, but also an asymptotic approximation regarding the magnetic-field-induced anisotropy: the spatial variation along the field lines is taken much slower than across them. Together, these assumptions result in a drastic reduction in the dimensionality of the linearized problem, which becomes a set of two nested one-dimensional problems: an integral equation in the poloidal variable, followed by a one-dimensional complex Schrödinger equation in the radial variable. We show here that the operator associated to the poloidal variable is meromorphic in the eigenparameter, the pulsation frequency. We also prove that, for all but a countable set of real pulsation frequencies, the operator is compact and thus behaves mostly as a finite-dimensional one. The numerical algorithms based on such ideas have been implemented in a companion paper [D. Coulette and N. Besse, “Numerical resolution of the global eigenvalue problem for gyrokinetic-waterbag model in toroidal geometry” (submitted)] and were found to be surprisingly close to those for the original gyrokinetic-Vlasov equations. The purpose of the present paper is to make these new ideas accessible to two readerships: applied mathematicians and plasma physicists.« less

Numerical methods for axisymmetric and 3D nonlinear beams

NASA Astrophysics Data System (ADS)

Pinton, Gianmarco F.; Trahey, Gregg E.

2005-04-01

Time domain algorithms that solve the Khokhlov--Zabolotzskaya--Kuznetsov (KZK) equation are described and implemented. This equation represents the propagation of finite amplitude sound beams in a homogenous thermoviscous fluid for axisymmetric and fully three dimensional geometries. In the numerical solution each of the terms is considered separately and the numerical methods are compared with known solutions. First and second order operator splitting are used to combine the separate terms in the KZK equation and their convergence is examined.

Advanced development of BEM for elastic and inelastic dynamic analysis of solids

NASA Technical Reports Server (NTRS)

Banerjee, P. K.; Ahmad, S.; Wang, H. C.

1989-01-01

Direct Boundary Element formulations and their numerical implementation for periodic and transient elastic as well as inelastic transient dynamic analyses of two-dimensional, axisymmetric and three-dimensional solids are presented. The inelastic formulation is based on an initial stress approach and is the first of its kind in the field of Boundary Element Methods. This formulation employs the Navier-Cauchy equation of motion, Graffi's dynamic reciprocal theorem, Stokes' fundamental solution, and the divergence theorem, together with kinematical and constitutive equations to obtain the pertinent integral equations of the problem in the time domain within the context of the small displacement theory of elastoplasticity. The dynamic (periodic, transient as well as nonlinear transient) formulations have been applied to a range of problems. The numerical formulations presented here are included in the BEST3D and GPBEST systems.

Numerical simulation of a turbulent flame stabilized behind a rearward-facing step

NASA Technical Reports Server (NTRS)

Hsiao, C. C.; Oppenheim, A. K.; Chorin, A. J.; Ghoniem, A. F.

1985-01-01

Flow of combustible mixtures in a plane channel past a smooth contraction followed by an abrupt expansion, in a typical dump combustor configuration, is modeled by a two-dimensional numerical technique based on the random vortex method. Both the inert and the reacting case are considered. In the latter, the flame is treated as an interface, self-advancing at a prescribed normal burning speed, while the dynamic effects of expansion due to the exothermicity of combustion are expressed by volumetric source lines delineated by its front. Solutions are shown to be in satisfactory agreement with experimental results, especially with respect to global properties such as the average velocity profiles and the reattachment length. The stochastic turbulent velocity components manifest interesting differences, especially near the walls where three-dimensional effects of turbulence are expected to be of importance.

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

Tipireddy, R.; Stinis, P.; Tartakovsky, A. M.

In this paper, we present a novel approach for solving steady-state stochastic partial differential equations (PDEs) with high-dimensional random parameter space. The proposed approach combines spatial domain decomposition with basis adaptation for each subdomain. The basis adaptation is used to address the curse of dimensionality by constructing an accurate low-dimensional representation of the stochastic PDE solution (probability density function and/or its leading statistical moments) in each subdomain. Restricting the basis adaptation to a specific subdomain affords finding a locally accurate solution. Then, the solutions from all of the subdomains are stitched together to provide a global solution. We support ourmore » construction with numerical experiments for a steady-state diffusion equation with a random spatially dependent coefficient. Lastly, our results show that highly accurate global solutions can be obtained with significantly reduced computational costs.« less

An efficient finite element method for simulation of droplet spreading on a topologically rough surface

NASA Astrophysics Data System (ADS)

Luo, Li; Wang, Xiao-Ping; Cai, Xiao-Chuan

2017-11-01

We study numerically the dynamics of a three-dimensional droplet spreading on a rough solid surface using a phase-field model consisting of the coupled Cahn-Hilliard and Navier-Stokes equations with a generalized Navier boundary condition (GNBC). An efficient finite element method on unstructured meshes is introduced to cope with the complex geometry of the solid surfaces. We extend the GNBC to surfaces with complex geometry by including its weak form along different normal and tangential directions in the finite element formulation. The semi-implicit time discretization scheme results in a decoupled system for the phase function, the velocity, and the pressure. In addition, a mass compensation algorithm is introduced to preserve the mass of the droplet. To efficiently solve the decoupled systems, we present a highly parallel solution strategy based on domain decomposition techniques. We validate the newly developed solution method through extensive numerical experiments, particularly for those phenomena that can not be achieved by two-dimensional simulations. On a surface with circular posts, we study how wettability of the rough surface depends on the geometry of the posts. The contact line motion for a droplet spreading over some periodic rough surfaces are also efficiently computed. Moreover, we study the spreading process of an impacting droplet on a microstructured surface, a qualitative agreement is achieved between the numerical and experimental results. The parallel performance suggests that the proposed solution algorithm is scalable with over 4,000 processors cores with tens of millions of unknowns.

The Space-Time Conservation Element and Solution Element Method: A New High-Resolution and Genuinely Multidimensional Paradigm for Solving Conservation Laws. 1; The Two Dimensional Time Marching Schemes

NASA Technical Reports Server (NTRS)

Chang, Sin-Chung; Wang, Xiao-Yen; Chow, Chuen-Yen

1998-01-01

A new high resolution and genuinely multidimensional numerical method for solving conservation laws is being, developed. It was designed to avoid the limitations of the traditional methods. and was built from round zero with extensive physics considerations. Nevertheless, its foundation is mathmatically simple enough that one can build from it a coherent, robust. efficient and accurate numerical framework. Two basic beliefs that set the new method apart from the established methods are at the core of its development. The first belief is that, in order to capture physics more efficiently and realistically, the modeling, focus should be placed on the original integral form of the physical conservation laws, rather than the differential form. The latter form follows from the integral form under the additional assumption that the physical solution is smooth, an assumption that is difficult to realize numerically in a region of rapid chance. such as a boundary layer or a shock. The second belief is that, with proper modeling of the integral and differential forms themselves, the resulting, numerical solution should automatically be consistent with the properties derived front the integral and differential forms, e.g., the jump conditions across a shock and the properties of characteristics. Therefore a much simpler and more robust method can be developed by not using the above derived properties explicitly.

Effective equations and the inverse cascade theory for Kolmogorov flows

NASA Technical Reports Server (NTRS)

Weinan, E.; Shu, Chi-Wang

1992-01-01

We study the two dimensional Kolmogorov flows in the limit as the forcing frequency goes to infinity. Direct numerical simulation indicates that the low frequency energy spectrum evolves to a universal kappa (exp -4) decay law. We derive effective equations governing the behavior of the large scale flow quantities. We then present numerical evidence that with smooth initial data, the solution to the effective equation develops a kappa (exp -4) type singularity at a finite time. This gives a convenient explanation for the kappa (exp -4) decay law exhibited by the original Kolmogorov flows.

Least-squares Legendre spectral element solutions to sound propagation problems.

PubMed

Lin, W H

2001-02-01

This paper presents a novel algorithm and numerical results of sound wave propagation. The method is based on a least-squares Legendre spectral element approach for spatial discretization and the Crank-Nicolson [Proc. Cambridge Philos. Soc. 43, 50-67 (1947)] and Adams-Bashforth [D. Gottlieb and S. A. Orszag, Numerical Analysis of Spectral Methods: Theory and Applications (CBMS-NSF Monograph, Siam 1977)] schemes for temporal discretization to solve the linearized acoustic field equations for sound propagation. Two types of NASA Computational Aeroacoustics (CAA) Workshop benchmark problems [ICASE/LaRC Workshop on Benchmark Problems in Computational Aeroacoustics, edited by J. C. Hardin, J. R. Ristorcelli, and C. K. W. Tam, NASA Conference Publication 3300, 1995a] are considered: a narrow Gaussian sound wave propagating in a one-dimensional space without flows, and the reflection of a two-dimensional acoustic pulse off a rigid wall in the presence of a uniform flow of Mach 0.5 in a semi-infinite space. The first problem was used to examine the numerical dispersion and dissipation characteristics of the proposed algorithm. The second problem was to demonstrate the capability of the algorithm in treating sound propagation in a flow. Comparisons were made of the computed results with analytical results and results obtained by other methods. It is shown that all results computed by the present method are in good agreement with the analytical solutions and results of the first problem agree very well with those predicted by other schemes.

General design method for 3-dimensional, potential flow fields. Part 2: Computer program DIN3D1 for simple, unbranched ducts

NASA Technical Reports Server (NTRS)

Stanitz, J. D.

1985-01-01

The general design method for three-dimensional, potential, incompressible or subsonic-compressible flow developed in part 1 of this report is applied to the design of simple, unbranched ducts. A computer program, DIN3D1, is developed and five numerical examples are presented: a nozzle, two elbows, an S-duct, and the preliminary design of a side inlet for turbomachines. The two major inputs to the program are the upstream boundary shape and the lateral velocity distribution on the duct wall. As a result of these inputs, boundary conditions are overprescribed and the problem is ill posed. However, it appears that there are degrees of compatibility between these two major inputs and that, for reasonably compatible inputs, satisfactory solutions can be obtained. By not prescribing the shape of the upstream boundary, the problem presumably becomes well posed, but it is not clear how to formulate a practical design method under this circumstance. Nor does it appear desirable, because the designer usually needs to retain control over the upstream (or downstream) boundary shape. The problem is further complicated by the fact that, unlike the two-dimensional case, and irrespective of the upstream boundary shape, some prescribed lateral velocity distributions do not have proper solutions.

Time-periodic solutions of the Benjamin-Ono equation

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

Ambrose , D.M.; Wilkening, Jon

2008-04-01

We present a spectrally accurate numerical method for finding non-trivial time-periodic solutions of non-linear partial differential equations. The method is based on minimizing a functional (of the initial condition and the period) that is positive unless the solution is periodic, in which case it is zero. We solve an adjoint PDE to compute the gradient of this functional with respect to the initial condition. We include additional terms in the functional to specify the free parameters, which, in the case of the Benjamin-Ono equation, are the mean, a spatial phase, a temporal phase and the real part of one ofmore » the Fourier modes at t = 0. We use our method to study global paths of non-trivial time-periodic solutions connecting stationary and traveling waves of the Benjamin-Ono equation. As a starting guess for each path, we compute periodic solutions of the linearized problem by solving an infinite dimensional eigenvalue problem in closed form. We then use our numerical method to continue these solutions beyond the realm of linear theory until another traveling wave is reached (or until the solution blows up). By experimentation with data fitting, we identify the analytical form of the solutions on the path connecting the one-hump stationary solution to the two-hump traveling wave. We then derive exact formulas for these solutions by explicitly solving the system of ODE's governing the evolution of solitons using the ansatz suggested by the numerical simulations.« less

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.

Constraint treatment techniques and parallel algorithms for multibody dynamic analysis. Ph.D. Thesis

NASA Technical Reports Server (NTRS)

Chiou, Jin-Chern

1990-01-01

Computational procedures for kinematic and dynamic analysis of three-dimensional multibody dynamic (MBD) systems are developed from the differential-algebraic equations (DAE's) viewpoint. Constraint violations during the time integration process are minimized and penalty constraint stabilization techniques and partitioning schemes are developed. The governing equations of motion, a two-stage staggered explicit-implicit numerical algorithm, are treated which takes advantage of a partitioned solution procedure. A robust and parallelizable integration algorithm is developed. This algorithm uses a two-stage staggered central difference algorithm to integrate the translational coordinates and the angular velocities. The angular orientations of bodies in MBD systems are then obtained by using an implicit algorithm via the kinematic relationship between Euler parameters and angular velocities. It is shown that the combination of the present solution procedures yields a computationally more accurate solution. To speed up the computational procedures, parallel implementation of the present constraint treatment techniques, the two-stage staggered explicit-implicit numerical algorithm was efficiently carried out. The DAE's and the constraint treatment techniques were transformed into arrowhead matrices to which Schur complement form was derived. By fully exploiting the sparse matrix structural analysis techniques, a parallel preconditioned conjugate gradient numerical algorithm is used to solve the systems equations written in Schur complement form. A software testbed was designed and implemented in both sequential and parallel computers. This testbed was used to demonstrate the robustness and efficiency of the constraint treatment techniques, the accuracy of the two-stage staggered explicit-implicit numerical algorithm, and the speed up of the Schur-complement-based parallel preconditioned conjugate gradient algorithm on a parallel computer.

a Speculative Study on Negative-Dimensional Potential and Wave Problems by Implicit Calculus Modeling Approach

NASA Astrophysics Data System (ADS)

Chen, Wen; Wang, Fajie

Based on the implicit calculus equation modeling approach, this paper proposes a speculative concept of the potential and wave operators on negative dimensionality. Unlike the standard partial differential equation (PDE) modeling, the implicit calculus modeling approach does not require the explicit expression of the PDE governing equation. Instead the fundamental solution of physical problem is used to implicitly define the differential operator and to implement simulation in conjunction with the appropriate boundary conditions. In this study, we conjecture an extension of the fundamental solution of the standard Laplace and Helmholtz equations to negative dimensionality. And then by using the singular boundary method, a recent boundary discretization technique, we investigate the potential and wave problems using the fundamental solution on negative dimensionality. Numerical experiments reveal that the physics behaviors on negative dimensionality may differ on positive dimensionality. This speculative study might open an unexplored territory in research.

Splash singularity for water waves.

PubMed

Castro, Angel; Córdoba, Diego; Fefferman, Charles L; Gancedo, Francisco; Gómez-Serrano, Javier

2012-01-17

We exhibit smooth initial data for the two-dimensional (2D) water-wave equation for which we prove that smoothness of the interface breaks down in finite time. Moreover, we show a stability result together with numerical evidence that there exist solutions of the 2D water-wave equation that start from a graph, turn over, and collapse in a splash singularity (self-intersecting curve in one point) in finite time.

Splash singularity for water waves

PubMed Central

Castro, Angel; Córdoba, Diego; Fefferman, Charles L.; Gancedo, Francisco; Gómez-Serrano, Javier

2012-01-01

We exhibit smooth initial data for the two-dimensional (2D) water-wave equation for which we prove that smoothness of the interface breaks down in finite time. Moreover, we show a stability result together with numerical evidence that there exist solutions of the 2D water-wave equation that start from a graph, turn over, and collapse in a splash singularity (self-intersecting curve in one point) in finite time. PMID:22219372

Series expansions of rotating two and three dimensional sound fields.

PubMed

Poletti, M A

2010-12-01

The cylindrical and spherical harmonic expansions of oscillating sound fields rotating at a constant rate are derived. These expansions are a generalized form of the stationary sound field expansions. The derivations are based on the representation of interior and exterior sound fields using the simple source approach and determination of the simple source solutions with uniform rotation. Numerical simulations of rotating sound fields are presented to verify the theory.

An adaptive method for a model of two-phase reactive flow on overlapping grids

NASA Astrophysics Data System (ADS)

Schwendeman, D. W.

2008-11-01

A two-phase model of heterogeneous explosives is handled computationally by a new numerical approach that is a modification of the standard Godunov scheme. The approach generates well-resolved and accurate solutions using adaptive mesh refinement on overlapping grids, and treats rationally the nozzling terms that render the otherwise hyperbolic model incapable of a conservative representation. The evolution and structure of detonation waves for a variety of one and two-dimensional configurations will be discussed with a focus given to problems of detonation diffraction and failure.

A robust, efficient equidistribution 2D grid generation method

NASA Astrophysics Data System (ADS)

Chacon, Luis; Delzanno, Gian Luca; Finn, John; Chung, Jeojin; Lapenta, Giovanni

2007-11-01

We present a new cell-area equidistribution method for two- dimensional grid adaptation [1]. The method is able to satisfy the equidistribution constraint to arbitrary precision while optimizing desired grid properties (such as isotropy and smoothness). The method is based on the minimization of the grid smoothness integral, constrained to producing a given positive-definite cell volume distribution. The procedure gives rise to a single, non-linear scalar equation with no free-parameters. We solve this equation numerically with the Newton-Krylov technique. The ellipticity property of the linearized scalar equation allows multigrid preconditioning techniques to be effectively used. We demonstrate a solution exists and is unique. Therefore, once the solution is found, the adapted grid cannot be folded due to the positivity of the constraint on the cell volumes. We present several challenging tests to show that our new method produces optimal grids in which the constraint is satisfied numerically to arbitrary precision. We also compare the new method to the deformation method [2] and show that our new method produces better quality grids. [1] G.L. Delzanno, L. Chac'on, J.M. Finn, Y. Chung, G. Lapenta, A new, robust equidistribution method for two-dimensional grid generation, in preparation. [2] G. Liao and D. Anderson, A new approach to grid generation, Appl. Anal. 44, 285--297 (1992).

Solutions of the two-dimensional Hubbard model: Benchmarks and results from a wide range of numerical algorithms

DOE PAGES

LeBlanc, J. P. F.; Antipov, Andrey E.; Becca, Federico; ...

2015-12-14

Numerical results for ground-state and excited-state properties (energies, double occupancies, and Matsubara-axis self-energies) of the single-orbital Hubbard model on a two-dimensional square lattice are presented, in order to provide an assessment of our ability to compute accurate results in the thermodynamic limit. Many methods are employed, including auxiliary-field quantum Monte Carlo, bare and bold-line diagrammatic Monte Carlo, method of dual fermions, density matrix embedding theory, density matrix renormalization group, dynamical cluster approximation, diffusion Monte Carlo within a fixed-node approximation, unrestricted coupled cluster theory, and multireference projected Hartree-Fock methods. Comparison of results obtained by different methods allows for the identification ofmore » uncertainties and systematic errors. The importance of extrapolation to converged thermodynamic-limit values is emphasized. Furthermore, cases where agreement between different methods is obtained establish benchmark results that may be useful in the validation of new approaches and the improvement of existing methods.« less

Three-dimensional time-dependent computer modeling of the electrothermal atomizers for analytical spectrometry

NASA Astrophysics Data System (ADS)

Tsivilskiy, I. V.; Nagulin, K. Yu.; Gilmutdinov, A. Kh.

2016-02-01

A full three-dimensional nonstationary numerical model of graphite electrothermal atomizers of various types is developed. The model is based on solution of a heat equation within solid walls of the atomizer with a radiative heat transfer and numerical solution of a full set of Navier-Stokes equations with an energy equation for a gas. Governing equations for the behavior of a discrete phase, i.e., atomic particles suspended in a gas (including gas-phase processes of evaporation and condensation), are derived from the formal equations molecular kinetics by numerical solution of the Hertz-Langmuir equation. The following atomizers test the model: a Varian standard heated electrothermal vaporizer (ETV), a Perkin Elmer standard THGA transversely heated graphite tube with integrated platform (THGA), and the original double-stage tube-helix atomizer (DSTHA). The experimental verification of computer calculations is carried out by a method of shadow spectral visualization of the spatial distributions of atomic and molecular vapors in an analytical space of an atomizer.

On the construction of a direct numerical simulation of a breaking inertia-gravity wave in the upper mesosphere

NASA Astrophysics Data System (ADS)

Fruman, Mark D.; Remmler, Sebastian; Achatz, Ulrich; Hickel, Stefan

2014-10-01

A systematic approach to the direct numerical simulation (DNS) of breaking upper mesospheric inertia-gravity waves of amplitude close to or above the threshold for static instability is presented. Normal mode or singular vector analysis applied in a frame of reference moving with the phase velocity of the wave (in which the wave is a steady solution) is used to determine the most likely scale and structure of the primary instability and to initialize nonlinear "2.5-D" simulations (with three-dimensional velocity and vorticity fields but depending only on two spatial coordinates). Singular vector analysis is then applied to the time-dependent 2.5-D solution to predict the transition of the breaking event to three-dimensional turbulence and to initialize three-dimensional DNS. The careful choice of the computational domain and the relatively low Reynolds numbers, on the order of 25,000, relevant to breaking waves in the upper mesosphere, makes the three-dimensional DNS tractable with present-day computing clusters. Three test cases are presented: a statically unstable low-frequency inertia-gravity wave, a statically and dynamically stable inertia-gravity wave, and a statically unstable high-frequency gravity wave. The three-dimensional DNS are compared to ensembles of 2.5-D simulations. In general, the decay of the wave and generation of turbulence is faster in three dimensions, but the results are otherwise qualitatively and quantitatively similar, suggesting that results of 2.5-D simulations are meaningful if the domain and initial condition are chosen properly.

Effects of ethanol addition on micellar solubilization and plume migration during surfactant enhanced recovery of tetrachloroethene.

PubMed

Taylor, Tammy P; Rathfelder, Klaus M; Pennell, Kurt D; Abriola, Linda M

2004-03-01

Alcohol addition has been suggested for use in combination with surfactant flushing to enhance solubilization kinetics and permit density control of dense non-aqueous phase liquid (DNAPL)-laden surfactant plumes. This study examined the effects of adding ethanol (EtOH) to a 4% Tween 80 (polyoxyethylene (20) sorbitan monooleate) solution used to flush tetrachloroethene (PCE)-contaminated porous media. The influence of EtOH concentration, subsurface layering and scale on flushing solution delivery and PCE recovery was investigated through a combination of experimental and mathematical modeling studies. Results of batch experiments demonstrated that the addition of 2.5%, 5% and 10% (wt.) EtOH incrementally increased the PCE solubilization capacity and viscosity of the surfactant solution, while reducing solution density from 1.002 to 0.986 g/cm3. Effluent concentration data obtained from one-dimensional (1-D) column experiments were used to characterize rate-limited micellar solubilization of residual PCE, which was strongly dependent upon flow velocity and weakly dependent upon EtOH concentration. Two-dimensional (2-D) box studies illustrated that minor differences (0.008 g/cm3) between flushing and resident solution density can strongly influence surfactant front propagation. A two-dimensional multiphase simulator, MISER, was used to model the influence of EtOH composition on the aqueous flow field and PCE mass recovery. The ability of the numerical simulator to predict effluent concentrations and front propagation was demonstrated for both 1-D columns and 2-D boxes flushed with EtOH-amended Tween 80 solutions. Results of this study quantify the potential influence of alcohol addition on surfactant solution properties and solubilization capacity, and demonstrate the importance of considering small density variations in remedial design.

The stability of freak waves with regard to external impact and perturbation of initial data

NASA Astrophysics Data System (ADS)

Smirnova, Anna; Shamin, Roman

2014-05-01

We investigate solutions of the equations, describing freak waves, in perspective of stability with regard to external impact and perturbation of initial data. The modeling of freak waves is based on numerical solution of equations describing a non-stationary potential flow of the ideal fluid with a free surface. We consider the two-dimensional infinitely deep flow. For waves modeling we use the equations in conformal variables. The variant of these equations is offered in [1]. Mathematical correctness of these equations was discussed in [2]. These works establish the uniqueness of solutions, offer the effective numerical solution calculation methods, prove the numerical convergence of these methods. The important aspect of numerical modeling of freak waves is the stability of solutions, describing these waves. In this work we study the questions of stability with regards to external impact and perturbation of initial data. We showed the stability of freak waves numerical model, corresponding to the external impact. We performed series of computational experiments with various freak wave initial data and random external impact. This impact means the power density on free surface. In each experiment examine two waves: the wave that was formed by external impact and without one. In all the experiments we see the stability of equation`s solutions. The random external impact practically does not change the time of freak wave formation and its form. Later our work progresses to the investigation of solution's stability under perturbations of initial data. We take the initial data that provide a freak wave and get the numerical solution. In common we take the numerical solution of equation with perturbation of initial data. The computing experiments showed that the freak waves equations solutions are stable under perturbations of initial data.So we can make a conclusion that freak waves are stable relatively external perturbation and perturbation of initial data both. 1. Zakharov V.E., Dyachenko A.I., Vasilyev O.A. New method for numerical simulation of a nonstationary potential flow of incompressible fluid with a free surface// Eur. J.~Mech. B Fluids. 2002. V. 21. P. 283-291. 2. R.V. Shamin. Dynamics of an Ideal Liquid with a Free Surface in Conformal Variables // Journal of Mathematical Sciences, Vol. 160, No. 5, 2009. P. 537-678. 3. R.V. Shamin, V.E. Zakharov, A.I. Dyachenko. How probability for freak wave formation can be found // THE EUROPEAN PHYSICAL JOURNAL - SPECIAL TOPICS Volume 185, Number 1, 113-124, DOI: 10.1140/epjst/e2010-01242-y

On the equivalence of dynamically orthogonal and bi-orthogonal methods: Theory and numerical simulations

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

Choi, Minseok; Sapsis, Themistoklis P.; Karniadakis, George Em, E-mail: george_karniadakis@brown.edu

2014-08-01

The Karhunen–Lòeve (KL) decomposition provides a low-dimensional representation for random fields as it is optimal in the mean square sense. Although for many stochastic systems of practical interest, described by stochastic partial differential equations (SPDEs), solutions possess this low-dimensional character, they also have a strongly time-dependent form and to this end a fixed-in-time basis may not describe the solution in an efficient way. Motivated by this limitation of standard KL expansion, Sapsis and Lermusiaux (2009) [26] developed the dynamically orthogonal (DO) field equations which allow for the simultaneous evolution of both the spatial basis where uncertainty ‘lives’ but also themore » stochastic characteristics of uncertainty. Recently, Cheng et al. (2013) [28] introduced an alternative approach, the bi-orthogonal (BO) method, which performs the exact same tasks, i.e. it evolves the spatial basis and the stochastic characteristics of uncertainty. In the current work we examine the relation of the two approaches and we prove theoretically and illustrate numerically their equivalence, in the sense that one method is an exact reformulation of the other. We show this by deriving a linear and invertible transformation matrix described by a matrix differential equation that connects the BO and the DO solutions. We also examine a pathology of the BO equations that occurs when two eigenvalues of the solution cross, resulting in an instantaneous, infinite-speed, internal rotation of the computed spatial basis. We demonstrate that despite the instantaneous duration of the singularity this has important implications on the numerical performance of the BO approach. On the other hand, it is observed that the BO is more stable in nonlinear problems involving a relatively large number of modes. Several examples, linear and nonlinear, are presented to illustrate the DO and BO methods as well as their equivalence.« less

Finite difference methods for reducing numerical diffusion in TEACH-type calculations. [Teaching Elliptic Axisymmetric Characteristics Heuristically

NASA Technical Reports Server (NTRS)

Syed, S. A.; Chiappetta, L. M.

1985-01-01

A methodological evaluation for two-finite differencing schemes for computer-aided gas turbine design is presented. The two computational schemes include; a Bounded Skewed Finite Differencing Scheme (BSUDS); and a Quadratic Upwind Differencing Scheme (QSDS). In the evaluation, the derivations of the schemes were incorporated into two-dimensional and three-dimensional versions of the Teaching Axisymmetric Characteristics Heuristically (TEACH) computer code. Assessments were made according to performance criteria for the solution of problems of turbulent, laminar, and coannular turbulent flow. The specific performance criteria used in the evaluation were simplicity, accuracy, and computational economy. It is found that the BSUDS scheme performed better with respect to the criteria than the QUDS. Some of the reasons for the more successful performance BSUDS are discussed.

A two-state hysteresis model from high-dimensional friction

PubMed Central

Biswas, Saurabh; Chatterjee, Anindya

2015-01-01

In prior work (Biswas & Chatterjee 2014 Proc. R. Soc. A 470, 20130817 (doi:10.1098/rspa.2013.0817)), we developed a six-state hysteresis model from a high-dimensional frictional system. Here, we use a more intuitively appealing frictional system that resembles one studied earlier by Iwan. The basis functions now have simple analytical description. The number of states required decreases further, from six to the theoretical minimum of two. The number of fitted parameters is reduced by an order of magnitude, to just six. An explicit and faster numerical solution method is developed. Parameter fitting to match different specified hysteresis loops is demonstrated. In summary, a new two-state model of hysteresis is presented that is ready for practical implementation. Essential Matlab code is provided. PMID:26587279

Particle paths and phase plane for time-dependent similarity solutions of the one-dimensional Vlasov-Maxwell equations

NASA Technical Reports Server (NTRS)

Roberts, Dana Aaron; Abraham-Shrauner, Barbara

1987-01-01

The phase trajectories of particles in a plasma described by the one-dimensional Vlasov-Maxwell equations are determined qualitatively, analyzing exact general similarity solutions for the cases of temporally damped and growing (sinusoidal or localized) electric fields. The results of numerical integration in both untransformed and Lie-group point-transformed coordinates are presented in extensive graphs and characterized in detail. The implications of the present analysis for the stability of BGK equilibria are explored, and the existence of nonlinear solutions arbitrarily close to and significantly different from the BGK solutions is demonstrated.

Computer program for calculating full potential transonic, quasi-three-dimensional flow through a rotating turbomachinery blade row

NASA Technical Reports Server (NTRS)

Farrell, C. A.

1982-01-01

A fast, reliable computer code is described for calculating the flow field about a cascade of arbitrary two dimensional airfoils. The method approximates the three dimensional flow in a turbomachinery blade row by correcting for stream tube convergence and radius change in the throughflow direction. A fully conservative solution of the full potential equation is combined with the finite volume technique on a body-fitted periodic mesh, with an artificial density imposed in the transonic region to insure stability and the capture of shock waves. The instructions required to set up and use the code are included. The name of the code is QSONIC. A numerical example is also given to illustrate the output of the program.

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.

Computer analysis of multicircuit shells of revolution by the field method

NASA Technical Reports Server (NTRS)

Cohen, G. A.

1975-01-01

The field method, presented previously for the solution of even-order linear boundary value problems defined on one-dimensional open branch domains, is extended to boundary value problems defined on one-dimensional domains containing circuits. This method converts the boundary value problem into two successive numerically stable initial value problems, which may be solved by standard forward integration techniques. In addition, a new method for the treatment of singular boundary conditions is presented. This method, which amounts to a partial interchange of the roles of force and displacement variables, is problem independent with respect to both accuracy and speed of execution. This method was implemented in a computer program to calculate the static response of ring stiffened orthotropic multicircuit shells of revolution to asymmetric loads. Solutions are presented for sample problems which illustrate the accuracy and efficiency of the method.

Intercomparison of 3D pore-scale flow and solute transport simulation methods

DOE PAGES

Mehmani, Yashar; Schoenherr, Martin; Pasquali, Andrea; ...

2015-09-28

Multiple numerical approaches have been developed to simulate porous media fluid flow and solute transport at the pore scale. These include 1) methods that explicitly model the three-dimensional geometry of pore spaces and 2) methods that conceptualize the pore space as a topologically consistent set of stylized pore bodies and pore throats. In previous work we validated a model of the first type, using computational fluid dynamics (CFD) codes employing a standard finite volume method (FVM), against magnetic resonance velocimetry (MRV) measurements of pore-scale velocities. Here we expand that validation to include additional models of the first type based onmore » the lattice Boltzmann method (LBM) and smoothed particle hydrodynamics (SPH), as well as a model of the second type, a pore-network model (PNM). The PNM approach used in the current study was recently improved and demonstrated to accurately simulate solute transport in a two-dimensional experiment. While the PNM approach is computationally much less demanding than direct numerical simulation methods, the effect of conceptualizing complex three-dimensional pore geometries on solute transport in the manner of PNMs has not been fully determined. We apply all four approaches (FVM-based CFD, LBM, SPH and PNM) to simulate pore-scale velocity distributions and (for capable codes) nonreactive solute transport, and intercompare the model results. Comparisons are drawn both in terms of macroscopic variables (e.g., permeability, solute breakthrough curves) and microscopic variables (e.g., local velocities and concentrations). Generally good agreement was achieved among the various approaches, but some differences were observed depending on the model context. The intercomparison work was challenging because of variable capabilities of the codes, and inspired some code enhancements to allow consistent comparison of flow and transport simulations across the full suite of methods. This paper provides support for confidence in a variety of pore-scale modeling methods and motivates further development and application of pore-scale simulation methods.« less

Intercomparison of 3D pore-scale flow and solute transport simulation methods

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

Yang, Xiaofan; Mehmani, Yashar; Perkins, William A.

2016-09-01

Multiple numerical approaches have been developed to simulate porous media fluid flow and solute transport at the pore scale. These include 1) methods that explicitly model the three-dimensional geometry of pore spaces and 2) methods that conceptualize the pore space as a topologically consistent set of stylized pore bodies and pore throats. In previous work we validated a model of the first type, using computational fluid dynamics (CFD) codes employing a standard finite volume method (FVM), against magnetic resonance velocimetry (MRV) measurements of pore-scale velocities. Here we expand that validation to include additional models of the first type based onmore » the lattice Boltzmann method (LBM) and smoothed particle hydrodynamics (SPH), as well as a model of the second type, a pore-network model (PNM). The PNM approach used in the current study was recently improved and demonstrated to accurately simulate solute transport in a two-dimensional experiment. While the PNM approach is computationally much less demanding than direct numerical simulation methods, the effect of conceptualizing complex three-dimensional pore geometries on solute transport in the manner of PNMs has not been fully determined. We apply all four approaches (FVM-based CFD, LBM, SPH and PNM) to simulate pore-scale velocity distributions and (for capable codes) nonreactive solute transport, and intercompare the model results. Comparisons are drawn both in terms of macroscopic variables (e.g., permeability, solute breakthrough curves) and microscopic variables (e.g., local velocities and concentrations). Generally good agreement was achieved among the various approaches, but some differences were observed depending on the model context. The intercomparison work was challenging because of variable capabilities of the codes, and inspired some code enhancements to allow consistent comparison of flow and transport simulations across the full suite of methods. This study provides support for confidence in a variety of pore-scale modeling methods and motivates further development and application of pore-scale simulation methods.« less

Modeling electrokinetic flows by consistent implicit incompressible smoothed particle hydrodynamics

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

Pan, Wenxiao; Kim, Kyungjoo; Perego, Mauro

2017-04-01

We present an efficient implicit incompressible smoothed particle hydrodynamics (I2SPH) discretization of Navier-Stokes, Poisson-Boltzmann, and advection-diffusion equations subject to Dirichlet or Robin boundary conditions. It is applied to model various two and three dimensional electrokinetic flows in simple or complex geometries. The I2SPH's accuracy and convergence are examined via comparison with analytical solutions, grid-based numerical solutions, or empirical models. The new method provides a framework to explore broader applications of SPH in microfluidics and complex fluids with charged objects, such as colloids and biomolecules, in arbitrary complex geometries.

Endogenous population growth may imply chaos.

PubMed

Prskawetz, A; Feichtinger, G

1995-01-01

The authors consider a discrete-time neoclassical growth model with an endogenous rate of population growth. The resulting one-dimensional map for the capital intensity has a tilted z-shape. Using the theory of nonlinear dynamical systems, they obtain numerical results on the qualitative behavior of time paths for changing parameter values. Besides stable and periodic solutions, erratic time paths may result. In particular, myopic and far-sighted economies--assumed to be characterized by low and high savings rate respectively--are characterized by stable per capita capital stocks, while solutions with chaotic windows exist between these two extremes.

Validation of a numerical method for interface-resolving simulation of multicomponent gas-liquid mass transfer and evaluation of multicomponent diffusion models

NASA Astrophysics Data System (ADS)

Woo, Mino; Wörner, Martin; Tischer, Steffen; Deutschmann, Olaf

2018-03-01

The multicomponent model and the effective diffusivity model are well established diffusion models for numerical simulation of single-phase flows consisting of several components but are seldom used for two-phase flows so far. In this paper, a specific numerical model for interfacial mass transfer by means of a continuous single-field concentration formulation is combined with the multicomponent model and effective diffusivity model and is validated for multicomponent mass transfer. For this purpose, several test cases for one-dimensional physical or reactive mass transfer of ternary mixtures are considered. The numerical results are compared with analytical or numerical solutions of the Maxell-Stefan equations and/or experimental data. The composition-dependent elements of the diffusivity matrix of the multicomponent and effective diffusivity model are found to substantially differ for non-dilute conditions. The species mole fraction or concentration profiles computed with both diffusion models are, however, for all test cases very similar and in good agreement with the analytical/numerical solutions or measurements. For practical computations, the effective diffusivity model is recommended due to its simplicity and lower computational costs.

A free energy satisfying discontinuous Galerkin method for one-dimensional Poisson-Nernst-Planck systems

NASA Astrophysics Data System (ADS)

Liu, Hailiang; Wang, Zhongming

2017-01-01

We design an arbitrary-order free energy satisfying discontinuous Galerkin (DG) method for solving time-dependent Poisson-Nernst-Planck systems. Both the semi-discrete and fully discrete DG methods are shown to satisfy the corresponding discrete free energy dissipation law for positive numerical solutions. Positivity of numerical solutions is enforced by an accuracy-preserving limiter in reference to positive cell averages. Numerical examples are presented to demonstrate the high resolution of the numerical algorithm and to illustrate the proven properties of mass conservation, free energy dissipation, as well as the preservation of steady states.

Nonnegative methods for bilinear discontinuous differencing of the S N equations on quadrilaterals

DOE PAGES

Maginot, Peter G.; Ragusa, Jean C.; Morel, Jim E.

2016-12-22

Historically, matrix lumping and ad hoc flux fixups have been the only methods used to eliminate or suppress negative angular flux solutions associated with the unlumped bilinear discontinuous (UBLD) finite element spatial discretization of the two-dimensional S N equations. Though matrix lumping inhibits negative angular flux solutions of the S N equations, it does not guarantee strictly positive solutions. In this paper, we develop and define a strictly nonnegative, nonlinear, Petrov-Galerkin finite element method that fully preserves the bilinear discontinuous spatial moments of the transport equation. Additionally, we define two ad hoc fixups that maintain particle balance and explicitly setmore » negative nodes of the UBLD finite element solution to zero but use different auxiliary equations to fully define their respective solutions. We assess the ability to inhibit negative angular flux solutions and the accuracy of every spatial discretization that we consider using a glancing void test problem with a discontinuous solution known to stress numerical methods. Though significantly more computationally intense, the nonlinear Petrov-Galerkin scheme results in a strictly nonnegative solution and is a more accurate solution than all the other methods considered. One fixup, based on shape preserving, results in a strictly nonnegative final solution but has increased numerical diffusion relative to the Petrov-Galerkin scheme and is less accurate than the UBLD solution. The second fixup, which preserves as many spatial moments as possible while setting negative values of the unlumped solution to zero, is less accurate than the Petrov-Galerkin scheme but is more accurate than the other fixup. However, it fails to guarantee a strictly nonnegative final solution. As a result, the fully lumped bilinear discontinuous finite element solution is the least accurate method, with significantly more numerical diffusion than the Petrov-Galerkin scheme and both fixups.« less

Nonnegative methods for bilinear discontinuous differencing of the S N equations on quadrilaterals

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

Maginot, Peter G.; Ragusa, Jean C.; Morel, Jim E.

Historically, matrix lumping and ad hoc flux fixups have been the only methods used to eliminate or suppress negative angular flux solutions associated with the unlumped bilinear discontinuous (UBLD) finite element spatial discretization of the two-dimensional S N equations. Though matrix lumping inhibits negative angular flux solutions of the S N equations, it does not guarantee strictly positive solutions. In this paper, we develop and define a strictly nonnegative, nonlinear, Petrov-Galerkin finite element method that fully preserves the bilinear discontinuous spatial moments of the transport equation. Additionally, we define two ad hoc fixups that maintain particle balance and explicitly setmore » negative nodes of the UBLD finite element solution to zero but use different auxiliary equations to fully define their respective solutions. We assess the ability to inhibit negative angular flux solutions and the accuracy of every spatial discretization that we consider using a glancing void test problem with a discontinuous solution known to stress numerical methods. Though significantly more computationally intense, the nonlinear Petrov-Galerkin scheme results in a strictly nonnegative solution and is a more accurate solution than all the other methods considered. One fixup, based on shape preserving, results in a strictly nonnegative final solution but has increased numerical diffusion relative to the Petrov-Galerkin scheme and is less accurate than the UBLD solution. The second fixup, which preserves as many spatial moments as possible while setting negative values of the unlumped solution to zero, is less accurate than the Petrov-Galerkin scheme but is more accurate than the other fixup. However, it fails to guarantee a strictly nonnegative final solution. As a result, the fully lumped bilinear discontinuous finite element solution is the least accurate method, with significantly more numerical diffusion than the Petrov-Galerkin scheme and both fixups.« less

The role of the dendritic growth model dimensionality in predicting the Columnar to Equiaxed Transition (CET)

NASA Astrophysics Data System (ADS)

Seredyński, M.; Rebow, M.; Banaszek, J.

2017-06-01

The dendrite tip kinetics model accuracy relies on the reliability of the stability constant used, which is usually experimentally determined for 3D situations and applied to 2D models. The paper reports authors` attempts to cure the situation by deriving 2D dendritic tip scaling parameter for aluminium-based alloy: Al-4wt%Cu. The obtained parameter is then incorporated into the KGT dendritic growth model in order to compare it with the original 3D KGT counterpart and to derive two-dimensional and three-dimensional versions of the modified Hunt's analytical model for the columnar-to-equiaxed transition (CET). The conclusions drawn from the above analysis are further confirmed through numerical calculations of the two cases of Al-4wt%Cu metallic alloy solidification using the front tracking technique. Results, including the porous zone-under-cooled liquid front position, the calculated solutal under-cooling, the average temperature gradient at a front of the dendrite tip envelope and a new predictor of the relative tendency to form an equiaxed zone, are shown, compared and discussed for two numerical cases. The necessity to calculate sufficiently precise values of the tip scaling parameter in 2D and 3D is stressed.

Aerodynamic interaction between vortical wakes and lifting two-dimensional bodies

NASA Technical Reports Server (NTRS)

Stremel, Paul M.

1989-01-01

Unsteady rotor wake interactions with the empennage, tail boom, and other aerodynamic surfaces of a helicopter have a significant influence on its aerodynamic performance, the ride quality, and vibration. A numerical method for computing the aerodynamic interaction between an interacting vortex wake and the viscous flow about arbitrary two-dimensional bodies was developed to address this helicopter problem. The method solves for the flow field velocities on a body-fitted computational mesh using finite-difference techniques. The interacting vortex wake is represented by an array of discrete vortices which, in turn, are represented by a finite-core model. The evolution of the interacting vortex wake is calculated by Lagrangian techniques. The viscous flow field of the two-dimensional body is calculated on an Eulerian grid. The flow around circular and elliptic cylinders in the absence of an interacting vortex wake was calculated. These results compare very well with other numerical results and with results obtained from experiment and thereby demonstrate the accuracy of the viscous solution. The interaction of a rotor wake with the flow about a 4 to 1 elliptic cylinder at 45 degree incidence was calculated for a Reynolds number of 3000. The results demonstrate the significant variations in the lift and drag on the elliptic cylinder in the presence of the interacting rotor wake.

Computer simulation of two-dimensional unsteady flows in estuaries and embayments by the method of characteristics : basic theory and the formulation of the numerical method

USGS Publications Warehouse

Lai, Chintu

1977-01-01

Two-dimensional unsteady flows of homogeneous density in estuaries and embayments can be described by hyperbolic, quasi-linear partial differential equations involving three dependent and three independent variables. A linear combination of these equations leads to a parametric equation of characteristic form, which consists of two parts: total differentiation along the bicharacteristics and partial differentiation in space. For its numerical solution, the specified-time-interval scheme has been used. The unknown, partial space-derivative terms can be eliminated first by suitable combinations of difference equations, converted from the corresponding differential forms and written along four selected bicharacteristics and a streamline. Other unknowns are thus made solvable from the known variables on the current time plane. The computation is carried to the second-order accuracy by using trapezoidal rule of integration. Means to handle complex boundary conditions are developed for practical application. Computer programs have been written and a mathematical model has been constructed for flow simulation. The favorable computer outputs suggest further exploration and development of model worthwhile. (Woodard-USGS)

Dynamical Connections in a Turbulent Fluid: Experiment and Simulation

NASA Astrophysics Data System (ADS)

Kageorge, Logan; Suri, Balachandra; Tithof, Jeff; Grigoriev, Roman; Schatz, Michael

2017-11-01

Embedded in the state space of a turbulent flow there exist invariant solutions to the Navier-Stokes equation called Exact Coherent Structures (ECS). Recent studies have demonstrated that the geometry of the ECS locally describes the evolution of the turbulent flow. Theory suggests that global connections may serve to guide the flow from the neighborhood of one ECS to that of another. We present here a numerical model of a Kolmogorov-like two-dimensional flow in which such connections have been calculated. Moreover, we present an experimental quasi-two-dimensional realization of this flow in which these connections prove dynamically relevant. This work is supported by the National Science Foundation (NSF CMMI 12-34436) and DARPA (HR0011-16-2-0033 subcontract to Georgia Tech).

Stable and Spectrally Accurate Schemes for the Navier-Stokes Equations

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

Jia, Jun; Liu, Jie

2011-01-01

In this paper, we present an accurate, efficient and stable numerical method for the incompressible Navier-Stokes equations (NSEs). The method is based on (1) an equivalent pressure Poisson equation formulation of the NSE with proper pressure boundary conditions, which facilitates the design of high-order and stable numerical methods, and (2) the Krylov deferred correction (KDC) accelerated method of lines transpose (mbox MoL{sup T}), which is very stable, efficient, and of arbitrary order in time. Numerical tests with known exact solutions in three dimensions show that the new method is spectrally accurate in time, and a numerical order of convergence 9more » was observed. Two-dimensional computational results of flow past a cylinder and flow in a bifurcated tube are also reported.« less

Numerical simulation of three dimensional transonic flows

NASA Technical Reports Server (NTRS)

Sahu, Jubaraj; Steger, Joseph L.

1987-01-01

The three-dimensional flow over a projectile has been computed using an implicit, approximately factored, partially flux-split algorithm. A simple composite grid scheme has been developed in which a single grid is partitioned into a series of smaller grids for applications which require an external large memory device such as the SSD of the CRAY X-MP/48, or multitasking. The accuracy and stability of the composite grid scheme has been tested by numerically simulating the flow over an ellipsoid at angle of attack and comparing the solution with a single grid solution. The flowfield over a projectile at M = 0.96 and 4 deg angle-of-attack has been computed using a fine grid, and compared with experiment.

The development of a revised version of multi-center molecular Ornstein-Zernike equation

NASA Astrophysics Data System (ADS)

Kido, Kentaro; Yokogawa, Daisuke; Sato, Hirofumi

2012-04-01

Ornstein-Zernike (OZ)-type theory is a powerful tool to obtain 3-dimensional solvent distribution around solute molecule. Recently, we proposed multi-center molecular OZ method, which is suitable for parallel computing of 3D solvation structure. The distribution function in this method consists of two components, namely reference and residue parts. Several types of the function were examined as the reference part to investigate the numerical robustness of the method. As the benchmark, the method is applied to water, benzene in aqueous solution and single-walled carbon nanotube in chloroform solution. The results indicate that fully-parallelization is achieved by utilizing the newly proposed reference functions.

Continuous surface force based lattice Boltzmann equation method for simulating thermocapillary flow

NASA Astrophysics Data System (ADS)

Zheng, Lin; Zheng, Song; Zhai, Qinglan

2016-02-01

In this paper, we extend a lattice Boltzmann equation (LBE) with continuous surface force (CSF) to simulate thermocapillary flows. The model is designed on our previous CSF LBE for athermal two phase flow, in which the interfacial tension forces and the Marangoni stresses as the results of the interface interactions between different phases are described by a conception of CSF. In this model, the sharp interfaces between different phases are separated by a narrow transition layers, and the kinetics and morphology evolution of phase separation would be characterized by an order parameter via Cahn-Hilliard equation which is solved in the frame work of LBE. The scalar convection-diffusion equation for temperature field is resolved by thermal LBE. The models are validated by thermal two layered Poiseuille flow, and two superimposed planar fluids at negligibly small Reynolds and Marangoni numbers for the thermocapillary driven convection, which have analytical solutions for the velocity and temperature. Then thermocapillary migration of two/three dimensional deformable droplet are simulated. Numerical results show that the predictions of present LBE agreed with the analytical solution/other numerical results.

A k-Omega Turbulence Model for Quasi-Three-Dimensional Turbomachinery Flows

NASA Technical Reports Server (NTRS)

Chima, Rodrick V.

1995-01-01

A two-equation k-omega turbulence model has been developed and applied to a quasi-three-dimensional viscous analysis code for blade-to-blade flows in turbomachinery. the code includes the effects of rotation, radius change, and variable stream sheet thickness. The flow equations are given and the explicit runge-Kutta solution scheme is described. the k-omega model equations are also given and the upwind implicit approximate-factorization solution scheme is described. Three cases were calculated: transitional flow over a flat plate, a transonic compressor rotor, and transonic turbine vane with heat transfer. Results were compared to theory, experimental data, and to results using the Baldwin-Lomax turbulence model. The two models compared reasonably well with the data and surprisingly well with each other. Although the k-omega model behaves well numerically and simulates effects of transition, freestream turbulence, and wall roughness, it was not decisively better than the Baldwin-Lomax model for the cases considered here.

Discontinuous dual-primal mixed finite elements for elliptic problems

NASA Technical Reports Server (NTRS)

Bottasso, Carlo L.; Micheletti, Stefano; Sacco, Riccardo

2000-01-01

We propose a novel discontinuous mixed finite element formulation for the solution of second-order elliptic problems. Fully discontinuous piecewise polynomial finite element spaces are used for the trial and test functions. The discontinuous nature of the test functions at the element interfaces allows to introduce new boundary unknowns that, on the one hand enforce the weak continuity of the trial functions, and on the other avoid the need to define a priori algorithmic fluxes as in standard discontinuous Galerkin methods. Static condensation is performed at the element level, leading to a solution procedure based on the sole interface unknowns. The resulting family of discontinuous dual-primal mixed finite element methods is presented in the one and two-dimensional cases. In the one-dimensional case, we show the equivalence of the method with implicit Runge-Kutta schemes of the collocation type exhibiting optimal behavior. Numerical experiments in one and two dimensions demonstrate the order accuracy of the new method, confirming the results of the analysis.

Origin of vertical orientation in two-dimensional metal halide perovskites and its effect on photovoltaic performance.

PubMed

Chen, Alexander Z; Shiu, Michelle; Ma, Jennifer H; Alpert, Matthew R; Zhang, Depei; Foley, Benjamin J; Smilgies, Detlef-M; Lee, Seung-Hun; Choi, Joshua J

2018-04-06

Thin films based on two-dimensional metal halide perovskites have achieved exceptional performance and stability in numerous optoelectronic device applications. Simple solution processing of the 2D perovskite provides opportunities for manufacturing devices at drastically lower cost compared to current commercial technologies. A key to high device performance is to align the 2D perovskite layers, during the solution processing, vertical to the electrodes to achieve efficient charge transport. However, it is yet to be understood how the counter-intuitive vertical orientations of 2D perovskite layers on substrates can be obtained. Here we report a formation mechanism of such vertically orientated 2D perovskite in which the nucleation and growth arise from the liquid-air interface. As a consequence, choice of substrates can be liberal from polymers to metal oxides depending on targeted application. We also demonstrate control over the degree of preferential orientation of the 2D perovskite layers and its drastic impact on device performance.

Two-dimensional fracture analysis of piezoelectric material based on the scaled boundary node method

NASA Astrophysics Data System (ADS)

Shen-Shen, Chen; Juan, Wang; Qing-Hua, Li

2016-04-01

A scaled boundary node method (SBNM) is developed for two-dimensional fracture analysis of piezoelectric material, which allows the stress and electric displacement intensity factors to be calculated directly and accurately. As a boundary-type meshless method, the SBNM employs the moving Kriging (MK) interpolation technique to an approximate unknown field in the circumferential direction and therefore only a set of scattered nodes are required to discretize the boundary. As the shape functions satisfy Kronecker delta property, no special techniques are required to impose the essential boundary conditions. In the radial direction, the SBNM seeks analytical solutions by making use of analytical techniques available to solve ordinary differential equations. Numerical examples are investigated and satisfactory solutions are obtained, which validates the accuracy and simplicity of the proposed approach. Project supported by the National Natural Science Foundation of China (Grant Nos. 11462006 and 21466012), the Foundation of Jiangxi Provincial Educational Committee, China (Grant No. KJLD14041), and the Foundation of East China Jiaotong University, China (Grant No. 09130020).

The calculation of a turbulent diffusion flame in a free shear flow with a statistical turbulence model

NASA Astrophysics Data System (ADS)

Bywater, R. J.

1980-01-01

Solutions are presented for the turbulent diffusion flame in a two-dimensional shear layer based upon a kinetic theory of turbulence (KTT). The fuel and oxidizer comprising the two streams are considered to react infinitely fast according to a one-step, irreversible kinetic mechanism. The solutions are obtained by direct numerical calculation of the transverse velocity probability density function (PDF) and the associated species distributions. The mean reactant profiles calculated from the solutions display the characteristic thick, turbulent flame zone. The phenomena result from the fact that in the context of the KTT, species react only when in the same velocity cell. This coincides with the known physical requirement that molecular mixing precedes reaction. The solutions demonstrate this behavior by showing how reactants can coexist in the mean, even when infinite reaction rates are enforced at each point (t,x,u) of velocity space.

Charged rotating black holes in Einstein-Maxwell-Chern-Simons theory with a negative cosmological constant

NASA Astrophysics Data System (ADS)

Blázquez-Salcedo, Jose Luis; Kunz, Jutta; Navarro-Lérida, Francisco; Radu, Eugen

2017-03-01

We consider rotating black hole solutions in five-dimensional Einstein-Maxwell-Chern-Simons theory with a negative cosmological constant and a generic value of the Chern-Simons coupling constant λ . Using both analytical and numerical techniques, we focus on cohomogeneity-1 configurations, with two equal-magnitude angular momenta, which approach at infinity a globally anti-de Sitter background. We find that the generic solutions share a number of basic properties with the known Cvetič, Lü, and Pope black holes which have λ =1 . New features occur as well; for example, when the Chern-Simons coupling constant exceeds a critical value, the solutions are no longer uniquely determined by their global charges. Moreover, the black holes possess radial excitations which can be labelled by the node number of the magnetic gauge potential function. Solutions with small values of λ possess other distinct features. For instance, the extremal black holes there form two disconnected branches, while not all near-horizon solutions are associated with global solutions.

A third-order computational method for numerical fluxes to guarantee nonnegative difference coefficients for advection-diffusion equations in a semi-conservative form

NASA Astrophysics Data System (ADS)

Sakai, K.; Watabe, D.; Minamidani, T.; Zhang, G. S.

2012-10-01

According to Godunov theorem for numerical calculations of advection equations, there exist no higher-order schemes with constant positive difference coefficients in a family of polynomial schemes with an accuracy exceeding the first-order. We propose a third-order computational scheme for numerical fluxes to guarantee the non-negative difference coefficients of resulting finite difference equations for advection-diffusion equations in a semi-conservative form, in which there exist two kinds of numerical fluxes at a cell surface and these two fluxes are not always coincident in non-uniform velocity fields. The present scheme is optimized so as to minimize truncation errors for the numerical fluxes while fulfilling the positivity condition of the difference coefficients which are variable depending on the local Courant number and diffusion number. The feature of the present optimized scheme consists in keeping the third-order accuracy anywhere without any numerical flux limiter. We extend the present method into multi-dimensional equations. Numerical experiments for advection-diffusion equations showed nonoscillatory solutions.

Scalable posterior approximations for large-scale Bayesian inverse problems via likelihood-informed parameter and state reduction

NASA Astrophysics Data System (ADS)

Cui, Tiangang; Marzouk, Youssef; Willcox, Karen

2016-06-01

Two major bottlenecks to the solution of large-scale Bayesian inverse problems are the scaling of posterior sampling algorithms to high-dimensional parameter spaces and the computational cost of forward model evaluations. Yet incomplete or noisy data, the state variation and parameter dependence of the forward model, and correlations in the prior collectively provide useful structure that can be exploited for dimension reduction in this setting-both in the parameter space of the inverse problem and in the state space of the forward model. To this end, we show how to jointly construct low-dimensional subspaces of the parameter space and the state space in order to accelerate the Bayesian solution of the inverse problem. As a byproduct of state dimension reduction, we also show how to identify low-dimensional subspaces of the data in problems with high-dimensional observations. These subspaces enable approximation of the posterior as a product of two factors: (i) a projection of the posterior onto a low-dimensional parameter subspace, wherein the original likelihood is replaced by an approximation involving a reduced model; and (ii) the marginal prior distribution on the high-dimensional complement of the parameter subspace. We present and compare several strategies for constructing these subspaces using only a limited number of forward and adjoint model simulations. The resulting posterior approximations can rapidly be characterized using standard sampling techniques, e.g., Markov chain Monte Carlo. Two numerical examples demonstrate the accuracy and efficiency of our approach: inversion of an integral equation in atmospheric remote sensing, where the data dimension is very high; and the inference of a heterogeneous transmissivity field in a groundwater system, which involves a partial differential equation forward model with high dimensional state and parameters.

Structural, thermodynamic, and electrical properties of polar fluids and ionic solutions on a hypersphere: Theoretical aspects

NASA Astrophysics Data System (ADS)

Caillol, J. M.

1992-01-01

We generalize previous work [J. Chem. Phys. 94, 597 (1991)] on an alternative to the Ewald method for the numerical simulations of Coulomb fluids. This new method consists in using as a simulation cell the three-dimensional surface of a four-dimensional sphere, or hypersphere. Here, we consider the case of polar fluids and electrolyte solutions. We derive all the formal expressions which are needed for numerical simulations of such systems. It includes a derivation of the multipolar interactions on a hypersphere, the expansion of the pair-correlation functions on rotational invariants, the expression of the static dielectric constant of a polar liquid, the expressions of the frequency-dependent conductivity and dielectric constant of an ionic solution, and the derivation of the Stillinger-Lovett sum rules for conductive systems.

On mixed derivatives type high dimensional multi-term fractional partial differential equations approximate solutions

NASA Astrophysics Data System (ADS)

Talib, Imran; Belgacem, Fethi Bin Muhammad; Asif, Naseer Ahmad; Khalil, Hammad

2017-01-01

In this research article, we derive and analyze an efficient spectral method based on the operational matrices of three dimensional orthogonal Jacobi polynomials to solve numerically the mixed partial derivatives type multi-terms high dimensions generalized class of fractional order partial differential equations. We transform the considered fractional order problem to an easily solvable algebraic equations with the aid of the operational matrices. Being easily solvable, the associated algebraic system leads to finding the solution of the problem. Some test problems are considered to confirm the accuracy and validity of the proposed numerical method. The convergence of the method is ensured by comparing our Matlab software simulations based obtained results with the exact solutions in the literature, yielding negligible errors. Moreover, comparative results discussed in the literature are extended and improved in this study.

Numerical modeling of the exterior-to-interior transmission of impulsive sound through three-dimensional, thin-walled elastic structures

NASA Astrophysics Data System (ADS)

Remillieux, Marcel C.; Pasareanu, Stephanie M.; Svensson, U. Peter

2013-12-01

Exterior propagation of impulsive sound and its transmission through three-dimensional, thin-walled elastic structures, into enclosed cavities, are investigated numerically in the framework of linear dynamics. A model was developed in the time domain by combining two numerical tools: (i) exterior sound propagation and induced structural loading are computed using the image-source method for the reflected field (specular reflections) combined with an extension of the Biot-Tolstoy-Medwin method for the diffracted field, (ii) the fully coupled vibro-acoustic response of the interior fluid-structure system is computed using a truncated modal-decomposition approach. In the model for exterior sound propagation, it is assumed that all surfaces are acoustically rigid. Since coupling between the structure and the exterior fluid is not enforced, the model is applicable to the case of a light exterior fluid and arbitrary interior fluid(s). The structural modes are computed with the finite-element method using shell elements. Acoustic modes are computed analytically assuming acoustically rigid boundaries and rectangular geometries of the enclosed cavities. This model is verified against finite-element solutions for the cases of rectangular structures containing one and two cavities, respectively.

The simulation of a two-dimensional (2D) transport problem in a rectangular region with Lattice Boltzmann method with two-relaxation-time

NASA Astrophysics Data System (ADS)

Sugiyanto, S.; Hardyanto, W.; Marwoto, P.

2018-03-01

Transport phenomena are found in many problems in many engineering and industrial sectors. We analyzed a Lattice Boltzmann method with Two-Relaxation Time (LTRT) collision operators for simulation of pollutant moving through the medium as a two-dimensional (2D) transport problem in a rectangular region model. This model consists of a 2D rectangular region with 54 length (x), 27 width (y), and it has isotropic homogeneous medium. Initially, the concentration is zero and is distributed evenly throughout the region of interest. A concentration of 1 is maintained at 9 < y < 18, whereas the concentration of zero is maintained at 0 < y < 9 and 18 < y < 27. A specific discharge (Darcy velocity) of 1.006 is assumed. A diffusion coefficient of 0.8333 is distributed uniformly with a uniform porosity of 0.35. A computer program is written in MATLAB to compute the concentration of pollutant at any specified place and time. The program shows that LTRT solution with quadratic equilibrium distribution functions (EDFs) and relaxation time τa=1.0 are in good agreement result with other numerical solutions methods such as 3DLEWASTE (Hybrid Three-dimensional Lagrangian-Eulerian Finite Element Model of Waste Transport Through Saturated-Unsaturated Media) obtained by Yeh and 3DFEMWATER-LHS (Three-dimensional Finite Element Model of Water Flow Through Saturated-Unsaturated Media with Latin Hypercube Sampling) obtained by Hardyanto.

Bifurcation analysis and phase diagram of a spin-string model with buckled states.

PubMed

Ruiz-Garcia, M; Bonilla, L L; Prados, A

2017-12-01

We analyze a one-dimensional spin-string model, in which string oscillators are linearly coupled to their two nearest neighbors and to Ising spins representing internal degrees of freedom. String-spin coupling induces a long-range ferromagnetic interaction among spins that competes with a spin-spin antiferromagnetic coupling. As a consequence, the complex phase diagram of the system exhibits different flat rippled and buckled states, with first or second order transition lines between states. This complexity translates to the two-dimensional version of the model, whose numerical solution has been recently used to explain qualitatively the rippled to buckled transition observed in scanning tunneling microscopy experiments with suspended graphene sheets. Here we describe in detail the phase diagram of the simpler one-dimensional model and phase stability using bifurcation theory. This gives additional insight into the physical mechanisms underlying the different phases and the behavior observed in experiments.

Bifurcation analysis and phase diagram of a spin-string model with buckled states

NASA Astrophysics Data System (ADS)

Ruiz-Garcia, M.; Bonilla, L. L.; Prados, A.

2017-12-01

We analyze a one-dimensional spin-string model, in which string oscillators are linearly coupled to their two nearest neighbors and to Ising spins representing internal degrees of freedom. String-spin coupling induces a long-range ferromagnetic interaction among spins that competes with a spin-spin antiferromagnetic coupling. As a consequence, the complex phase diagram of the system exhibits different flat rippled and buckled states, with first or second order transition lines between states. This complexity translates to the two-dimensional version of the model, whose numerical solution has been recently used to explain qualitatively the rippled to buckled transition observed in scanning tunneling microscopy experiments with suspended graphene sheets. Here we describe in detail the phase diagram of the simpler one-dimensional model and phase stability using bifurcation theory. This gives additional insight into the physical mechanisms underlying the different phases and the behavior observed in experiments.

On the Sensitivity of Atmospheric Ensembles to Cloud Microphysics in Long-Term Cloud-Resolving Model Simulations

NASA Technical Reports Server (NTRS)

Zeng, Xiping; Tao, Wei-Kuo; Lang, Stephen; Hou, Arthur Y.; Zhang, Minghua; Simpson, Joanne

2008-01-01

Month-long large-scale forcing data from two field campaigns are used to drive a cloud-resolving model (CRM) and produce ensemble simulations of clouds and precipitation. Observational data are then used to evaluate the model results. To improve the model results, a new parameterization of the Bergeron process is proposed that incorporates the number concentration of ice nuclei (IN). Numerical simulations reveal that atmospheric ensembles are sensitive to IN concentration and ice crystal multiplication. Two- (2D) and three-dimensional (3D) simulations are carried out to address the sensitivity of atmospheric ensembles to model dimensionality. It is found that the ensembles with high IN concentration are more sensitive to dimensionality than those with low IN concentration. Both the analytic solutions of linear dry models and the CRM output show that there are more convective cores with stronger updrafts in 3D simulations than in 2D, which explains the differing sensitivity of the ensembles to dimensionality at different IN concentrations.

Computational technique and performance of Transient Inundation Model for Rivers--2 Dimensional (TRIM2RD) : a depth-averaged two-dimensional flow model

USGS Publications Warehouse

Fulford, Janice M.

2003-01-01

A numerical computer model, Transient Inundation Model for Rivers -- 2 Dimensional (TrimR2D), that solves the two-dimensional depth-averaged flow equations is documented and discussed. The model uses a semi-implicit, semi-Lagrangian finite-difference method. It is a variant of the Trim model and has been used successfully in estuarine environments such as San Francisco Bay. The abilities of the model are documented for three scenarios: uniform depth flows, laboratory dam-break flows, and large-scale riverine flows. The model can start computations from a ?dry? bed and converge to accurate solutions. Inflows are expressed as source terms, which limits the use of the model to sufficiently long reaches where the flow reaches equilibrium with the channel. The data sets used by the investigation demonstrate that the model accurately propagates flood waves through long river reaches and simulates dam breaks with abrupt water-surface changes.

Mechanisms behind injecting the combination of nano-clay particles and polymer solution for enhanced oil recovery

NASA Astrophysics Data System (ADS)

Khalili Nezhad, Seyyed Shahram; Cheraghian, Goshtasp

2016-08-01

Laboratory investigations and field applications have proved injection of polymer solution to be an effective means to improve oil recovery for reservoirs of medium oil viscosity. The incremental oil produced in this case is the result of an increase in areal and vertical sweep efficiencies. Biopolymers and synthetic polymers are the major categories used in the petroleum industry for specific reasons. Biopolymers like xanthan are limited in their application as they are more susceptible to biodegradation. Synthetic polymers like Hydrolyzed PolyAcrylaMide (HPAM) have a much wider application as they are less susceptible to biodegradation. Furthermore, development of nanotechnology has successfully provided technical and economical viable alternatives for present materials. The objective of this study is to investigate the effect of combining clay nanoparticles with polymer solution on oil recovery. This paper includes a history match of both one-dimensional and two-dimensional polymer floods using a three-dimensional numerical model for fluid flow and mass transport. Results indicated that the amount of polymer adsorption decreased when clay nanoparticles were added to the PolyAcrylaMide solution; however, mobility ratio improvement is believed to be the main contributor for the proposed method in order to enhance much oil recovery compared to xanthan flood and HPAM flood.

The three-dimensional compressible flow in a radial inflow turbine scroll

NASA Technical Reports Server (NTRS)

Hamed, A.; Tabakoff, W.; Malak, M.

1984-01-01

This work presents the results of an analytical study and an experimental investigation of the three-dimensional flow in a turbine scroll. The finite element method is used in the iterative numerical solution of the locally linearized governing equations for the three-dimensional velocity potential field. The results of the numerical computations are compared with the experimental measurements in the scroll cross sections, which were obtained using laser Doppler velocimetry and hot wire techniques. The results of the computations show a variation in the flow conditions around the rotor periphery which was found to depend on the scroll geometry.

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

Stubos, A.K.; Caseiras, C.P.; Buchlin, J.M.

The transient two-phase flow and phase change heat transfer processes in porous media are investigated. Based on an enthalpic approach, a one-domain formulation of the problem is developed, avoiding explicit internal boundary tracking between single- and two-phase regions. An efficient numerical scheme is applied to obtain the solution on a fixed two-dimensional grid. The transient response of a liquid-saturated, self-heated porous bed is examined in detail. A physical interpretation of a liquid-saturated, self-heated porous bed is examined in detail. A physical interpretation of the computed response to fast power transients is attempted. Comparisons with experimental data are made regarding themore » average void fraction and the limiting dryout heat flux. The numerical approach is extended, keeping the one-domain formulation, to include the surrounding wall structure in the calculation.« less

Pattern formation for NO+N H3 on Pt(100): Two-dimensional numerical results

NASA Astrophysics Data System (ADS)

Uecker, Hannes

2005-01-01

The Lombardo-Fink-Imbihl model of the NO+NH3 reaction on a Pt(100) surface consists of seven coupled ordinary differential equations (ODE) and shows stable relaxation oscillations with sharp transitions in the relevant temperature range. Here we study numerically the effect of coupling of these oscillators by surface diffusion in two dimensions. We find different types of patterns, in particular phase clusters and standing waves. In models of related surface reactions such clustered solutions are known to exist only under a global coupling through the gas phase. This global coupling is replaced here by relatively fast diffusion of two variables which are kinetically slaved in the ODE. We also compare our simulations with experimental results and discuss some shortcomings of the model.

Solution of elliptic PDEs by fast Poisson solvers using a local relaxation factor

NASA Technical Reports Server (NTRS)

Chang, Sin-Chung

1986-01-01

A large class of two- and three-dimensional, nonseparable elliptic partial differential equations (PDEs) is presently solved by means of novel one-step (D'Yakanov-Gunn) and two-step (accelerated one-step) iterative procedures, using a local, discrete Fourier analysis. In addition to being easily implemented and applicable to a variety of boundary conditions, these procedures are found to be computationally efficient on the basis of the results of numerical comparison with other established methods, which lack the present one's: (1) insensitivity to grid cell size and aspect ratio, and (2) ease of convergence rate estimation by means of the coefficient of the PDE being solved. The two-step procedure is numerically demonstrated to outperform the one-step procedure in the case of PDEs with variable coefficients.

Solitary wave solutions of two-dimensional nonlinear Kadomtsev-Petviashvili dynamic equation in dust-acoustic plasmas

NASA Astrophysics Data System (ADS)

Seadawy, Aly R.

2017-09-01

Nonlinear two-dimensional Kadomtsev-Petviashvili (KP) equation governs the behaviour of nonlinear waves in dusty plasmas with variable dust charge and two temperature ions. By using the reductive perturbation method, the two-dimensional dust-acoustic solitary waves (DASWs) in unmagnetized cold plasma consisting of dust fluid, ions and electrons lead to a KP equation. We derived the solitary travelling wave solutions of the two-dimensional nonlinear KP equation by implementing sech-tanh, sinh-cosh, extended direct algebraic and fraction direct algebraic methods. We found the electrostatic field potential and electric field in the form travelling wave solutions for two-dimensional nonlinear KP equation. The solutions for the KP equation obtained by using these methods can be demonstrated precisely and efficiency. As an illustration, we used the readymade package of Mathematica program 10.1 to solve the original problem. These solutions are in good agreement with the analytical one.

Nonlinear feedback in a six-dimensional Lorenz Model: impact of an additional heating term

NASA Astrophysics Data System (ADS)

Shen, B.-W.

2015-03-01

In this study, a six-dimensional Lorenz model (6DLM) is derived, based on a recent study using a five-dimensional (5-D) Lorenz model (LM), in order to examine the impact of an additional mode and its accompanying heating term on solution stability. The new mode added to improve the representation of the steamfunction is referred to as a secondary streamfunction mode, while the two additional modes, that appear in both the 6DLM and 5DLM but not in the original LM, are referred to as secondary temperature modes. Two energy conservation relationships of the 6DLM are first derived in the dissipationless limit. The impact of three additional modes on solution stability is examined by comparing numerical solutions and ensemble Lyapunov exponents of the 6DLM and 5DLM as well as the original LM. For the onset of chaos, the critical value of the normalized Rayleigh number (rc) is determined to be 41.1. The critical value is larger than that in the 3DLM (rc ~ 24.74), but slightly smaller than the one in the 5DLM (rc ~ 42.9). A stability analysis and numerical experiments obtained using generalized LMs, with or without simplifications, suggest the following: (1) negative nonlinear feedback in association with the secondary temperature modes, as first identified using the 5DLM, plays a dominant role in providing feedback for improving the solution's stability of the 6DLM, (2) the additional heating term in association with the secondary streamfunction mode may destabilize the solution, and (3) overall feedback due to the secondary streamfunction mode is much smaller than the feedback due to the secondary temperature modes; therefore, the critical Rayleigh number of the 6DLM is comparable to that of the 5DLM. The 5DLM and 6DLM collectively suggest different roles for small-scale processes (i.e., stabilization vs. destabilization), consistent with the following statement by Lorenz (1972): If the flap of a butterfly's wings can be instrumental in generating a tornado, it can equally well be instrumental in preventing a tornado. The implications of this and previous work, as well as future work, are also discussed.

Nonlinear feedback in a six-dimensional Lorenz model: impact of an additional heating term

NASA Astrophysics Data System (ADS)

Shen, B.-W.

2015-12-01

In this study, a six-dimensional Lorenz model (6DLM) is derived, based on a recent study using a five-dimensional (5-D) Lorenz model (LM), in order to examine the impact of an additional mode and its accompanying heating term on solution stability. The new mode added to improve the representation of the streamfunction is referred to as a secondary streamfunction mode, while the two additional modes, which appear in both the 6DLM and 5DLM but not in the original LM, are referred to as secondary temperature modes. Two energy conservation relationships of the 6DLM are first derived in the dissipationless limit. The impact of three additional modes on solution stability is examined by comparing numerical solutions and ensemble Lyapunov exponents of the 6DLM and 5DLM as well as the original LM. For the onset of chaos, the critical value of the normalized Rayleigh number (rc) is determined to be 41.1. The critical value is larger than that in the 3DLM (rc ~ 24.74), but slightly smaller than the one in the 5DLM (rc ~ 42.9). A stability analysis and numerical experiments obtained using generalized LMs, with or without simplifications, suggest the following: (1) negative nonlinear feedback in association with the secondary temperature modes, as first identified using the 5DLM, plays a dominant role in providing feedback for improving the solution's stability of the 6DLM, (2) the additional heating term in association with the secondary streamfunction mode may destabilize the solution, and (3) overall feedback due to the secondary streamfunction mode is much smaller than the feedback due to the secondary temperature modes; therefore, the critical Rayleigh number of the 6DLM is comparable to that of the 5DLM. The 5DLM and 6DLM collectively suggest different roles for small-scale processes (i.e., stabilization vs. destabilization), consistent with the following statement by Lorenz (1972): "If the flap of a butterfly's wings can be instrumental in generating a tornado, it can equally well be instrumental in preventing a tornado." The implications of this and previous work, as well as future work, are also discussed.

Nonlocal electron transport: direct and Greens function solution and comparison of our model with the SNB model

NASA Astrophysics Data System (ADS)

Colombant, Denis; Manheimer, Wallace; Schmitt, Andrew J.

2013-10-01

At least two models, ours and SNB (Schurtz-Nicolai-Busquet), and two methods of solution, direct numerical solution (DS) and Greens function (GF) are being used in multi-dimensional radiation hydrodynamics codes. We present results of a laser target implosion using both methods of solution. Although our model and SNB differ in some physical content, direct comparisons have been non-existent up to now. However a paper by Marocchino et al. has recently presented the results of two nanosecond-time-scale test problems, showing that the preheat calculated by the two models are different by about three orders of magnitude. We have rerun these problems and we find much less difference between the two than they do. One can show analytically that the results should be quite similar and are about an order of magnitude less than the maximum, and two orders of magnitude more than the minimum preheating in. We have been able to trace the somewhat different results back to the different physical assumptions made in each model. Work supported by DoE-NNSA and ONR.

A two-column formalism for time-dependent modelling of stellar convection. I. Description of the method

NASA Astrophysics Data System (ADS)

Stökl, A.

2008-11-01

Context: In spite of all the advances in multi-dimensional hydrodynamics, investigations of stellar evolution and stellar pulsations still depend on one-dimensional computations. This paper devises an alternative to the mixing-length theory or turbulence models usually adopted in modelling convective transport in such studies. Aims: The present work attempts to develop a time-dependent description of convection, which reflects the essential physics of convection and that is only moderately dependent on numerical parameters and far less time consuming than existing multi-dimensional hydrodynamics computations. Methods: Assuming that the most extensive convective patterns generate the majority of convective transport, the convective velocity field is described using two parallel, radial columns to represent up- and downstream flows. Horizontal exchange, in the form of fluid flow and radiation, over their connecting interface couples the two columns and allows a simple circulating motion. The main parameters of this convective description have straightforward geometrical meanings, namely the diameter of the columns (corresponding to the size of the convective cells) and the ratio of the cross-section between up- and downdrafts. For this geometrical setup, the time-dependent solution of the equations of radiation hydrodynamics is computed from an implicit scheme that has the advantage of being unaffected by the Courant-Friedrichs-Lewy time-step limit. This implementation is part of the TAPIR-Code (short for The adaptive, implicit RHD-Code). Results: To demonstrate the approach, results for convection zones in Cepheids are presented. The convective energy transport and convective velocities agree with expectations for Cepheids and the scheme reproduces both the kinetic energy flux and convective overshoot. A study of the parameter influence shows that the type of solution derived for these stars is in fact fairly robust with respect to the constitutive numerical parameters.

Analysis of the fluid flow and heat transfer in a thin liquid film in the presence and absence of gravity

NASA Technical Reports Server (NTRS)

Rahman, M. M.; Hankey, W. L.; Faghri, A.

1991-01-01

The hydrodynamic and thermal behavior of a thin liquid film flowing over a solid horizontal surface is analyzed for both plane and radially spreading flows. The situations where the gravitational force is completely absent and where it is significant are analyzed separately and their practical relevance to a micro-gravity environment is discussed. In the presence of gravity, in addition to Reynolds number, the Froude number of the film is found to be an important parameter that determines the supercritical and subcritical flow regimes and any associated hydraulic jump. A closed-form solution is possible under some flow situations, whereas others require numerical integration of ordinary differential equations. The approximate analytical results are found to compare well with the available two-dimensional numerical solutions.

Numerical scheme approximating solution and parameters in a beam equation

NASA Astrophysics Data System (ADS)

Ferdinand, Robert R.

2003-12-01

We present a mathematical model which describes vibration in a metallic beam about its equilibrium position. This model takes the form of a nonlinear second-order (in time) and fourth-order (in space) partial differential equation with boundary and initial conditions. A finite-element Galerkin approximation scheme is used to estimate model solution. Infinite-dimensional model parameters are then estimated numerically using an inverse method procedure which involves the minimization of a least-squares cost functional. Numerical results are presented and future work to be done is discussed.

Current status of one- and two-dimensional numerical models: Successes and limitations

NASA Technical Reports Server (NTRS)

Schwartz, R. J.; Gray, J. L.; Lundstrom, M. S.

1985-01-01

The capabilities of one and two-dimensional numerical solar cell modeling programs (SCAP1D and SCAP2D) are described. The occasions when a two-dimensional model is required are discussed. The application of the models to design, analysis, and prediction are presented along with a discussion of problem areas for solar cell modeling.

Similarity solutions of some two-space-dimensional nonlinear wave evolution equations

NASA Technical Reports Server (NTRS)

Redekopp, L. G.

1980-01-01

Similarity reductions of the two-space-dimensional versions of the Korteweg-de Vries, modified Korteweg-de Vries, Benjamin-Davis-Ono, and nonlinear Schroedinger equations are presented, and some solutions of the reduced equations are discussed. Exact dispersive solutions of the two-dimensional Korteweg-de Vries equation are obtained, and the similarity solution of this equation is shown to be reducible to the second Painleve transcendent.

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.

Modeling variably saturated subsurface solute transport with MODFLOW-UZF and MT3DMS

USGS Publications Warehouse

Morway, Eric D.; Niswonger, Richard G.; Langevin, Christian D.; Bailey, Ryan T.; Healy, Richard W.

2013-01-01

The MT3DMS groundwater solute transport model was modified to simulate solute transport in the unsaturated zone by incorporating the unsaturated-zone flow (UZF1) package developed for MODFLOW. The modified MT3DMS code uses a volume-averaged approach in which Lagrangian-based UZF1 fluid fluxes and storage changes are mapped onto a fixed grid. Referred to as UZF-MT3DMS, the linked model was tested against published benchmarks solved analytically as well as against other published codes, most frequently the U.S. Geological Survey's Variably-Saturated Two-Dimensional Flow and Transport Model. Results from a suite of test cases demonstrate that the modified code accurately simulates solute advection, dispersion, and reaction in the unsaturated zone. Two- and three-dimensional simulations also were investigated to ensure unsaturated-saturated zone interaction was simulated correctly. Because the UZF1 solution is analytical, large-scale flow and transport investigations can be performed free from the computational and data burdens required by numerical solutions to Richards' equation. Results demonstrate that significant simulation runtime savings can be achieved with UZF-MT3DMS, an important development when hundreds or thousands of model runs are required during parameter estimation and uncertainty analysis. Three-dimensional variably saturated flow and transport simulations revealed UZF-MT3DMS to have runtimes that are less than one tenth of the time required by models that rely on Richards' equation. Given its accuracy and efficiency, and the wide-spread use of both MODFLOW and MT3DMS, the added capability of unsaturated-zone transport in this familiar modeling framework stands to benefit a broad user-ship.

Modeling variably saturated subsurface solute transport with MODFLOW-UZF and MT3DMS.

PubMed

Morway, Eric D; Niswonger, Richard G; Langevin, Christian D; Bailey, Ryan T; Healy, Richard W

2013-03-01

The MT3DMS groundwater solute transport model was modified to simulate solute transport in the unsaturated zone by incorporating the unsaturated-zone flow (UZF1) package developed for MODFLOW. The modified MT3DMS code uses a volume-averaged approach in which Lagrangian-based UZF1 fluid fluxes and storage changes are mapped onto a fixed grid. Referred to as UZF-MT3DMS, the linked model was tested against published benchmarks solved analytically as well as against other published codes, most frequently the U.S. Geological Survey's Variably-Saturated Two-Dimensional Flow and Transport Model. Results from a suite of test cases demonstrate that the modified code accurately simulates solute advection, dispersion, and reaction in the unsaturated zone. Two- and three-dimensional simulations also were investigated to ensure unsaturated-saturated zone interaction was simulated correctly. Because the UZF1 solution is analytical, large-scale flow and transport investigations can be performed free from the computational and data burdens required by numerical solutions to Richards' equation. Results demonstrate that significant simulation runtime savings can be achieved with UZF-MT3DMS, an important development when hundreds or thousands of model runs are required during parameter estimation and uncertainty analysis. Three-dimensional variably saturated flow and transport simulations revealed UZF-MT3DMS to have runtimes that are less than one tenth of the time required by models that rely on Richards' equation. Given its accuracy and efficiency, and the wide-spread use of both MODFLOW and MT3DMS, the added capability of unsaturated-zone transport in this familiar modeling framework stands to benefit a broad user-ship. Published 2012. This article is a U.S. Government work and is in the public domain in the USA.

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.

Numerical Simulations of Particle Deposition in Metal Foam Heat Exchangers

NASA Astrophysics Data System (ADS)

Sauret, Emilie; Saha, Suvash C.; Gu, Yuantong

2013-01-01

Australia is a high-potential country for geothermal power with reserves currently estimated in the tens of millions of petajoules, enough to power the nation for at least 1000 years at current usage. However, these resources are mainly located in isolated arid regions where water is scarce. Therefore, wet cooling systems for geothermal plants in Australia are the least attractive solution and thus air-cooled heat exchangers are preferred. In order to increase the efficiency of such heat exchangers, metal foams have been used. One issue raised by this solution is the fouling caused by dust deposition. In this case, the heat transfer characteristics of the metal foam heat exchanger can dramatically deteriorate. Exploring the particle deposition property in the metal foam exchanger becomes crucial. This paper is a numerical investigation aimed to address this issue. Two-dimensional (2D) numerical simulations of a standard one-row tube bundle wrapped with metal foam in cross-flow are performed and highlight preferential particle deposition areas.

Numerical studies of the fluid and optical fields associated with complex cavity flows

NASA Technical Reports Server (NTRS)

Atwood, Christopher A.

1992-01-01

Numerical solutions for the flowfield about several cavity configurations have been computed using the Reynolds averaged Navier-Stokes equations. Comparisons between numerical and experimental results are made in two dimensions for free shear layers and a rectangular cavity, and in three dimensions for the transonic aero-window problem of the Stratospheric Observatory for Infrared Astronomy (SOFIA). Results show that dominant acoustic frequencies and magnitudes of the self excited resonant cavity flows compare well with the experiment. In addition, solution sensitivity to artificial dissipation and grid resolution levels are determined. Optical path distortion due to the flow field is modelled geometrically and is found to match the experiment. The fluid field was computed using a diagonalized scheme within an overset mesh framework. An existing code, OVERFLOW, was utilized with the additions of characteristic boundary condition and output routines required for reduction of the unsteady data. The newly developed code is directly applicable to a generalized three dimensional structured grid zone. Details are provided in a paper included in Appendix A.

Immersed Boundary Methods for High-Resolution Simulation of Atmospheric Boundary-Layer Flow Over Complex Terrain

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

Lundquist, K A

Mesoscale models, such as the Weather Research and Forecasting (WRF) model, are increasingly used for high resolution simulations, particularly in complex terrain, but errors associated with terrain-following coordinates degrade the accuracy of the solution. Use of an alternative Cartesian gridding technique, known as an immersed boundary method (IBM), alleviates coordinate transformation errors and eliminates restrictions on terrain slope which currently limit mesoscale models to slowly varying terrain. In this dissertation, an immersed boundary method is developed for use in numerical weather prediction. Use of the method facilitates explicit resolution of complex terrain, even urban terrain, in the WRF mesoscale model.more » First, the errors that arise in the WRF model when complex terrain is present are presented. This is accomplished using a scalar advection test case, and comparing the numerical solution to the analytical solution. Results are presented for different orders of advection schemes, grid resolutions and aspect ratios, as well as various degrees of terrain slope. For comparison, results from the same simulation are presented using the IBM. Both two-dimensional and three-dimensional immersed boundary methods are then described, along with details that are specific to the implementation of IBM in the WRF code. Our IBM is capable of imposing both Dirichlet and Neumann boundary conditions. Additionally, a method for coupling atmospheric physics parameterizations at the immersed boundary is presented, making IB methods much more functional in the context of numerical weather prediction models. The two-dimensional IB method is verified through comparisons of solutions for gentle terrain slopes when using IBM and terrain-following grids. The canonical case of flow over a Witch of Agnesi hill provides validation of the basic no-slip and zero gradient boundary conditions. Specified diurnal heating in a valley, producing anabatic winds, is used to validate the use of flux (non-zero) boundary conditions. This anabatic flow set-up is further coupled to atmospheric physics parameterizations, which calculate surface fluxes, demonstrating that the IBM can be coupled to various land-surface parameterizations in atmospheric models. Additionally, the IB method is extended to three dimensions, using both trilinear and inverse distance weighted interpolations. Results are presented for geostrophic flow over a three-dimensional hill. It is found that while the IB method using trilinear interpolation works well for simple three-dimensional geometries, a more flexible and robust method is needed for extremely complex geometries, as found in three-dimensional urban environments. A second, more flexible, immersed boundary method is devised using inverse distance weighting, and results are compared to the first IBM approach. Additionally, the functionality to nest a domain with resolved complex geometry inside of a parent domain without resolved complex geometry is described. The new IBM approach is used to model urban terrain from Oklahoma City in a one-way nested configuration, where lateral boundary conditions are provided by the parent domain. Finally, the IB method is extended to include wall model parameterizations for rough surfaces. Two possible implementations are presented, one which uses the log law to reconstruct velocities exterior to the solid domain, and one which reconstructs shear stress at the immersed boundary, rather than velocity. These methods are tested on the three-dimensional canonical case of neutral atmospheric boundary layer flow over flat terrain.« less

A two-dimensional analytical model of vapor intrusion involving vertical heterogeneity.

PubMed

Yao, Yijun; Verginelli, Iason; Suuberg, Eric M

2017-05-01

In this work, we present an analytical chlorinated vapor intrusion (CVI) model that can estimate source-to-indoor air concentration attenuation by simulating two-dimensional (2-D) vapor concentration profile in vertically heterogeneous soils overlying a homogenous vapor source. The analytical solution describing the 2-D soil gas transport was obtained by applying a modified Schwarz-Christoffel mapping method. A partial field validation showed that the developed model provides results (especially in terms of indoor emission rates) in line with the measured data from a case involving a building overlying a layered soil. In further testing, it was found that the new analytical model can very closely replicate the results of three-dimensional (3-D) numerical models at steady state in scenarios involving layered soils overlying homogenous groundwater sources. By contrast, by adopting a two-layer approach (capillary fringe and vadose zone) as employed in the EPA implementation of the Johnson and Ettinger model, the spatially and temporally averaged indoor concentrations in the case of groundwater sources can be higher than the ones estimated by the numerical model up to two orders of magnitude. In short, the model proposed in this work can represent an easy-to-use tool that can simulate the subsurface soil gas concentration in layered soils overlying a homogenous vapor source while keeping the simplicity of an analytical approach that requires much less computational effort.

Three-Dimensional Simulations of Marangoni-Benard Convection in Small Containers by the Least-Squares Finite Element Method

NASA Technical Reports Server (NTRS)

Yu, Sheng-Tao; Jiang, Bo-Nan; Wu, Jie; Duh, J. C.

1996-01-01

This paper reports a numerical study of the Marangoni-Benard (MB) convection in a planar fluid layer. The least-squares finite element method (LSFEM) is employed to solve the three-dimensional Stokes equations and the energy equation. First, the governing equations are reduced to be first-order by introducing variables such as vorticity and heat fluxes. The resultant first-order system is then cast into a div-curl-grad formulation, and its ellipticity and permissible boundary conditions are readily proved. This numerical approach provides an equal-order discretization for velocity, pressure, vorticity, temperature, and heat conduction fluxes, and therefore can provide high fidelity solutions for the complex flow physics of the MB convection. Numerical results reported include the critical Marangoni numbers (M(sub ac)) for the onset of the convection in containers with various aspect ratios, and the planforms of supercritical MB flows. The numerical solutions compared favorably with the experimental results reported by Koschmieder et al..

Exact solution of the hidden Markov processes.

PubMed

Saakian, David B

2017-11-01

We write a master equation for the distributions related to hidden Markov processes (HMPs) and solve it using a functional equation. Thus the solution of HMPs is mapped exactly to the solution of the functional equation. For a general case the latter can be solved only numerically. We derive an exact expression for the entropy of HMPs. Our expression for the entropy is an alternative to the ones given before by the solution of integral equations. The exact solution is possible because actually the model can be considered as a generalized random walk on a one-dimensional strip. While we give the solution for the two second-order matrices, our solution can be easily generalized for the L values of the Markov process and M values of observables: We should be able to solve a system of L functional equations in the space of dimension M-1.

Exact solution of the hidden Markov processes

NASA Astrophysics Data System (ADS)

Saakian, David B.

2017-11-01

We write a master equation for the distributions related to hidden Markov processes (HMPs) and solve it using a functional equation. Thus the solution of HMPs is mapped exactly to the solution of the functional equation. For a general case the latter can be solved only numerically. We derive an exact expression for the entropy of HMPs. Our expression for the entropy is an alternative to the ones given before by the solution of integral equations. The exact solution is possible because actually the model can be considered as a generalized random walk on a one-dimensional strip. While we give the solution for the two second-order matrices, our solution can be easily generalized for the L values of the Markov process and M values of observables: We should be able to solve a system of L functional equations in the space of dimension M -1 .

A Fixed-point Scheme for the Numerical Construction of Magnetohydrostatic Atmospheres in Three Dimensions

NASA Astrophysics Data System (ADS)

Gilchrist, S. A.; Braun, D. C.; Barnes, G.

2016-12-01

Magnetohydrostatic models of the solar atmosphere are often based on idealized analytic solutions because the underlying equations are too difficult to solve in full generality. Numerical approaches, too, are often limited in scope and have tended to focus on the two-dimensional problem. In this article we develop a numerical method for solving the nonlinear magnetohydrostatic equations in three dimensions. Our method is a fixed-point iteration scheme that extends the method of Grad and Rubin ( Proc. 2nd Int. Conf. on Peaceful Uses of Atomic Energy 31, 190, 1958) to include a finite gravity force. We apply the method to a test case to demonstrate the method in general and our implementation in code in particular.

Two-dimensional numerical modeling and solution of convection heat transfer in turbulent He II

NASA Technical Reports Server (NTRS)

Zhang, Burt X.; Karr, Gerald R.

1991-01-01

Numerical schemes are employed to investigate heat transfer in the turbulent flow of He II. FEM is used to solve a set of equations governing the heat transfer and hydrodynamics of He II in the turbulent regime. Numerical results are compared with available experimental data and interpreted in terms of conventional heat transfer parameters such as the Prandtl number, the Peclet number, and the Nusselt number. Within the prescribed Reynolds number domain, the Gorter-Mellink thermal counterflow mechanism becomes less significant, and He II acts like an ordinary fluid. The convection heat transfer characteristics of He II in the highly turbulent regime can be successfully described by using the conventional turbulence and heat transfer theories.

Low Reynolds number numerical solutions of chaotic flow

NASA Technical Reports Server (NTRS)

Pulliam, Thomas H.

1989-01-01

Numerical computations of two-dimensional flow past an airfoil at low Mach number, large angle of attack, and low Reynolds number are reported which show a sequence of flow states leading from single-period vortex shedding to chaos via the period-doubling mechanism. Analysis of the flow in terms of phase diagrams, Poincare sections, and flowfield variables are used to substantiate these results. The critical Reynolds number for the period-doubling bifurcations is shown to be sensitive to mesh refinement and the influence of large amounts of numerical dissipation. In extreme cases, large amounts of added dissipation can delay or completely eliminate the chaotic response. The effect of artificial dissipation at these low Reynolds numbers is to produce a new effective Reynolds number for the computations.

Stochastic Least-Squares Petrov--Galerkin Method for Parameterized Linear Systems

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

Lee, Kookjin; Carlberg, Kevin; Elman, Howard C.

Here, we consider the numerical solution of parameterized linear systems where the system matrix, the solution, and the right-hand side are parameterized by a set of uncertain input parameters. We explore spectral methods in which the solutions are approximated in a chosen finite-dimensional subspace. It has been shown that the stochastic Galerkin projection technique fails to minimize any measure of the solution error. As a remedy for this, we propose a novel stochatic least-squares Petrov--Galerkin (LSPG) method. The proposed method is optimal in the sense that it produces the solution that minimizes a weightedmore » $$\\ell^2$$-norm of the residual over all solutions in a given finite-dimensional subspace. Moreover, the method can be adapted to minimize the solution error in different weighted $$\\ell^2$$-norms by simply applying a weighting function within the least-squares formulation. In addition, a goal-oriented seminorm induced by an output quantity of interest can be minimized by defining a weighting function as a linear functional of the solution. We establish optimality and error bounds for the proposed method, and extensive numerical experiments show that the weighted LSPG method outperforms other spectral methods in minimizing corresponding target weighted norms.« less

A comparison of solute-transport solution techniques and their effect on sensitivity analysis and inverse modeling results

USGS Publications Warehouse

Mehl, S.; Hill, M.C.

2001-01-01

Five common numerical techniques for solving the advection-dispersion equation (finite difference, predictor corrector, total variation diminishing, method of characteristics, and modified method of characteristics) were tested using simulations of a controlled conservative tracer-test experiment through a heterogeneous, two-dimensional sand tank. The experimental facility was constructed using discrete, randomly distributed, homogeneous blocks of five sand types. This experimental model provides an opportunity to compare the solution techniques: the heterogeneous hydraulic-conductivity distribution of known structure can be accurately represented by a numerical model, and detailed measurements can be compared with simulated concentrations and total flow through the tank. The present work uses this opportunity to investigate how three common types of results - simulated breakthrough curves, sensitivity analysis, and calibrated parameter values - change in this heterogeneous situation given the different methods of simulating solute transport. The breakthrough curves show that simulated peak concentrations, even at very fine grid spacings, varied between the techniques because of different amounts of numerical dispersion. Sensitivity-analysis results revealed: (1) a high correlation between hydraulic conductivity and porosity given the concentration and flow observations used, so that both could not be estimated; and (2) that the breakthrough curve data did not provide enough information to estimate individual values of dispersivity for the five sands. This study demonstrates that the choice of assigned dispersivity and the amount of numerical dispersion present in the solution technique influence estimated hydraulic conductivity values to a surprising degree.

Development of an upwind, finite-volume code with finite-rate chemistry

NASA Technical Reports Server (NTRS)

Molvik, Gregory A.

1994-01-01

Under this grant, two numerical algorithms were developed to predict the flow of viscous, hypersonic, chemically reacting gases over three-dimensional bodies. Both algorithms take advantage of the benefits of upwind differencing, total variation diminishing techniques, and a finite-volume framework, but obtain their solution in two separate manners. The first algorithm is a zonal, time-marching scheme, and is generally used to obtain solutions in the subsonic portions of the flow field. The second algorithm is a much less expensive, space-marching scheme and can be used for the computation of the larger, supersonic portion of the flow field. Both codes compute their interface fluxes with a temporal Riemann solver and the resulting schemes are made fully implicit including the chemical source terms and boundary conditions. Strong coupling is used between the fluid dynamic, chemical, and turbulence equations. These codes have been validated on numerous hypersonic test cases and have provided excellent comparison with existing data.

A comprehensive one-dimensional numerical model for solute transport in rivers

NASA Astrophysics Data System (ADS)

Barati Moghaddam, Maryam; Mazaheri, Mehdi; MohammadVali Samani, Jamal

2017-01-01

One of the mechanisms that greatly affect the pollutant transport in rivers, especially in mountain streams, is the effect of transient storage zones. The main effect of these zones is to retain pollutants temporarily and then release them gradually. Transient storage zones indirectly influence all phenomena related to mass transport in rivers. This paper presents the TOASTS (third-order accuracy simulation of transient storage) model to simulate 1-D pollutant transport in rivers with irregular cross-sections under unsteady flow and transient storage zones. The proposed model was verified versus some analytical solutions and a 2-D hydrodynamic model. In addition, in order to demonstrate the model applicability, two hypothetical examples were designed and four sets of well-established frequently cited tracer study data were used. These cases cover different processes governing transport, cross-section types and flow regimes. The results of the TOASTS model, in comparison with two common contaminant transport models, shows better accuracy and numerical stability.

Benchmark problems in computational aeroacoustics

NASA Technical Reports Server (NTRS)

Porter-Locklear, Freda

1994-01-01

A recent directive at NASA Langley is aimed at numerically predicting principal noise sources. During my summer stay, I worked with high-order ENO code, developed by Dr. Harold Atkins, for solving the unsteady compressible Navier-Stokes equations, as it applies to computational aeroacoustics (CAA). A CAA workshop, composed of six categories of benchmark problems, has been organized to test various numerical properties of code. My task was to determine the robustness of Atkins' code for these test problems. In one category, we tested the nonlinear wave propagation of the code for the one-dimensional Euler equations, with initial pressure, density, and velocity conditions. Using freestream boundary conditions, our results were plausible. In another category, we solved the linearized two-dimensional Euler equations to test the effectiveness of radiation boundary conditions. Here we utilized MAPLE to compute eigenvalues and eigenvectors of the Jacobian given variable and flux vectors. We experienced a minor problem with inflow and outflow boundary conditions. Next, we solved the quasi one dimensional unsteady flow equations with an incoming acoustic wave of amplitude 10(exp -6). The small amplitude sound wave was incident on a convergent-divergent nozzle. After finding a steady-state solution and then marching forward, our solution indicated that after 30 periods the acoustic wave had dissipated (a period is time required for sound wave to traverse one end of nozzle to other end).

Distributed Two-Dimensional Fourier Transforms on DSPs with an Application for Phase Retrieval

NASA Technical Reports Server (NTRS)

Smith, Jeffrey Scott

2006-01-01

Many applications of two-dimensional Fourier Transforms require fixed timing as defined by system specifications. One example is image-based wavefront sensing. The image-based approach has many benefits, yet it is a computational intensive solution for adaptive optic correction, where optical adjustments are made in real-time to correct for external (atmospheric turbulence) and internal (stability) aberrations, which cause image degradation. For phase retrieval, a type of image-based wavefront sensing, numerous two-dimensional Fast Fourier Transforms (FFTs) are used. To meet the required real-time specifications, a distributed system is needed, and thus, the 2-D FFT necessitates an all-to-all communication among the computational nodes. The 1-D floating point FFT is very efficient on a digital signal processor (DSP). For this study, several architectures and analysis of such are presented which address the all-to-all communication with DSPs. Emphasis of this research is on a 64-node cluster of Analog Devices TigerSharc TS-101 DSPs.

The point explosion with radiation transport

NASA Astrophysics Data System (ADS)

Lin, Zhiwei; Zhang, Lu; Kuang, Longyu; Jiang, Shaoen

2017-10-01

Some amount of energy is released instantaneously at the origin to generate simultaneously a spherical radiative heat wave and a spherical shock wave in the point explosion with radiation transport, which is a complicated problem due to the competition between these two waves. The point explosion problem possesses self-similar solutions when only hydrodynamic motion or only heat conduction is considered, which are Sedov solution and Barenblatt solution respectively. The point explosion problem wherein both physical mechanisms of hydrodynamic motion and heat conduction are included has been studied by P. Reinicke and A.I. Shestakov. In this talk we numerically investigate the point explosion problem wherein both physical mechanisms of hydrodynamic motion and radiation transport are taken into account. The radiation transport equation in one dimensional spherical geometry has to be solved for this problem since the ambient medium is optically thin with respect to the initially extremely high temperature at the origin. The numerical results reveal a high compression of medium and a bi-peak structure of density, which are further theoretically analyzed at the end.

A time dependent difference theory for sound propagation in ducts with flow. [characteristic of inlet and exhaust ducts of turbofan engines

NASA Technical Reports Server (NTRS)

Baumeister, K. J.

1979-01-01

A time dependent numerical solution of the linearized continuity and momentum equation was developed for sound propagation in a two dimensional straight hard or soft wall duct with a sheared mean flow. The time dependent governing acoustic difference equations and boundary conditions were developed along with a numerical determination of the maximum stable time increments. A harmonic noise source radiating into a quiescent duct was analyzed. This explicit iteration method then calculated stepwise in real time to obtain the transient as well as the steady state solution of the acoustic field. Example calculations were presented for sound propagation in hard and soft wall ducts, with no flow and plug flow. Although the problem with sheared flow was formulated and programmed, sample calculations were not examined. The time dependent finite difference analysis was found to be superior to the steady state finite difference and finite element techniques because of shorter solution times and the elimination of large matrix storage requirements.

A numerical method for computing unsteady 2-D boundary layer flows

NASA Technical Reports Server (NTRS)

Krainer, Andreas

1988-01-01

A numerical method for computing unsteady two-dimensional boundary layers in incompressible laminar and turbulent flows is described and applied to a single airfoil changing its incidence angle in time. The solution procedure adopts a first order panel method with a simple wake model to solve for the inviscid part of the flow, and an implicit finite difference method for the viscous part of the flow. Both procedures integrate in time in a step-by-step fashion, in the course of which each step involves the solution of the elliptic Laplace equation and the solution of the parabolic boundary layer equations. The Reynolds shear stress term of the boundary layer equations is modeled by an algebraic eddy viscosity closure. The location of transition is predicted by an empirical data correlation originating from Michel. Since transition and turbulence modeling are key factors in the prediction of viscous flows, their accuracy will be of dominant influence to the overall results.

Numerical methods for engine-airframe integration

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

Murthy, S.N.B.; Paynter, G.C.

1986-01-01

Various papers on numerical methods for engine-airframe integration are presented. The individual topics considered include: scientific computing environment for the 1980s, overview of prediction of complex turbulent flows, numerical solutions of the compressible Navier-Stokes equations, elements of computational engine/airframe integrations, computational requirements for efficient engine installation, application of CAE and CFD techniques to complete tactical missile design, CFD applications to engine/airframe integration, and application of a second-generation low-order panel methods to powerplant installation studies. Also addressed are: three-dimensional flow analysis of turboprop inlet and nacelle configurations, application of computational methods to the design of large turbofan engine nacelles, comparison ofmore » full potential and Euler solution algorithms for aeropropulsive flow field computations, subsonic/transonic, supersonic nozzle flows and nozzle integration, subsonic/transonic prediction capabilities for nozzle/afterbody configurations, three-dimensional viscous design methodology of supersonic inlet systems for advanced technology aircraft, and a user's technology assessment.« less

An alternative to FASTSIM for tangential solution of the wheel-rail contact

NASA Astrophysics Data System (ADS)

Sichani, Matin Sh.; Enblom, Roger; Berg, Mats

2016-06-01

In most rail vehicle dynamics simulation packages, tangential solution of the wheel-rail contact is gained by means of Kalker's FASTSIM algorithm. While 5-25% error is expected for creep force estimation, the errors of shear stress distribution, needed for wheel-rail damage analysis, may rise above 30% due to the parabolic traction bound. Therefore, a novel algorithm named FaStrip is proposed as an alternative to FASTSIM. It is based on the strip theory which extends the two-dimensional rolling contact solution to three-dimensional contacts. To form FaStrip, the original strip theory is amended to obtain accurate estimations for any contact ellipse size and it is combined by a numerical algorithm to handle spin. The comparison between the two algorithms shows that using FaStrip improves the accuracy of the estimated shear stress distribution and the creep force estimation in all studied cases. In combined lateral creepage and spin cases, for instance, the error in force estimation reduces from 18% to less than 2%. The estimation of the slip velocities in the slip zone, needed for wear analysis, is also studied. Since FaStrip is as fast as FASTSIM, it can be an alternative for tangential solution of the wheel-rail contact in simulation packages.

Distortion of liquid film discharging from twin-fluid atomizer

NASA Astrophysics Data System (ADS)

Mehring, C.; Sirignano, W. A.

2001-11-01

The nonlinear distortion and disintegration of a thin liquid film exiting from a two-dimensional twin-fluid atomizer is analyzed numerically. Pulsed gas jets impacting on both sides of the discharging liquid film at the atomizer exit generate dilational and/or sinuous deformations of the film. Both liquid phase and gas phase are inviscid and incompressible. For the liquid phase the so-called long-wavelength approximation is employed yielding a system of unsteady one-dimensional equations for the planar film. Solution of Laplace's equation for the velocity potential yields the gas-phase velocity field on both sides of the liquid stream. Coupling between both phases is described through kinematic and dynamic boundary conditions at the phase interfaces, and includes the solution of the unsteady Bernoulli equation to determine the gas-phase pressure along the interfaces. Both gas- and liquid-phase equations are solved simultaneously. Solution of Laplace's equation for the gas streams is obtained by means of a boundary-element method. Numerical solutions for the liquid phase use the Lax-Wendroff method with Richtmyer splitting. Sheet distortion resulting from the stagnation pressure of the impacting gas jets and subsequent disturbance amplification due to Kelvin-Helmholtz effects are studied for various combinations of gas-pulse timing, gas-jet impact angles, gas-to-liquid-density ratio, liquid-phase Weber number and gas-jet-to-liquid-jet-momentum ratio. Dilational and sinuous oscillations of the liquid are examined and film pinch-off is predicted.

An analysis of the vapor flow and the heat conduction through the liquid-wick and pipe wall in a heat pipe with single or multiple heat sources

NASA Technical Reports Server (NTRS)

Chen, Ming-Ming; Faghri, Amir

1990-01-01

A numerical analysis is presented for the overall performance of heat pipes with single or multiple heat sources. The analysis includes the heat conduction in the wall and liquid-wick regions as well as the compressibility effect of the vapor inside the heat pipe. The two-dimensional elliptic governing equations in conjunction with the thermodynamic equilibrium relation and appropriate boundary conditions are solved numerically. The solutions are in agreement with existing experimental data for the vapor and wall temperatures at both low and high operating temperatures.

Stochastic analysis of multiphase flow in porous media: II. Numerical simulations

NASA Astrophysics Data System (ADS)

Abin, A.; Kalurachchi, J. J.; Kemblowski, M. W.; Chang, C.-M.

1996-08-01

The first paper (Chang et al., 1995b) of this two-part series described the stochastic analysis using spectral/perturbation approach to analyze steady state two-phase (water and oil) flow in a, liquid-unsaturated, three fluid-phase porous medium. In this paper, the results between the numerical simulations and closed-form expressions obtained using the perturbation approach are compared. We present the solution to the one-dimensional, steady-state oil and water flow equations. The stochastic input processes are the spatially correlated logk where k is the intrinsic permeability and the soil retention parameter, α. These solutions are subsequently used in the numerical simulations to estimate the statistical properties of the key output processes. The comparison between the results of the perturbation analysis and numerical simulations showed a good agreement between the two methods over a wide range of logk variability with three different combinations of input stochastic processes of logk and soil parameter α. The results clearly demonstrated the importance of considering the spatial variability of key subsurface properties under a variety of physical scenarios. The variability of both capillary pressure and saturation is affected by the type of input stochastic process used to represent the spatial variability. The results also demonstrated the applicability of perturbation theory in predicting the system variability and defining effective fluid properties through the ergodic assumption.

Analytical Solution for Time-drawdown Response to Constant Pumping from a Homogeneous, Confined Horizontal Aquifer with Unidirectional Flow

NASA Astrophysics Data System (ADS)

Parrish, K. E.; Zhang, J.; Teasdale, E.

2007-12-01

An exact analytical solution to the ordinary one-dimensional partial differential equation is derived for transient groundwater flow in a homogeneous, confined, horizontal aquifer using Laplace transformation. The theoretical analysis is based on the assumption that the aquifer is homogeneous and one-dimensional (horizontal); confined between impermeable formations on top and bottom; and of infinite horizontal extent and constant thickness. It is also assumed that there is only a single pumping well penetrating the entire aquifer; flow is everywhere horizontal within the aquifer to the well; the well is pumping with a constant discharge rate; the well diameter is infinitesimally small; and the hydraulic head is uniform throughout the aquifer before pumping. Similar to the Theis solution, this solution is suited to determine transmissivity and storativity for a two- dimensional, vertically confined aquifer, such as a long vertically fractured zone of high permeability within low permeable rocks or a long, high-permeability trench inside a low-permeability porous media. In addition, it can be used to analyze time-drawdown responses to pumping and injection in similar settings. The solution can also be used to approximate the groundwater flow for unconfined conditions if (1) the variation of transmissivity is negligible (groundwater table variation is small in comparison to the saturated thickness); and (2) the unsaturated flow is negligible. The errors associated with the use of the solution to unconfined conditions depend on the accuracies of the above two assumptions. The solution can also be used to assess the impacts of recharge from a seasonal river or irrigation canal on the groundwater system by assuming uniform, time- constant recharge along the river or canal. This paper presents the details for derivation of the analytical solution. The analytical solution is compared to numerical simulation results with example cases. Its accuracy is also assessed and discussed for confined and unconfined conditions.