Mathematical problems arising in interfacial electrohydrodynamics

NASA Astrophysics Data System (ADS)

Tseluiko, Dmitri

In this work we consider the nonlinear stability of thin films in the presence of electric fields. We study a perfectly conducting thin film flow down an inclined plane in the presence of an electric field which is uniform in its undisturbed state, and normal to the plate at infinity. In addition, the effect of normal electric fields on films lying above, or hanging from, horizontal substrates is considered. Systematic asymptotic expansions are used to derive fully nonlinear long wave model equations for the scaled interface motion and corresponding flow fields. For the case of an inclined plane, higher order terms are need to be retained to regularize the problem in the sense that the long wave approximation remains valid for long times. For the case of a horizontal plane the fully nonlinear evolution equation which is derived at the leading order, is asymptotically correct and no regularization procedure is required. In both physical situations, the effect of the electric field is to introduce a non-local term which arises from the potential region above the liquid film, and enters through the electric Maxwell stresses at the interface. This term is always linearly destabilizing and produces growth rates proportional to the cubic power of the wavenumber - surface tension is included and provides a short wavelength cut-off, that is, all sufficiently short waves are linearly stable. For the case of film flow down an inclined plane, the fully nonlinear equation can produce singular solutions (for certain parameter values) after a finite time, even in the absence of an electric field. This difficulty is avoided at smaller amplitudes where the weakly nonlinear evolution is governed by an extension of the Kuramoto-Sivashinsky (KS) equation. Global existence and uniqueness results are proved, and refined estimates of the radius of the absorbing ball in L2 are obtained in terms of the parameters of the equations for a generalized class of modified KS equations. The established estimates are compared with numerical solutions of the equations which in turn suggest an optimal upper bound for the radius of the absorbing ball. A scaling argument is used to explain this, and a general conjecture is made based on extensive computations. We also carry out a complete study of the nonlinear behavior of competing physical mechanisms: long wave instability above a critical Reynolds number, short wave damping due to surface tension and intermediate growth due to the electric field. Through a combination of analysis and extensive numerical experiments, we elucidate parameter regimes that support non-uniform travelling waves, time-periodic travelling waves and complex nonlinear dynamics including chaotic interfacial oscillations. It is established that a sufficiently high electric field will drive the system to chaotic oscillations, even when the Reynolds number is smaller than the critical value below which the non-electrified problem is linearly stable. A particular case of this is Stokes flow, which is known to be stable for this class of problems (an analogous statement holds for horizontally supported films also). Our theoretical results indicate that such highly stable flows can be rendered unstable by using electric fields. This opens the way for possible heat and mass transfer applications which can benefit significantly from interfacial oscillations and interfacial turbulence. For the case of a horizontal plane, a weakly nonlinear theory is not possible due to the absence of the shear flow generated by the gravitational force along the plate when the latter is inclined. We study the fully nonlinear equation, which in this case is asymptotically correct and is obtained at the leading order. The model equation describes both overlying and hanging films - in the former case gravity is stabilizing while in the latter it is destabilizing. The numerical and theoretical analysis of the fully nonlinear evolution is complicated by the fact that the coefficients of the highest order terms (surface tension in this instance) are nonlinear. We implement a fully implicit two level numerical scheme and perform numerical experiments. We also prove global boundedness of positive periodic smooth solutions, using an appropriate energy functional. This global boundedness result is seen in all our numerical results. Through a combination of analysis and extensive numerical experiments we present evidence for global existence of positive smooth solutions. This means, in turn, that the film does not touch the wall in finite time but asymptotically at infinite time. Numerical solutions are presented to support such phenomena.

Physical uniqueness of higher-order Korteweg-de Vries theory for continuously stratified fluids without background shear

NASA Astrophysics Data System (ADS)

Shimizu, Kenji

2017-10-01

The 2nd-order Korteweg-de Vries (KdV) equation and the Gardner (or extended KdV) equation are often used to investigate internal solitary waves, commonly observed in oceans and lakes. However, application of these KdV-type equations for continuously stratified fluids to geophysical problems is hindered by nonuniqueness of the higher-order coefficients and the associated correction functions to the wave fields. This study proposes to reduce arbitrariness of the higher-order KdV theory by considering its uniqueness in the following three physical senses: (i) consistency of the nonlinear higher-order coefficients and correction functions with the corresponding phase speeds, (ii) wavenumber-independence of the vertically integrated available potential energy, and (iii) its positive definiteness. The spectral (or generalized Fourier) approach based on vertical modes in the isopycnal coordinate is shown to enable an alternative derivation of the 2nd-order KdV equation, without encountering nonuniqueness. Comparison with previous theories shows that Parseval's theorem naturally yields a unique set of special conditions for (ii) and (iii). Hydrostatic fully nonlinear solutions, derived by combining the spectral approach and simple-wave analysis, reveal that both proposed and previous 2nd-order theories satisfy (i), provided that consistent definitions are used for the wave amplitude and the nonlinear correction. This condition reduces the arbitrariness when higher-order KdV-type theories are compared with observations or numerical simulations. The coefficients and correction functions that satisfy (i)-(iii) are given by explicit formulae to 2nd order and by algebraic recurrence relationships to arbitrary order for hydrostatic fully nonlinear and linear fully nonhydrostatic effects.

Turbulence and deterministic chaos. [computational fluid dynamics

NASA Technical Reports Server (NTRS)

Deissler, Robert G.

1992-01-01

Several turbulent and nonturbulent solutions of the Navier-Stokes equations are obtained. The unaveraged equations are used numerically in conjunction with tools and concepts from nonlinear dynamics, including time series, phase portraits, Poincare sections, largest Liapunov exponents, power spectra, and strange attractors. Initially neighboring solutions for a low Reynolds number fully developed turbulence are compared. Several flows are noted: fully chaotic, complex periodic, weakly chaotic, simple periodic, and fixed-point. Of these, only fully chaotic is classified as turbulent. Besides the sustained flows, a flow which decays as it becomes turbulent is examined. For the finest grid, 128(exp 3) points, the spatial resolution appears to be quite good. As a final note, the variation of the velocity derivatives skewness of a Navier-Stokes flow as the Reynolds number goes to zero is calculated numerically. The value of the skewness is shown to become small at low Reynolds numbers, in agreement with intuitive arguments that nonlinear terms should be negligible.

Fully Nonlinear Modeling and Analysis of Precision Membranes

NASA Technical Reports Server (NTRS)

Pai, P. Frank; Young, Leyland G.

2003-01-01

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

Two-dimensional evolution equation of finite-amplitude internal gravity waves in a uniformly stratified fluid

PubMed

Kataoka; Tsutahara; Akuzawa

2000-02-14

We derive a fully nonlinear evolution equation that can describe the two-dimensional motion of finite-amplitude long internal waves in a uniformly stratified three-dimensional fluid of finite depth. The derived equation is the two-dimensional counterpart of the evolution equation obtained by Grimshaw and Yi [J. Fluid Mech. 229, 603 (1991)]. In the small-amplitude limit, our equation is reduced to the celebrated Kadomtsev-Petviashvili equation.

A Maximum Likelihood Approach for Multisample Nonlinear Structural Equation Models with Missing Continuous and Dichotomous Data

ERIC Educational Resources Information Center

Song, Xin-Yuan; Lee, Sik-Yum

2006-01-01

Structural equation models are widely appreciated in social-psychological research and other behavioral research to model relations between latent constructs and manifest variables and to control for measurement error. Most applications of SEMs are based on fully observed continuous normal data and models with a linear structural equation.…

Building Flexible User Interfaces for Solving PDEs

NASA Astrophysics Data System (ADS)

Logg, Anders; Wells, Garth N.

2010-09-01

FEniCS is a collection of software tools for the automated solution of differential equations by finite element methods. In this note, we describe how FEniCS can be used to solve a simple nonlinear model problem with varying levels of automation. At one extreme, FEniCS provides tools for the fully automated and adaptive solution of nonlinear partial differential equations. At the other extreme, FEniCS provides a range of tools that allow the computational scientist to experiment with novel solution algorithms.

Dissipative behavior of some fully non-linear KdV-type equations

NASA Astrophysics Data System (ADS)

Brenier, Yann; Levy, Doron

2000-03-01

The KdV equation can be considered as a special case of the general equation u t+f(u) x-δg(u xx) x=0, δ>0, where f is non-linear and g is linear, namely f( u)= u2/2 and g( v)= v. As the parameter δ tends to 0, the dispersive behavior of the KdV equation has been throughly investigated (see, e.g., [P.G. Drazin, Solitons, London Math. Soc. Lect. Note Ser. 85, Cambridge University Press, Cambridge, 1983; P.D. Lax, C.D. Levermore, The small dispersion limit of the Korteweg-de Vries equation, III, Commun. Pure Appl. Math. 36 (1983) 809-829; G.B. Whitham, Linear and Nonlinear Waves, Wiley/Interscience, New York, 1974] and the references therein). We show through numerical evidence that a completely different, dissipative behavior occurs when g is non-linear, namely when g is an even concave function such as g( v)=-∣ v∣ or g( v)=- v2. In particular, our numerical results hint that as δ→0 the solutions strongly converge to the unique entropy solution of the formal limit equation, in total contrast with the solutions of the KdV equation.

A discontinuous Galerkin approach for conservative modeling of fully nonlinear and weakly dispersive wave transformations

NASA Astrophysics Data System (ADS)

Sharifian, Mohammad Kazem; Kesserwani, Georges; Hassanzadeh, Yousef

2018-05-01

This work extends a robust second-order Runge-Kutta Discontinuous Galerkin (RKDG2) method to solve the fully nonlinear and weakly dispersive flows, within a scope to simultaneously address accuracy, conservativeness, cost-efficiency and practical needs. The mathematical model governing such flows is based on a variant form of the Green-Naghdi (GN) equations decomposed as a hyperbolic shallow water system with an elliptic source term. Practical features of relevance (i.e. conservative modeling over irregular terrain with wetting and drying and local slope limiting) have been restored from an RKDG2 solver to the Nonlinear Shallow Water (NSW) equations, alongside new considerations to integrate elliptic source terms (i.e. via a fourth-order local discretization of the topography) and to enable local capturing of breaking waves (i.e. via adding a detector for switching off the dispersive terms). Numerical results are presented, demonstrating the overall capability of the proposed approach in achieving realistic prediction of nearshore wave processes involving both nonlinearity and dispersion effects within a single model.

Dynamical systems theory for nonlinear evolution equations.

PubMed

Choudhuri, Amitava; Talukdar, B; Das, Umapada

2010-09-01

We observe that the fully nonlinear evolution equations of Rosenau and Hymann, often abbreviated as K(n,m) equations, can be reduced to Hamiltonian form only on a zero-energy hypersurface belonging to some potential function associated with the equations. We treat the resulting Hamiltonian equations by the dynamical systems theory and present a phase-space analysis of their stable points. The results of our study demonstrate that the equations can, in general, support both compacton and soliton solutions. For the K(2,2) and K(3,3) cases one type of solutions can be obtained from the other by continuously varying a parameter of the equations. This is not true for the K(3,2) equation for which the parameter can take only negative values. The K(2,3) equation does not have any stable point and, in the language of mechanics, represents a particle moving with constant acceleration.

Gaussian variational ansatz in the problem of anomalous sea waves: Comparison with direct numerical simulation

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

Ruban, V. P., E-mail: ruban@itp.ac.ru

2015-05-15

The nonlinear dynamics of an obliquely oriented wave packet on a sea surface is analyzed analytically and numerically for various initial parameters of the packet in relation to the problem of the so-called rogue waves. Within the Gaussian variational ansatz applied to the corresponding (1+2)-dimensional hyperbolic nonlinear Schrödinger equation (NLSE), a simplified Lagrangian system of differential equations is derived that describes the evolution of the coefficients of the real and imaginary quadratic forms appearing in the Gaussian. This model provides a semi-quantitative description of the process of nonlinear spatiotemporal focusing, which is one of the most probable mechanisms of roguemore » wave formation in random wave fields. The system of equations is integrated in quadratures, which allows one to better understand the qualitative differences between linear and nonlinear focusing regimes of a wave packet. Predictions of the Gaussian model are compared with the results of direct numerical simulation of fully nonlinear long-crested waves.« less

Solitary wave solutions and their interactions for fully nonlinear water waves with surface tension in the generalized Serre equations

NASA Astrophysics Data System (ADS)

Dutykh, Denys; Hoefer, Mark; Mitsotakis, Dimitrios

2018-04-01

Some effects of surface tension on fully nonlinear, long, surface water waves are studied by numerical means. The differences between various solitary waves and their interactions in subcritical and supercritical surface tension regimes are presented. Analytical expressions for new peaked traveling wave solutions are presented in the dispersionless case of critical surface tension. Numerical experiments are performed using a high-accurate finite element method based on smooth cubic splines and the four-stage, classical, explicit Runge-Kutta method of order 4.

New Similarity Reductions and Compacton Solutions for Boussinesq-Like Equations with Fully Nonlinear Dispersion

NASA Astrophysics Data System (ADS)

Yan, Zhen-Ya

2001-10-01

In this paper, similarity reductions of Boussinesq-like equations with nonlinear dispersion (simply called B(m,n) equations) utt=(u^n)xx+(u^m)xxxx, which is a generalized model of Boussinesq equation utt=(u^2)xx+uxxxx and modified Bousinesq equation utt=(u^3)xx+uxxxx, are considered by using the direct reduction method. As a result, several new types of similarity reductions are found. Based on the reduction equations and some simple transformations, we obtain the solitary wave solutions and compacton solutions (which are solitary waves with the property that after colliding with other compacton solutions, they re-emerge with the same coherent shape) of B(1,n) equations and B(m,m) equations, respectively. The project supported by National Key Basic Research Development Project Program of China under Grant No. G1998030600 and Doctoral Foundation of China under Grant No. 98014119

A problem in non-linear Diophantine approximation

NASA Astrophysics Data System (ADS)

Harrap, Stephen; Hussain, Mumtaz; Kristensen, Simon

2018-05-01

In this paper we obtain the Lebesgue and Hausdorff measure results for the set of vectors satisfying infinitely many fully non-linear Diophantine inequalities. The set is associated with a class of linear inhomogeneous partial differential equations whose solubility depends on a certain Diophantine condition. The failure of the Diophantine condition guarantees the existence of a smooth solution.

Nonlinear dynamics; Proceedings of the International Conference, New York, NY, December 17-21, 1979

NASA Technical Reports Server (NTRS)

Helleman, R. H. G.

1980-01-01

Papers were presented on turbulence, ergodic and integrable behavior, chaotic maps and flows, chemical and fully developed turbulence, and strange attractors. Specific attention was given to measures describing a turbulent flow, stochastization and collapse of vortex systems, a subharmonic route to turbulent convection, and weakly nonlinear turbulence in a rotating convection layer. The Korteweg-de Vries and Hill equations, plasma transport in three dimensions, a horseshoe in the dynamics of a forced beam, and the explosion of strange attractors exhibited by Duffing's equation were also considered.

MOOSE: A PARALLEL COMPUTATIONAL FRAMEWORK FOR COUPLED SYSTEMS OF NONLINEAR EQUATIONS.

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

G. Hansen; C. Newman; D. Gaston

Systems of coupled, nonlinear partial di?erential equations often arise in sim- ulation of nuclear processes. MOOSE: Multiphysics Ob ject Oriented Simulation Environment, a parallel computational framework targeted at solving these systems is presented. As opposed to traditional data / ?ow oriented com- putational frameworks, MOOSE is instead founded on mathematics based on Jacobian-free Newton Krylov (JFNK). Utilizing the mathematical structure present in JFNK, physics are modularized into “Kernels” allowing for rapid production of new simulation tools. In addition, systems are solved fully cou- pled and fully implicit employing physics based preconditioning allowing for a large amount of ?exibility even withmore » large variance in time scales. Background on the mathematics, an inspection of the structure of MOOSE and several rep- resentative solutions from applications built on the framework are presented.« less

Weighted interior penalty discretization of fully nonlinear and weakly dispersive free surface shallow water flows

NASA Astrophysics Data System (ADS)

Di Pietro, Daniele A.; Marche, Fabien

2018-02-01

In this paper, we further investigate the use of a fully discontinuous Finite Element discrete formulation for the study of shallow water free surface flows in the fully nonlinear and weakly dispersive flow regime. We consider a decoupling strategy in which we approximate the solutions of the classical shallow water equations supplemented with a source term globally accounting for the non-hydrostatic effects. This source term can be computed through the resolution of elliptic second-order linear sub-problems, which only involve second order partial derivatives in space. We then introduce an associated Symmetric Weighted Internal Penalty discrete bilinear form, allowing to deal with the discontinuous nature of the elliptic problem's coefficients in a stable and consistent way. Similar discrete formulations are also introduced for several recent optimized fully nonlinear and weakly dispersive models. These formulations are validated again several benchmarks involving h-convergence, p-convergence and comparisons with experimental data, showing optimal convergence properties.

A Zonal Approach for Prediction of Jet Noise

NASA Technical Reports Server (NTRS)

Shih, S. H.; Hixon, D. R.; Mankbadi, Reda R.

1995-01-01

A zonal approach for direct computation of sound generation and propagation from a supersonic jet is investigated. The present work splits the computational domain into a nonlinear, acoustic-source regime and a linear acoustic wave propagation regime. In the nonlinear regime, the unsteady flow is governed by the large-scale equations, which are the filtered compressible Navier-Stokes equations. In the linear acoustic regime, the sound wave propagation is described by the linearized Euler equations. Computational results are presented for a supersonic jet at M = 2. 1. It is demonstrated that no spurious modes are generated in the matching region and the computational expense is reduced substantially as opposed to fully large-scale simulation.

Extremely Fast Numerical Integration of Ocean Surface Wave Dynamics

DTIC Science & Technology

2007-09-30

sub-processor must be added as shown in the blue box of Fig. 1. We first consider the Kadomtsev - Petviashvili (KP) equation ηt + coηx +αηηx + βη ...analytic integration of the so-called “soliton equations ,” I have discovered how the GFT can be used to solved higher order equations for which study...analytical study and extremely fast numerical integration of the extended nonlinear Schroedinger equation for fully three dimensional wave motion

Recurrence due to periodic multisoliton fission in the defocusing nonlinear Schrödinger equation

NASA Astrophysics Data System (ADS)

Deng, Guo; Li, Sitai; Biondini, Gino; Trillo, Stefano

2017-11-01

We address the degree of universality of the Fermi-Pasta-Ulam recurrence induced by multisoliton fission from a harmonic excitation by analyzing the case of the semiclassical defocusing nonlinear Schrödinger equation, which models nonlinear wave propagation in a variety of physical settings. Using a suitable Wentzel-Kramers-Brillouin approach to the solution of the associated scattering problem we accurately predict, in a fully analytical way, the number and the features (amplitude and velocity) of solitonlike excitations emerging post-breaking, as a function of the dispersion smallness parameter. This also permits us to predict and analyze the near-recurrences, thereby inferring the universal character of the mechanism originally discovered for the Korteweg-deVries equation. We show, however, that important differences exist between the two models, arising from the different scaling rules obeyed by the soliton velocities.

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

Gong, Jinn-Ouk; Hwang, Jai-chan; Noh, Hyerim

We present a complete set of exact and fully non-linear equations describing all three types of cosmological perturbations—scalar, vector and tensor perturbations. We derive the equations in a thoroughly gauge-ready manner, so that any spatial and temporal gauge conditions can be employed. The equations are completely general without any physical restriction except that we assume a flat homogeneous and isotropic universe as a background. We also comment briefly on the application of our formulation to the non-expanding Minkowski background.

Building Blocks for Reliable Complex Nonlinear Numerical Simulations

NASA Technical Reports Server (NTRS)

Yee, H. C.; Mansour, Nagi N. (Technical Monitor)

2002-01-01

This talk describes some of the building blocks to ensure a higher level of confidence in the predictability and reliability (PAR) of numerical simulation of multiscale complex nonlinear problems. The focus is on relating PAR of numerical simulations with complex nonlinear phenomena of numerics. To isolate sources of numerical uncertainties, the possible discrepancy between the chosen partial differential equation (PDE) model and the real physics and/or experimental data is set aside. The discussion is restricted to how well numerical schemes can mimic the solution behavior of the underlying PDE model for finite time steps and grid spacings. The situation is complicated by the fact that the available theory for the understanding of nonlinear behavior of numerics is not at a stage to fully analyze the nonlinear Euler and Navier-Stokes equations. The discussion is based on the knowledge gained for nonlinear model problems with known analytical solutions to identify and explain the possible sources and remedies of numerical uncertainties in practical computations. Examples relevant to turbulent flow computations are included.

Building Blocks for Reliable Complex Nonlinear Numerical Simulations

NASA Technical Reports Server (NTRS)

Yee, H. C.

2005-01-01

This chapter describes some of the building blocks to ensure a higher level of confidence in the predictability and reliability (PAR) of numerical simulation of multiscale complex nonlinear problems. The focus is on relating PAR of numerical simulations with complex nonlinear phenomena of numerics. To isolate sources of numerical uncertainties, the possible discrepancy between the chosen partial differential equation (PDE) model and the real physics and/or experimental data is set aside. The discussion is restricted to how well numerical schemes can mimic the solution behavior of the underlying PDE model for finite time steps and grid spacings. The situation is complicated by the fact that the available theory for the understanding of nonlinear behavior of numerics is not at a stage to fully analyze the nonlinear Euler and Navier-Stokes equations. The discussion is based on the knowledge gained for nonlinear model problems with known analytical solutions to identify and explain the possible sources and remedies of numerical uncertainties in practical computations.

Building Blocks for Reliable Complex Nonlinear Numerical Simulations. Chapter 2

NASA Technical Reports Server (NTRS)

Yee, H. C.; Mansour, Nagi N. (Technical Monitor)

2001-01-01

This chapter describes some of the building blocks to ensure a higher level of confidence in the predictability and reliability (PAR) of numerical simulation of multiscale complex nonlinear problems. The focus is on relating PAR of numerical simulations with complex nonlinear phenomena of numerics. To isolate sources of numerical uncertainties, the possible discrepancy between the chosen partial differential equation (PDE) model and the real physics and/or experimental data is set aside. The discussion is restricted to how well numerical schemes can mimic the solution behavior of the underlying PDE model for finite time steps and grid spacings. The situation is complicated by the fact that the available theory for the understanding of nonlinear behavior of numerics is not at a stage to fully analyze the nonlinear Euler and Navier-Stokes equations. The discussion is based on the knowledge gained for nonlinear model problems with known analytical solutions to identify and explain the possible sources and remedies of numerical uncertainties in practical computations. Examples relevant to turbulent flow computations are included.

Overdetermined elliptic problems in topological disks

NASA Astrophysics Data System (ADS)

Mira, Pablo

2018-06-01

We introduce a method, based on the Poincaré-Hopf index theorem, to classify solutions to overdetermined problems for fully nonlinear elliptic equations in domains diffeomorphic to a closed disk. Applications to some well-known nonlinear elliptic PDEs are provided. Our result can be seen as the analogue of Hopf's uniqueness theorem for constant mean curvature spheres, but for the general analytic context of overdetermined elliptic problems.

Fully nonlinear heavy ion-acoustic solitary waves in astrophysical degenerate relativistic quantum plasmas

NASA Astrophysics Data System (ADS)

Sultana, S.; Schlickeiser, R.

2018-05-01

Fully nonlinear features of heavy ion-acoustic solitary waves (HIASWs) have been investigated in an astrophysical degenerate relativistic quantum plasma (ADRQP) containing relativistically degenerate electrons and non-relativistically degenerate light ion species, and non-degenerate heavy ion species. The pseudo-energy balance equation is derived from the fluid dynamical equations by adopting the well-known Sagdeev-potential approach, and the properties of arbitrary amplitude HIASWs are examined. The small amplitude limit for the propagation of HIASWs is also recovered. The basic features (width, amplitude, polarity, critical Mach number, speed, etc.) of HIASWs are found to be significantly modified by the relativistic effect of the electron species, and also by the variation of the number density of electron, light ion, and heavy ion species. The basic properties of HIASWs, that may propagated in some realistic astrophysical plasma systems (e.g., in white dwarfs), are briefly discussed.

Fully electromagnetic nonlinear gyrokinetic equations for tokamak edge turbulence

NASA Astrophysics Data System (ADS)

Hahm, T. S.; Wang, Lu; Madsen, J.

2009-02-01

An energy conserving set of the fully electromagnetic nonlinear gyrokinetic Vlasov equation and Maxwell's equations, which is applicable to both L-mode turbulence with large amplitude and H-mode turbulence in the presence of high E ×B shear has been derived. The phase-space action variational Lie perturbation method ensures the preservation of the conservation laws of the underlying Vlasov-Maxwell system. Generalized ordering takes ρi≪ρθi˜LE˜Lp≪R [here ρi is the thermal ion Larmor radius and ρθi=B /(Bθρi)], as typically observed in the tokamak H-mode edge, with LE and Lp being the radial electric field and pressure gradient lengths. k⊥ρi˜1 is assumed for generality, and the relative fluctuation amplitudes eδϕ /Ti˜δB/B are kept up to the second order. Extending the electrostatic theory in the presence of high E ×B shear [Hahm, Phys. Plasmas 3, 4658 (1996)], contributions of electromagnetic fluctuations to the particle charge density and current are explicitly evaluated via pullback transformation from the gyrocenter distribution function in the gyrokinetic Maxwell's equation.

Performance of a parallel algebraic multilevel preconditioner for stabilized finite element semiconductor device modeling

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

Lin, Paul T.; Shadid, John N.; Sala, Marzio

In this study results are presented for the large-scale parallel performance of an algebraic multilevel preconditioner for solution of the drift-diffusion model for semiconductor devices. The preconditioner is the key numerical procedure determining the robustness, efficiency and scalability of the fully-coupled Newton-Krylov based, nonlinear solution method that is employed for this system of equations. The coupled system is comprised of a source term dominated Poisson equation for the electric potential, and two convection-diffusion-reaction type equations for the electron and hole concentration. The governing PDEs are discretized in space by a stabilized finite element method. Solution of the discrete system ismore » obtained through a fully-implicit time integrator, a fully-coupled Newton-based nonlinear solver, and a restarted GMRES Krylov linear system solver. The algebraic multilevel preconditioner is based on an aggressive coarsening graph partitioning of the nonzero block structure of the Jacobian matrix. Representative performance results are presented for various choices of multigrid V-cycles and W-cycles and parameter variations for smoothers based on incomplete factorizations. Parallel scalability results are presented for solution of up to 10{sup 8} unknowns on 4096 processors of a Cray XT3/4 and an IBM POWER eServer system.« less

Turbulent Fluid Motion 6: Turbulence, Nonlinear Dynamics, and Deterministic Chaos

NASA Technical Reports Server (NTRS)

Deissler, Robert G.

1996-01-01

Several turbulent and nonturbulent solutions of the Navier-Stokes equations are obtained. The unaveraged equations are used numerically in conjunction with tools and concepts from nonlinear dynamics, including time series, phase portraits, Poincare sections, Liapunov exponents, power spectra, and strange attractors. Initially neighboring solutions for a low-Reynolds-number fully developed turbulence are compared. The turbulence is sustained by a nonrandom time-independent external force. The solutions, on the average, separate exponentially with time, having a positive Liapunov exponent. Thus, the turbulence is characterized as chaotic. In a search for solutions which contrast with the turbulent ones, the Reynolds number (or strength of the forcing) is reduced. Several qualitatively different flows are noted. These are, respectively, fully chaotic, complex periodic, weakly chaotic, simple periodic, and fixed-point. Of these, we classify only the fully chaotic flows as turbulent. Those flows have both a positive Liapunov exponent and Poincare sections without pattern. By contrast, the weakly chaotic flows, although having positive Liapunov exponents, have some pattern in their Poincare sections. The fixed-point and periodic flows are nonturbulent, since turbulence, as generally understood, is both time-dependent and aperiodic.

A harmonic polynomial cell (HPC) method for 3D Laplace equation with application in marine hydrodynamics

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

Shao, Yan-Lin, E-mail: yanlin.shao@dnvgl.com; Faltinsen, Odd M.

2014-10-01

We propose a new efficient and accurate numerical method based on harmonic polynomials to solve boundary value problems governed by 3D Laplace equation. The computational domain is discretized by overlapping cells. Within each cell, the velocity potential is represented by the linear superposition of a complete set of harmonic polynomials, which are the elementary solutions of Laplace equation. By its definition, the method is named as Harmonic Polynomial Cell (HPC) method. The characteristics of the accuracy and efficiency of the HPC method are demonstrated by studying analytical cases. Comparisons will be made with some other existing boundary element based methods,more » e.g. Quadratic Boundary Element Method (QBEM) and the Fast Multipole Accelerated QBEM (FMA-QBEM) and a fourth order Finite Difference Method (FDM). To demonstrate the applications of the method, it is applied to some studies relevant for marine hydrodynamics. Sloshing in 3D rectangular tanks, a fully-nonlinear numerical wave tank, fully-nonlinear wave focusing on a semi-circular shoal, and the nonlinear wave diffraction of a bottom-mounted cylinder in regular waves are studied. The comparisons with the experimental results and other numerical results are all in satisfactory agreement, indicating that the present HPC method is a promising method in solving potential-flow problems. The underlying procedure of the HPC method could also be useful in other fields than marine hydrodynamics involved with solving Laplace equation.« less

Solving Nonlinear Euler Equations with Arbitrary Accuracy

NASA Technical Reports Server (NTRS)

Dyson, Rodger W.

2005-01-01

A computer program that efficiently solves the time-dependent, nonlinear Euler equations in two dimensions to an arbitrarily high order of accuracy has been developed. The program implements a modified form of a prior arbitrary- accuracy simulation algorithm that is a member of the class of algorithms known in the art as modified expansion solution approximation (MESA) schemes. Whereas millions of lines of code were needed to implement the prior MESA algorithm, it is possible to implement the present MESA algorithm by use of one or a few pages of Fortran code, the exact amount depending on the specific application. The ability to solve the Euler equations to arbitrarily high accuracy is especially beneficial in simulations of aeroacoustic effects in settings in which fully nonlinear behavior is expected - for example, at stagnation points of fan blades, where linearizing assumptions break down. At these locations, it is necessary to solve the full nonlinear Euler equations, and inasmuch as the acoustical energy is of the order of 4 to 5 orders of magnitude below that of the mean flow, it is necessary to achieve an overall fractional error of less than 10-6 in order to faithfully simulate entropy, vortical, and acoustical waves.

Trajectory following and stabilization control of fully actuated AUV using inverse kinematics and self-tuning fuzzy PID.

PubMed

Hammad, Mohanad M; Elshenawy, Ahmed K; El Singaby, M I

2017-01-01

In this work a design for self-tuning non-linear Fuzzy Proportional Integral Derivative (FPID) controller is presented to control position and speed of Multiple Input Multiple Output (MIMO) fully-actuated Autonomous Underwater Vehicles (AUV) to follow desired trajectories. Non-linearity that results from the hydrodynamics and the coupled AUV dynamics makes the design of a stable controller a very difficult task. In this study, the control scheme in a simulation environment is validated using dynamic and kinematic equations for the AUV model and hydrodynamic damping equations. An AUV configuration with eight thrusters and an inverse kinematic model from a previous work is utilized in the simulation. In the proposed controller, Mamdani fuzzy rules are used to tune the parameters of the PID. Nonlinear fuzzy Gaussian membership functions are selected to give better performance and response in the non-linear system. A control architecture with two feedback loops is designed such that the inner loop is for velocity control and outer loop is for position control. Several test scenarios are executed to validate the controller performance including different complex trajectories with and without injection of ocean current disturbances. A comparison between the proposed FPID controller and the conventional PID controller is studied and shows that the FPID controller has a faster response to the reference signal and more stable behavior in a disturbed non-linear environment.

Trajectory following and stabilization control of fully actuated AUV using inverse kinematics and self-tuning fuzzy PID

PubMed Central

Elshenawy, Ahmed K.; El Singaby, M.I.

2017-01-01

In this work a design for self-tuning non-linear Fuzzy Proportional Integral Derivative (FPID) controller is presented to control position and speed of Multiple Input Multiple Output (MIMO) fully-actuated Autonomous Underwater Vehicles (AUV) to follow desired trajectories. Non-linearity that results from the hydrodynamics and the coupled AUV dynamics makes the design of a stable controller a very difficult task. In this study, the control scheme in a simulation environment is validated using dynamic and kinematic equations for the AUV model and hydrodynamic damping equations. An AUV configuration with eight thrusters and an inverse kinematic model from a previous work is utilized in the simulation. In the proposed controller, Mamdani fuzzy rules are used to tune the parameters of the PID. Nonlinear fuzzy Gaussian membership functions are selected to give better performance and response in the non-linear system. A control architecture with two feedback loops is designed such that the inner loop is for velocity control and outer loop is for position control. Several test scenarios are executed to validate the controller performance including different complex trajectories with and without injection of ocean current disturbances. A comparison between the proposed FPID controller and the conventional PID controller is studied and shows that the FPID controller has a faster response to the reference signal and more stable behavior in a disturbed non-linear environment. PMID:28683071

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

NASA Technical Reports Server (NTRS)

Marchesin, D.

1984-01-01

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

Exact Solutions for Stokes' Flow of a Non-Newtonian Nanofluid Model: A Lie Similarity Approach

NASA Astrophysics Data System (ADS)

Aziz, Taha; Aziz, A.; Khalique, C. M.

2016-07-01

The fully developed time-dependent flow of an incompressible, thermodynamically compatible non-Newtonian third-grade nanofluid is investigated. The classical Stokes model is considered in which the flow is generated due to the motion of the plate in its own plane with an impulsive velocity. The Lie symmetry approach is utilised to convert the governing nonlinear partial differential equation into different linear and nonlinear ordinary differential equations. The reduced ordinary differential equations are then solved by using the compatibility and generalised group method. Exact solutions for the model equation are deduced in the form of closed-form exponential functions which are not available in the literature before. In addition, we also derived the conservation laws associated with the governing model. Finally, the physical features of the pertinent parameters are discussed in detail through several graphs.

Tripartite correlations over two octaves from cascaded harmonic generation

NASA Astrophysics Data System (ADS)

Olsen, M. K.

2018-03-01

We analyse the output quantum tripartite correlations from an intracavity nonlinear optical system which uses cascaded nonlinearities to produce both second and fourth harmonic outputs from an input field at the fundamental frequency. Using fully quantum equations of motion, we investigate two parameter regimes and show that the system produces tripartite inseparability, entanglement and EPR steering, with the detection of these depending on the correlations being considered.

Bifurcation Analysis of an Electrostatically Actuated Nano-Beam Based on Modified Couple Stress Theory

NASA Astrophysics Data System (ADS)

Rezaei Kivi, Araz; Azizi, Saber; Norouzi, Peyman

2017-12-01

In this paper, the nonlinear size-dependent static and dynamic behavior of an electrostatically actuated nano-beam is investigated. A fully clamped nano-beam is considered for the modeling of the deformable electrode of the NEMS. The governing differential equation of the motion is derived using Hamiltonian principle based on couple stress theory; a non-classical theory for considering length scale effects. The nonlinear partial differential equation of the motion is discretized to a nonlinear Duffing type ODE's using Galerkin method. Static and dynamic pull-in instabilities obtained by both classical theory and MCST are compared. At the second stage of analysis, shooting technique is utilized to obtain the frequency response curve, and to capture the periodic solutions of the motion; the stability of the periodic solutions are gained by Floquet theory. The nonlinear dynamic behavior of the deformable electrode due to the AC harmonic accompanied with size dependency is investigated.

Development of an Integrated Modeling Framework for Simulations of Coastal Processes in Deltaic Environments Using High-Performance Computing

DTIC Science & Technology

2008-01-01

exceeds the local water depth. The approximation eliminates the vertical dimension of the elliptic equation that is normally required for the fully non...used for vertical resolution. The shallow water equations (SWE) are a set of non-linear hyperbolic equations. As the equations are derived under...linear standing wave with a wavelength of 10 m in a square 10 m by 10 m basin. The still water depth is 0.5 m. In order to compare with the analytical

Nonlinear travelling waves in rotating Hagen–Poiseuille flow

NASA Astrophysics Data System (ADS)

Pier, Benoît; Govindarajan, Rama

2018-03-01

The dynamics of viscous flow through a rotating pipe is considered. Small-amplitude stability characteristics are obtained by linearizing the Navier–Stokes equations around the base flow and solving the resulting eigenvalue problems. For linearly unstable configurations, the dynamics leads to fully developed finite-amplitude perturbations that are computed by direct numerical simulations of the complete Navier–Stokes equations. By systematically investigating all linearly unstable combinations of streamwise wave number k and azimuthal mode number m, for streamwise Reynolds numbers {{Re}}z ≤slant 500 and rotational Reynolds numbers {{Re}}{{Ω }} ≤slant 500, the complete range of nonlinear travelling waves is obtained and the associated flow fields are characterized.

Effect of Second-Order and Fully Nonlinear Wave Kinematics on a Tension-Leg-Platform Wind Turbine in Extreme Wave Conditions: Preprint

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

Robertson, Amy N; Jonkman, Jason; Pegalajar-Jurado, Antonio

In this study, we assess the impact of different wave kinematics models on the dynamic response of a tension-leg-platform wind turbine. Aero-hydro-elastic simulations of the floating wind turbine are carried out employing linear, second-order, and fully nonlinear kinematics using the Morison equation for the hydrodynamic forcing. The wave kinematics are computed from either theoretical or measured signals of free-surface elevation. The numerical results from each model are compared to results from wave basin tests on a scaled prototype. The comparison shows that sub and superharmonic responses can be introduced by second-order and fully nonlinear wave kinematics. The response at themore » wave frequency range is better reproduced when kinematics are generated from the measured surface elevation. In the future, the numerical response may be further improved by replacing the global, constant damping coefficients in the model by a more detailed, customizable definition of the user-defined numerical damping.« less

Effect of Second-Order and Fully Nonlinear Wave Kinematics on a Tension-Leg-Platform Wind Turbine in Extreme Wave Conditions

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

Robertson, Amy N; Jonkman, Jason; Pegalajar-Jurado, Antonio

In this study, we assess the impact of different wave kinematics models on the dynamic response of a tension-leg-platform wind turbine. Aero-hydro-elastic simulations of the floating wind turbine are carried out employing linear, second-order, and fully nonlinear kinematics using the Morison equation for the hydrodynamic forcing. The wave kinematics are computed from either theoretical or measured signals of free-surface elevation. The numerical results from each model are compared to results from wave basin tests on a scaled prototype. The comparison shows that sub and superharmonic responses can be introduced by second-order and fully nonlinear wave kinematics. The response at themore » wave frequency range is better reproduced when kinematics are generated from the measured surface elevation. In the future, the numerical response may be further improved by replacing the global, constant damping coefficients in the model by a more detailed, customizable definition of the user-defined numerical damping.« 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.

A high-order relaxation method with projective integration for solving nonlinear systems of hyperbolic conservation laws

NASA Astrophysics Data System (ADS)

Lafitte, Pauline; Melis, Ward; Samaey, Giovanni

2017-07-01

We present a general, high-order, fully explicit relaxation scheme which can be applied to any system of nonlinear hyperbolic conservation laws in multiple dimensions. The scheme consists of two steps. In a first (relaxation) step, the nonlinear hyperbolic conservation law is approximated by a kinetic equation with stiff BGK source term. Then, this kinetic equation is integrated in time using a projective integration method. After taking a few small (inner) steps with a simple, explicit method (such as direct forward Euler) to damp out the stiff components of the solution, the time derivative is estimated and used in an (outer) Runge-Kutta method of arbitrary order. We show that, with an appropriate choice of inner step size, the time step restriction on the outer time step is similar to the CFL condition for the hyperbolic conservation law. Moreover, the number of inner time steps is also independent of the stiffness of the BGK source term. We discuss stability and consistency, and illustrate with numerical results (linear advection, Burgers' equation and the shallow water and Euler equations) in one and two spatial dimensions.

Effects of group velocity and multiplasmon resonances on the modulation of Langmuir waves in a degenerate plasma

NASA Astrophysics Data System (ADS)

Misra, Amar P.; Chatterjee, Debjani; Brodin, Gert

2017-11-01

We study the nonlinear wave modulation of Langmuir waves (LWs) in a fully degenerate plasma. Using the Wigner-Moyal equation coupled to the Poisson equation and the multiple scale expansion technique, a modified nonlocal nonlinear Schrödinger (NLS) equation is derived which governs the evolution of LW envelopes in degenerate plasmas. The nonlocal nonlinearity in the NLS equation appears due to the group velocity and multiplasmon resonances, i.e., resonances induced by the simultaneous particle absorption of multiple wave quanta. We focus on the regime where the resonant velocity of electrons is larger than the Fermi velocity and thereby the linear Landau damping is forbidden. As a result, the nonlinear wave-particle resonances due to the group velocity and multiplasmon processes are the dominant mechanisms for wave-particle interaction. It is found that in contrast to classical or semiclassical plasmas, the group velocity resonance does not necessarily give rise the wave damping in the strong quantum regime where ℏ k ˜m vF with ℏ denoting the reduced Planck's constant, m the electron mass, and vF the Fermi velocity; however, the three-plasmon process plays a dominant role in the nonlinear Landau damping of wave envelopes. In this regime, the decay rate of the wave amplitude is also found to be higher compared to that in the modest quantum regime where the multiplasmon effects are forbidden.

The Evolution of Finite Amplitude Wavetrains in Plane Channel Flow

NASA Technical Reports Server (NTRS)

Hewitt, R. E.; Hall, P.

1996-01-01

We consider a viscous incompressible fluid flow driven between two parallel plates by a constant pressure gradient. The flow is at a finite Reynolds number, with an 0(l) disturbance in the form of a traveling wave. A phase equation approach is used to discuss the evolution of slowly varying fully nonlinear two dimensional wavetrains. We consider uniform wavetrains in detail, showing that the development of a wavenumber perturbation is governed by Burgers equation in most cases. The wavenumber perturbation theory, constructed using the phase equation approach for a uniform wavetrain, is shown to be distinct from an amplitude perturbation expansion about the periodic flow. In fact we show that the amplitude equation contains only linear terms and is simply the heat equation. We review, briefly, the well known dynamics of Burgers equation, which imply that both shock structures and finite time singularities of the wavenumber perturbation can occur with respect to the slow scales. Numerical computations have been performed to identify areas of the (wavenumber, Reynolds number, energy) neutral surface for which each of these possibilities can occur. We note that the evolution equations will breakdown under certain circumstances, in particular for a weakly nonlinear secondary flow. Finally we extend the theory to three dimensions and discuss the limit of a weak spanwise dependence for uniform wavetrains, showing that two functions are required to describe the evolution. These unknowns are a phase and a pressure function which satisfy a pair of linearly coupled partial differential equations. The results obtained from applying the same analysis to the fully three dimensional problem are included as an appendix.

A coupled electro-thermal Discontinuous Galerkin method

NASA Astrophysics Data System (ADS)

Homsi, L.; Geuzaine, C.; Noels, L.

2017-11-01

This paper presents a Discontinuous Galerkin scheme in order to solve the nonlinear elliptic partial differential equations of coupled electro-thermal problems. In this paper we discuss the fundamental equations for the transport of electricity and heat, in terms of macroscopic variables such as temperature and electric potential. A fully coupled nonlinear weak formulation for electro-thermal problems is developed based on continuum mechanics equations expressed in terms of energetically conjugated pair of fluxes and fields gradients. The weak form can thus be formulated as a Discontinuous Galerkin method. The existence and uniqueness of the weak form solution are proved. The numerical properties of the nonlinear elliptic problems i.e., consistency and stability, are demonstrated under specific conditions, i.e. use of high enough stabilization parameter and at least quadratic polynomial approximations. Moreover the prior error estimates in the H1-norm and in the L2-norm are shown to be optimal in the mesh size with the polynomial approximation degree.

SIERRA Multimechanics Module: Aria User Manual Version 4.44

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

Sierra Thermal /Fluid Team

2017-04-01

Aria is a Galerkin fnite element based program for solving coupled-physics problems described by systems of PDEs and is capable of solving nonlinear, implicit, transient and direct-to-steady state problems in two and three dimensions on parallel architectures. The suite of physics currently supported by Aria includes thermal energy transport, species transport, and electrostatics as well as generalized scalar, vector and tensor transport equations. Additionally, Aria includes support for manufacturing process fows via the incompressible Navier-Stokes equations specialized to a low Reynolds number ( %3C 1 ) regime. Enhanced modeling support of manufacturing processing is made possible through use of eithermore » arbitrary Lagrangian- Eulerian (ALE) and level set based free and moving boundary tracking in conjunction with quasi-static nonlinear elastic solid mechanics for mesh control. Coupled physics problems are solved in several ways including fully-coupled Newton's method with analytic or numerical sensitivities, fully-coupled Newton- Krylov methods and a loosely-coupled nonlinear iteration about subsets of the system that are solved using combinations of the aforementioned methods. Error estimation, uniform and dynamic h -adaptivity and dynamic load balancing are some of Aria's more advanced capabilities. Aria is based upon the Sierra Framework.« less

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

Sierra Thermal/Fluid Team

Aria is a Galerkin fnite element based program for solving coupled-physics problems described by systems of PDEs and is capable of solving nonlinear, implicit, transient and direct-to-steady state problems in two and three dimensions on parallel architectures. The suite of physics currently supported by Aria includes thermal energy transport, species transport, and electrostatics as well as generalized scalar, vector and tensor transport equations. Additionally, Aria includes support for manufacturing process fows via the incompressible Navier-Stokes equations specialized to a low Reynolds number ( %3C 1 ) regime. Enhanced modeling support of manufacturing processing is made possible through use of eithermore » arbitrary Lagrangian- Eulerian (ALE) and level set based free and moving boundary tracking in conjunction with quasi-static nonlinear elastic solid mechanics for mesh control. Coupled physics problems are solved in several ways including fully-coupled Newton's method with analytic or numerical sensitivities, fully-coupled Newton- Krylov methods and a loosely-coupled nonlinear iteration about subsets of the system that are solved using combinations of the aforementioned methods. Error estimation, uniform and dynamic h -adaptivity and dynamic load balancing are some of Aria's more advanced capabilities. Aria is based upon the Sierra Framework.« less

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

Sierra Thermal /Fluid Team

Aria is a Galerkin finite element based program for solving coupled-physics problems described by systems of PDEs and is capable of solving nonlinear, implicit, transient and direct-to-steady state problems in two and three dimensions on parallel architectures. The suite of physics currently supported by Aria includes thermal energy transport, species transport, and electrostatics as well as generalized scalar, vector and tensor transport equations. Additionally, Aria includes support for manufacturing process flows via the incompressible Navier-Stokes equations specialized to a low Reynolds number (Re %3C 1) regime. Enhanced modeling support of manufacturing processing is made possible through use of either arbitrarymore » Lagrangian- Eulerian (ALE) and level set based free and moving boundary tracking in conjunction with quasi-static nonlinear elastic solid mechanics for mesh control. Coupled physics problems are solved in several ways including fully-coupled Newton's method with analytic or numerical sensitivities, fully-coupled Newton- Krylov methods and a loosely-coupled nonlinear iteration about subsets of the system that are solved using combinations of the aforementioned methods. Error estimation, uniform and dynamic h-adaptivity and dynamic load balancing are some of Aria's more advanced capabilities. Aria is based upon the Sierra Framework.« less

Exact closed-form solutions of a fully nonlinear asymptotic two-fluid model

NASA Astrophysics Data System (ADS)

Cheviakov, Alexei F.

2018-05-01

A fully nonlinear model of Choi and Camassa (1999) describing one-dimensional incompressible dynamics of two non-mixing fluids in a horizontal channel, under a shallow water approximation, is considered. An equivalence transformation is presented, leading to a special dimensionless form of the system, involving a single dimensionless constant physical parameter, as opposed to five parameters present in the original model. A first-order dimensionless ordinary differential equation describing traveling wave solutions is analyzed. Several multi-parameter families of physically meaningful exact closed-form solutions of the two-fluid model are derived, corresponding to periodic, solitary, and kink-type bidirectional traveling waves; specific examples are given, and properties of the exact solutions are analyzed.

Neoclassical simulation of tokamak plasmas using the continuum gyrokinetic code TEMPEST.

PubMed

Xu, X Q

2008-07-01

We present gyrokinetic neoclassical simulations of tokamak plasmas with a self-consistent electric field using a fully nonlinear (full- f ) continuum code TEMPEST in a circular geometry. A set of gyrokinetic equations are discretized on a five-dimensional computational grid in phase space. The present implementation is a method of lines approach where the phase-space derivatives are discretized with finite differences, and implicit backward differencing formulas are used to advance the system in time. The fully nonlinear Boltzmann model is used for electrons. The neoclassical electric field is obtained by solving the gyrokinetic Poisson equation with self-consistent poloidal variation. With a four-dimensional (psi,theta,micro) version of the TEMPEST code, we compute the radial particle and heat fluxes, the geodesic-acoustic mode, and the development of the neoclassical electric field, which we compare with neoclassical theory using a Lorentz collision model. The present work provides a numerical scheme for self-consistently studying important dynamical aspects of neoclassical transport and electric field in toroidal magnetic fusion devices.

Neoclassical simulation of tokamak plasmas using the continuum gyrokinetic code TEMPEST

NASA Astrophysics Data System (ADS)

Xu, X. Q.

2008-07-01

We present gyrokinetic neoclassical simulations of tokamak plasmas with a self-consistent electric field using a fully nonlinear (full- f ) continuum code TEMPEST in a circular geometry. A set of gyrokinetic equations are discretized on a five-dimensional computational grid in phase space. The present implementation is a method of lines approach where the phase-space derivatives are discretized with finite differences, and implicit backward differencing formulas are used to advance the system in time. The fully nonlinear Boltzmann model is used for electrons. The neoclassical electric field is obtained by solving the gyrokinetic Poisson equation with self-consistent poloidal variation. With a four-dimensional (ψ,θ,γ,μ) version of the TEMPEST code, we compute the radial particle and heat fluxes, the geodesic-acoustic mode, and the development of the neoclassical electric field, which we compare with neoclassical theory using a Lorentz collision model. The present work provides a numerical scheme for self-consistently studying important dynamical aspects of neoclassical transport and electric field in toroidal magnetic fusion devices.

Discrete approach to stochastic parametrization and dimension reduction in nonlinear dynamics.

PubMed

Chorin, Alexandre J; Lu, Fei

2015-08-11

Many physical systems are described by nonlinear differential equations that are too complicated to solve in full. A natural way to proceed is to divide the variables into those that are of direct interest and those that are not, formulate solvable approximate equations for the variables of greater interest, and use data and statistical methods to account for the impact of the other variables. In the present paper we consider time-dependent problems and introduce a fully discrete solution method, which simplifies both the analysis of the data and the numerical algorithms. The resulting time series are identified by a NARMAX (nonlinear autoregression moving average with exogenous input) representation familiar from engineering practice. The connections with the Mori-Zwanzig formalism of statistical physics are discussed, as well as an application to the Lorenz 96 system.

Fully- and weakly-nonlinear biperiodic traveling waves in shallow water

NASA Astrophysics Data System (ADS)

Hirakawa, Tomoaki; Okamura, Makoto

2018-04-01

We directly calculate fully nonlinear traveling waves that are periodic in two independent horizontal directions (biperiodic) in shallow water. Based on the Riemann theta function, we also calculate exact periodic solutions to the Kadomtsev-Petviashvili (KP) equation, which can be obtained by assuming weakly-nonlinear, weakly-dispersive, weakly-two-dimensional waves. To clarify how the accuracy of the biperiodic KP solution is affected when some of the KP approximations are not satisfied, we compare the fully- and weakly-nonlinear periodic traveling waves of various wave amplitudes, wave depths, and interaction angles. As the interaction angle θ decreases, the wave frequency and the maximum wave height of the biperiodic KP solution both increase, and the central peak sharpens and grows beyond the height of the corresponding direct numerical solutions, indicating that the biperiodic KP solution cannot qualitatively model direct numerical solutions for θ ≲ 45^\\circ . To remedy the weak two-dimensionality approximation, we apply the correction of Yeh et al (2010 Eur. Phys. J. Spec. Top. 185 97-111) to the biperiodic KP solution, which substantially improves the solution accuracy and results in wave profiles that are indistinguishable from most other cases.

Development and Verification of the Charring Ablating Thermal Protection Implicit System Solver

NASA Technical Reports Server (NTRS)

Amar, Adam J.; Calvert, Nathan D.; Kirk, Benjamin S.

2010-01-01

The development and verification of the Charring Ablating Thermal Protection Implicit System Solver is presented. This work concentrates on the derivation and verification of the stationary grid terms in the equations that govern three-dimensional heat and mass transfer for charring thermal protection systems including pyrolysis gas flow through the porous char layer. The governing equations are discretized according to the Galerkin finite element method with first and second order implicit time integrators. The governing equations are fully coupled and are solved in parallel via Newton's method, while the fully implicit linear system is solved with the Generalized Minimal Residual method. Verification results from exact solutions and the Method of Manufactured Solutions are presented to show spatial and temporal orders of accuracy as well as nonlinear convergence rates.

Development and Verification of the Charring, Ablating Thermal Protection Implicit System Simulator

NASA Technical Reports Server (NTRS)

Amar, Adam J.; Calvert, Nathan; Kirk, Benjamin S.

2011-01-01

The development and verification of the Charring Ablating Thermal Protection Implicit System Solver (CATPISS) is presented. This work concentrates on the derivation and verification of the stationary grid terms in the equations that govern three-dimensional heat and mass transfer for charring thermal protection systems including pyrolysis gas flow through the porous char layer. The governing equations are discretized according to the Galerkin finite element method (FEM) with first and second order fully implicit time integrators. The governing equations are fully coupled and are solved in parallel via Newton s method, while the linear system is solved via the Generalized Minimum Residual method (GMRES). Verification results from exact solutions and Method of Manufactured Solutions (MMS) are presented to show spatial and temporal orders of accuracy as well as nonlinear convergence rates.

A Fast Solver for Implicit Integration of the Vlasov--Poisson System in the Eulerian Framework

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

Garrett, C. Kristopher; Hauck, Cory D.

In this paper, we present a domain decomposition algorithm to accelerate the solution of Eulerian-type discretizations of the linear, steady-state Vlasov equation. The steady-state solver then forms a key component in the implementation of fully implicit or nearly fully implicit temporal integrators for the nonlinear Vlasov--Poisson system. The solver relies on a particular decomposition of phase space that enables the use of sweeping techniques commonly used in radiation transport applications. The original linear system for the phase space unknowns is then replaced by a smaller linear system involving only unknowns on the boundary between subdomains, which can then be solvedmore » efficiently with Krylov methods such as GMRES. Steady-state solves are combined to form an implicit Runge--Kutta time integrator, and the Vlasov equation is coupled self-consistently to the Poisson equation via a linearized procedure or a nonlinear fixed-point method for the electric field. Finally, numerical results for standard test problems demonstrate the efficiency of the domain decomposition approach when compared to the direct application of an iterative solver to the original linear system.« less

A Fast Solver for Implicit Integration of the Vlasov--Poisson System in the Eulerian Framework

DOE PAGES

Garrett, C. Kristopher; Hauck, Cory D.

2018-04-05

In this paper, we present a domain decomposition algorithm to accelerate the solution of Eulerian-type discretizations of the linear, steady-state Vlasov equation. The steady-state solver then forms a key component in the implementation of fully implicit or nearly fully implicit temporal integrators for the nonlinear Vlasov--Poisson system. The solver relies on a particular decomposition of phase space that enables the use of sweeping techniques commonly used in radiation transport applications. The original linear system for the phase space unknowns is then replaced by a smaller linear system involving only unknowns on the boundary between subdomains, which can then be solvedmore » efficiently with Krylov methods such as GMRES. Steady-state solves are combined to form an implicit Runge--Kutta time integrator, and the Vlasov equation is coupled self-consistently to the Poisson equation via a linearized procedure or a nonlinear fixed-point method for the electric field. Finally, numerical results for standard test problems demonstrate the efficiency of the domain decomposition approach when compared to the direct application of an iterative solver to the original linear system.« less

Coarse-grained description of cosmic structure from Szekeres models

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

Sussman, Roberto A.; Gaspar, I. Delgado; Hidalgo, Juan Carlos, E-mail: sussman@nucleares.unam.mx, E-mail: ismael.delgadog@uaem.edu.mx, E-mail: hidalgo@fis.unam.mx

2016-03-01

We show that the full dynamical freedom of the well known Szekeres models allows for the description of elaborated 3-dimensional networks of cold dark matter structures (over-densities and/or density voids) undergoing ''pancake'' collapse. By reducing Einstein's field equations to a set of evolution equations, which themselves reduce in the linear limit to evolution equations for linear perturbations, we determine the dynamics of such structures, with the spatial comoving location of each structure uniquely specified by standard early Universe initial conditions. By means of a representative example we examine in detail the density contrast, the Hubble flow and peculiar velocities ofmore » structures that evolved, from linear initial data at the last scattering surface, to fully non-linear 10–20 Mpc scale configurations today. To motivate further research, we provide a qualitative discussion on the connection of Szekeres models with linear perturbations and the pancake collapse of the Zeldovich approximation. This type of structure modelling provides a coarse grained—but fully relativistic non-linear and non-perturbative —description of evolving large scale cosmic structures before their virialisation, and as such it has an enormous potential for applications in cosmological research.« less

The method of projected characteristics for the evolution of magnetic arches

NASA Technical Reports Server (NTRS)

Nakagawa, Y.; Hu, Y. Q.; Wu, S. T.

1987-01-01

A numerical method of solving fully nonlinear MHD equation is described. In particular, the formulation based on the newly developed method of projected characteristics (Nakagawa, 1981) suitable to study the evolution of magnetic arches due to motions of their foot-points is presented. The final formulation is given in the form of difference equations; therefore, the analysis of numerical stability is also presented. Further, the most important derivation of physically self-consistent, time-dependent boundary conditions (i.e. the evolving boundary equations) is given in detail, and some results obtained with such boundary equations are reported.

Nonlinear excitations for the positron acoustic shock waves in dissipative nonextensive electron-positron-ion plasmas

NASA Astrophysics Data System (ADS)

Saha, Asit

2017-03-01

Positron acoustic shock waves (PASHWs) in unmagnetized electron-positron-ion (e-p-i) plasmas consisting of mobile cold positrons, immobile positive ions, q-nonextensive distributed electrons, and hot positrons are studied. The cold positron kinematic viscosity is considered and the reductive perturbation technique is used to derive the Burgers equation. Applying traveling wave transformation, the Burgers equation is transformed to a one dimensional dynamical system. All possible vector fields corresponding to the dynamical system are presented. We have analyzed the dynamical system with the help of potential energy, which helps to identify the stability and instability of the equilibrium points. It is found that the viscous force acting on cold mobile positron fluid is a source of dissipation and is responsible for the formation of the PASHWs. Furthermore, fully nonlinear arbitrary amplitude positron acoustic waves are also studied applying the theory of planar dynamical systems. It is also observed that the fundamental features of the small amplitude and arbitrary amplitude PASHWs are significantly affected by the effect of the physical parameters q e , q h , μ e , μ h , σ , η , and U. This work can be useful to understand the qualitative changes in the dynamics of nonlinear small amplitude and fully nonlinear arbitrary amplitude PASHWs in solar wind, ionosphere, lower part of magnetosphere, and auroral acceleration regions.

Documentation of computer program VS2D to solve the equations of fluid flow in variably saturated porous media

USGS Publications Warehouse

Lappala, E.G.; Healy, R.W.; Weeks, E.P.

1987-01-01

This report documents FORTRAN computer code for solving problems involving variably saturated single-phase flow in porous media. The flow equation is written with total hydraulic potential as the dependent variable, which allows straightforward treatment of both saturated and unsaturated conditions. The spatial derivatives in the flow equation are approximated by central differences, and time derivatives are approximated either by a fully implicit backward or by a centered-difference scheme. Nonlinear conductance and storage terms may be linearized using either an explicit method or an implicit Newton-Raphson method. Relative hydraulic conductivity is evaluated at cell boundaries by using either full upstream weighting, the arithmetic mean, or the geometric mean of values from adjacent cells. Nonlinear boundary conditions treated by the code include infiltration, evaporation, and seepage faces. Extraction by plant roots that is caused by atmospheric demand is included as a nonlinear sink term. These nonlinear boundary and sink terms are linearized implicitly. The code has been verified for several one-dimensional linear problems for which analytical solutions exist and against two nonlinear problems that have been simulated with other numerical models. A complete listing of data-entry requirements and data entry and results for three example problems are provided. (USGS)

Numerical simulation of solitary waves on deep water with constant vorticity

NASA Astrophysics Data System (ADS)

Dosaev, A. S.; Shishina, M. I.; Troitskaya, Yu I.

2018-01-01

Characteristics of solitary deep water waves on a flow with constant vorticity are investigated by numerical simulation within the framework of fully nonlinear equations of motion (Euler equations) using the method of surface-tracking conformal coordinates. To ensure that solutions observed are stable, soliton formation as a result of disintegration of an initial pulse-like disturbance is modeled. Evidence is obtained that solitary waves with height above a certain threshold are unstable.

Nearly spherical constant power detonation waves driven by focused radiation

NASA Technical Reports Server (NTRS)

George, Y. H.

1973-01-01

Shape and inner flow of a tridimensional spark are studied. The spark is created by focusing a laser beam in a gas. A second order fully non-linear equation is derived for the radial velocity on the axis of symmetry in the neighborhood of the origin. Solutions to that equation display the existence of a forbidden region near the focus, thus indicating the limits of applicability of a small perturbation solution.

Exploring New Physics Frontiers Through Numerical Relativity.

PubMed

Cardoso, Vitor; Gualtieri, Leonardo; Herdeiro, Carlos; Sperhake, Ulrich

2015-01-01

The demand to obtain answers to highly complex problems within strong-field gravity has been met with significant progress in the numerical solution of Einstein's equations - along with some spectacular results - in various setups. We review techniques for solving Einstein's equations in generic spacetimes, focusing on fully nonlinear evolutions but also on how to benchmark those results with perturbative approaches. The results address problems in high-energy physics, holography, mathematical physics, fundamental physics, astrophysics and cosmology.

SPHERES tethered formation flight testbed: application to NASA's SPECS mission

NASA Astrophysics Data System (ADS)

Chung, Soon-Jo; Kong, Edmund M.; Miller, David W.

2005-08-01

This paper elaborates on theory and experiment of the formation flight control for the future space-borne tethered interferometers. The nonlinear equations of multi-vehicle tethered spacecraft system are derived by Lagrange equations and decoupling method. The preliminary analysis predicts unstable dynamics depending on the direction of the tether motor. The controllability analysis indicates that both array resizing and spin-up are fully controllable only by the reaction wheels and the tether motor, thereby eliminating the need for thrusters. Linear and nonlinear decentralized control techniques have been implemented into the tethered SPHERES testbed, and tested at the NASA MSFC's flat floor facility using two and three SPHERES configurations. The nonlinear control using feedback linearization technique performed successfully in both two SPHERES in-line configuration and three triangular configuration while varying the tether length. The relative metrology system, using the ultra sound metrology system and the inertial sensors as well as the decentralized nonlinear estimator, is developed to provide necessary state information.

Coupled nonlinear aeroelasticity and flight dynamics of fully flexible aircraft

NASA Astrophysics Data System (ADS)

Su, Weihua

This dissertation introduces an approach to effectively model and analyze the coupled nonlinear aeroelasticity and flight dynamics of highly flexible aircraft. A reduced-order, nonlinear, strain-based finite element framework is used, which is capable of assessing the fundamental impact of structural nonlinear effects in preliminary vehicle design and control synthesis. The cross-sectional stiffness and inertia properties of the wings are calculated along the wing span, and then incorporated into the one-dimensional nonlinear beam formulation. Finite-state unsteady subsonic aerodynamics is used to compute airloads along lifting surfaces. Flight dynamic equations are then introduced to complete the aeroelastic/flight dynamic system equations of motion. Instead of merely considering the flexibility of the wings, the current work allows all members of the vehicle to be flexible. Due to their characteristics of being slender structures, the wings, tail, and fuselage of highly flexible aircraft can be modeled as beams undergoing three dimensional displacements and rotations. New kinematic relationships are developed to handle the split beam systems, such that fully flexible vehicles can be effectively modeled within the existing framework. Different aircraft configurations are modeled and studied, including Single-Wing, Joined-Wing, Blended-Wing-Body, and Flying-Wing configurations. The Lagrange Multiplier Method is applied to model the nodal displacement constraints at the joint locations. Based on the proposed models, roll response and stability studies are conducted on fully flexible and rigidized models. The impacts of the flexibility of different vehicle members on flutter with rigid body motion constraints, flutter in free flight condition, and roll maneuver performance are presented. Also, the static stability of the compressive member of the Joined-Wing configuration is studied. A spatially-distributed discrete gust model is incorporated into the time simulation of the framework. Gust responses of the Flying-Wing configuration subject to stall effects are investigated. A bilinear torsional stiffness model is introduced to study the skin wrinkling due to large bending curvature of the Flying-Wing. The numerical studies illustrate the improvements of the existing reduced-order formulation with new capabilities of both structural modeling and coupled aeroelastic and flight dynamic analysis of fully flexible aircraft.

Nonlinear interaction between underwater explosion bubble and structure based on fully coupled model

NASA Astrophysics Data System (ADS)

Zhang, A. M.; Wu, W. B.; Liu, Y. L.; Wang, Q. X.

2017-08-01

The interaction between an underwater explosion bubble and an elastic-plastic structure is a complex transient process, accompanying violent bubble collapsing, jet impact, penetration through the bubble, and large structural deformation. In the present study, the bubble dynamics are modeled using the boundary element method and the nonlinear transient structural response is modeled using the explicit finite element method. A new fully coupled 3D model is established through coupling the equations for the state variables of the fluid and structure and solving them as a set of coupled linear algebra equations. Based on the acceleration potential theory, the mutual dependence between the hydrodynamic load and the structural motion is decoupled. The pressure distribution in the flow field is calculated with the Bernoulli equation, where the partial derivative of the velocity potential in time is calculated using the boundary integral method to avoid numerical instabilities. To validate the present fully coupled model, the experiments of small-scale underwater explosion near a stiffened plate are carried out. High-speed imaging is used to capture the bubble behaviors and strain gauges are used to measure the strain response. The numerical results correspond well with the experimental data, in terms of bubble shapes and structural strain response. By both the loosely coupled model and the fully coupled model, the interaction between a bubble and a hollow spherical shell is studied. The bubble patterns vary with different parameters. When the fully coupled model and the loosely coupled model are advanced with the same time step, the error caused by the loosely coupled model becomes larger with the coupling effect becoming stronger. The fully coupled model is more stable than the loosely coupled model. Besides, the influences of the internal fluid on the dynamic response of the spherical shell are studied. At last, the case that the bubble interacts with an air-backed stiffened plate is simulated. The associated interesting physical phenomenon is obtained and expounded.

Accurate electrostatic and van der Waals pull-in prediction for fully clamped nano/micro-beams using linear universal graphs of pull-in instability

NASA Astrophysics Data System (ADS)

Tahani, Masoud; Askari, Amir R.

2014-09-01

In spite of the fact that pull-in instability of electrically actuated nano/micro-beams has been investigated by many researchers to date, no explicit formula has been presented yet which can predict pull-in voltage based on a geometrically non-linear and distributed parameter model. The objective of present paper is to introduce a simple and accurate formula to predict this value for a fully clamped electrostatically actuated nano/micro-beam. To this end, a non-linear Euler-Bernoulli beam model is employed, which accounts for the axial residual stress, geometric non-linearity of mid-plane stretching, distributed electrostatic force and the van der Waals (vdW) attraction. The non-linear boundary value governing equation of equilibrium is non-dimensionalized and solved iteratively through single-term Galerkin based reduced order model (ROM). The solutions are validated thorough direct comparison with experimental and other existing results reported in previous studies. Pull-in instability under electrical and vdW loads are also investigated using universal graphs. Based on the results of these graphs, non-dimensional pull-in and vdW parameters, which are defined in the text, vary linearly versus the other dimensionless parameters of the problem. Using this fact, some linear equations are presented to predict pull-in voltage, the maximum allowable length, the so-called detachment length, and the minimum allowable gap for a nano/micro-system. These linear equations are also reduced to a couple of universal pull-in formulas for systems with small initial gap. The accuracy of the universal pull-in formulas are also validated by comparing its results with available experimental and some previous geometric linear and closed-form findings published in the literature.

On Flexible Tubes Conveying Fluid: Geometric Nonlinear Theory, Stability and Dynamics

NASA Astrophysics Data System (ADS)

Gay-Balmaz, François; Putkaradze, Vakhtang

2015-08-01

We derive a fully three-dimensional, geometrically exact theory for flexible tubes conveying fluid. The theory also incorporates the change of the cross section available to the fluid motion during the dynamics. Our approach is based on the symmetry-reduced, exact geometric description for elastic rods, coupled with the fluid transport and subject to the volume conservation constraint for the fluid. We first derive the equations of motion directly, by using an Euler-Poincaré variational principle. We then justify this derivation with a more general theory elucidating the interesting mathematical concepts appearing in this problem, such as partial left (elastic) and right (fluid) invariance of the system, with the added holonomic constraint (volume). We analyze the fully nonlinear behavior of the model when the axis of the tube remains straight. We then proceed to the linear stability analysis and show that our theory introduces important corrections to previously derived results, both in the consistency at all wavelength and in the effects arising from the dynamical change of the cross section. Finally, we derive and analyze several analytical, fully nonlinear solutions of traveling wave type in two dimensions.

Laboratory tests of short intense envelope solitons

NASA Astrophysics Data System (ADS)

Slunyaev, A.; Clauss, G. F.; Klein, M.; Onorato, M.

2012-04-01

Stability of short intense nonlinear wave groups propagating over deep water is tested in laboratory runs which are performed in the facility of the Technical University of Berlin. The strongly nonlinear simulation of quasi-steady nonlinear wave groups within the framework of the Euler equations is used to generate the surface elevation time series at a border of the water tank. Besides, the exact analytic solution of the nonlinear Schrodinger equation is used for this purpose. The time series is then transformed to a wave maker signal with use of a designed transfer algorithm. Wave group propagation along the tank was recorded by 4 distant gauges and by an array of 6 densely situated gauges. This setup allows to consider the wave evolution from 10 to 85 m from the wave maker, and to obtain the wave envelope shape directly from the instrumental data. In the experiments wave groups were characterized by the steepness values up to kAcr < 0.32 and kAtr < 0.24, where k is the mean wavenumber, Acr is the crest amplitude, and Atr is the trough amplitude; and the maximum local wave slope was up to 0.34. Wave breaking phenomenon was not observed in the experiments. Different mean wave numbers and wave groups of different intensities were considered. In some cases the wave groups exhibit noticeable radiation in the course of propagation, though the groups are not dispersed fully. The effect of finite water depth is found to be significant on the wave group stability. Intense wave groups have shorter time of adjustment, what in some sense may help them to manifest their individuality clearer. The experimental tests confirm recent numerical simulations of fully nonlinear equations, where very steep stable single and interacting nonlinear wave groups were reported [1-3]. The quasi-stationary wave groups observed in numerical and laboratory experiments are strongly nonlinear analogues of the nonlinear Schrodinger envelope solitons. The results emphasize the importance of long-living nonlinear wave groups in dynamics of intense sea waves. [1] V.E. Zakharov, A.I. Dyachenko, A.O. Prokofiev, Eur. J. Mech. B / Fluids 25, 677 (2006). [2] A.I. Dyachenko, V.E. Zakharov, JETP Lett. 88, 307 (2008). [3] A.V. Slunyaev, JETP 109, 676 (2009).

Toward textbook multigrid efficiency for fully implicit resistive magnetohydrodynamics

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

Adams, Mark F.; Samtaney, Ravi, E-mail: samtaney@pppl.go; Brandt, Achi

2010-09-01

Multigrid methods can solve some classes of elliptic and parabolic equations to accuracy below the truncation error with a work-cost equivalent to a few residual calculations - so-called 'textbook' multigrid efficiency. We investigate methods to solve the system of equations that arise in time dependent magnetohydrodynamics (MHD) simulations with textbook multigrid efficiency. We apply multigrid techniques such as geometric interpolation, full approximate storage, Gauss-Seidel smoothers, and defect correction for fully implicit, nonlinear, second-order finite volume discretizations of MHD. We apply these methods to a standard resistive MHD benchmark problem, the GEM reconnection problem, and add a strong magnetic guide field,more » which is a critical characteristic of magnetically confined fusion plasmas. We show that our multigrid methods can achieve near textbook efficiency on fully implicit resistive MHD simulations.« less

Toward textbook multigrid efficiency for fully implicit resistive magnetohydrodynamics

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

Adams, Mark F.; Samtaney, Ravi; Brandt, Achi

2010-09-01

Multigrid methods can solve some classes of elliptic and parabolic equations to accuracy below the truncation error with a work-cost equivalent to a few residual calculations – so-called ‘‘textbook” multigrid efficiency. We investigate methods to solve the system of equations that arise in time dependent magnetohydrodynamics (MHD) simulations with textbook multigrid efficiency. We apply multigrid techniques such as geometric interpolation, full approximate storage, Gauss–Seidel smoothers, and defect correction for fully implicit, nonlinear, second-order finite volume discretizations of MHD. We apply these methods to a standard resistive MHD benchmark problem, the GEM reconnection problem, and add a strong magnetic guide field,more » which is a critical characteristic of magnetically confined fusion plasmas. We show that our multigrid methods can achieve near textbook efficiency on fully implicit resistive MHD simulations.« less

Toward textbook multigrid efficiency for fully implicit resistive magnetohydrodynamics

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

Adams, Mark F.; Samtaney, Ravi; Brandt, Achi

2013-12-14

Multigrid methods can solve some classes of elliptic and parabolic equations to accuracy below the truncation error with a work-cost equivalent to a few residual calculations – so-called “textbook” multigrid efficiency. We investigate methods to solve the system of equations that arise in time dependent magnetohydrodynamics (MHD) simulations with textbook multigrid efficiency. We apply multigrid techniques such as geometric interpolation, full approximate storage, Gauss-Seidel smoothers, and defect correction for fully implicit, nonlinear, second-order finite volume discretizations of MHD. We apply these methods to a standard resistive MHD benchmark problem, the GEM reconnection problem, and add a strong magnetic guide field,more » which is a critical characteristic of magnetically confined fusion plasmas. We show that our multigrid methods can achieve near textbook efficiency on fully implicit resistive MHD simulations.« less

One-dimensional nonlinear theory for rectangular helix traveling-wave tube

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

Fu, Chengfang, E-mail: fchffchf@126.com; Zhao, Bo; Yang, Yudong

A 1-D nonlinear theory of a rectangular helix traveling-wave tube (TWT) interacting with a ribbon beam is presented in this paper. The RF field is modeled by a transmission line equivalent circuit, the ribbon beam is divided into a sequence of thin rectangular electron discs with the same cross section as the beam, and the charges are assumed to be uniformly distributed over these discs. Then a method of computing the space-charge field by solving Green's Function in the Cartesian Coordinate-system is fully described. Nonlinear partial differential equations for field amplitudes and Lorentz force equations for particles are solved numericallymore » using the fourth-order Runge-Kutta technique. The tube's gain, output power, and efficiency of the above TWT are computed. The results show that increasing the cross section of the ribbon beam will improve a rectangular helix TWT's efficiency and reduce the saturated length.« less

Hierarchically Parallelized Constrained Nonlinear Solvers with Automated Substructuring

NASA Technical Reports Server (NTRS)

Padovan, Joe; Kwang, Abel

1994-01-01

This paper develops a parallelizable multilevel multiple constrained nonlinear equation solver. The substructuring process is automated to yield appropriately balanced partitioning of each succeeding level. Due to the generality of the procedure,_sequential, as well as partially and fully parallel environments can be handled. This includes both single and multiprocessor assignment per individual partition. Several benchmark examples are presented. These illustrate the robustness of the procedure as well as its capability to yield significant reductions in memory utilization and calculational effort due both to updating and inversion.

Propagation of Finite Amplitude Sound in Multiple Waveguide Modes.

NASA Astrophysics Data System (ADS)

van Doren, Thomas Walter

1993-01-01

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

Mechanical balance laws for fully nonlinear and weakly dispersive water waves

NASA Astrophysics Data System (ADS)

Kalisch, Henrik; Khorsand, Zahra; Mitsotakis, Dimitrios

2016-10-01

The Serre-Green-Naghdi system is a coupled, fully nonlinear system of dispersive evolution equations which approximates the full water wave problem. The system is known to describe accurately the wave motion at the surface of an incompressible inviscid fluid in the case when the fluid flow is irrotational and two-dimensional. The system is an extension of the well known shallow-water system to the situation where the waves are long, but not so long that dispersive effects can be neglected. In the current work, the focus is on deriving mass, momentum and energy densities and fluxes associated with the Serre-Green-Naghdi system. These quantities arise from imposing balance equations of the same asymptotic order as the evolution equations. In the case of an even bed, the conservation equations are satisfied exactly by the solutions of the Serre-Green-Naghdi system. The case of variable bathymetry is more complicated, with mass and momentum conservation satisfied exactly, and energy conservation satisfied only in a global sense. In all cases, the quantities found here reduce correctly to the corresponding counterparts in both the Boussinesq and the shallow-water scaling. One consequence of the present analysis is that the energy loss appearing in the shallow-water theory of undular bores is fully compensated by the emergence of oscillations behind the bore front. The situation is analyzed numerically by approximating solutions of the Serre-Green-Naghdi equations using a finite-element discretization coupled with an adaptive Runge-Kutta time integration scheme, and it is found that the energy is indeed conserved nearly to machine precision. As a second application, the shoaling of solitary waves on a plane beach is analyzed. It appears that the Serre-Green-Naghdi equations are capable of predicting both the shape of the free surface and the evolution of kinetic and potential energy with good accuracy in the early stages of shoaling.

Eulerian Dynamics with a Commutator Forcing

DTIC Science & Technology

2017-01-09

SIAM Review 56(4) (2014) 577–621. [Pes2015] J. Peszek. Discrete Cucker-Smale flocking model with a weakly singular weight. SIAM J. Math . Anal., to...viscosities in bounded domains. J. Math . Pures Appl. (9), 87(2):227– 235, 2007. [CV2010] L. Caffarelli, A. Vasseur, Drift diffusion equations with...Further time regularity for fully non-linear parabolic equations. Math . Res. Lett., 22(6):1749–1766, 2015. [CCTT2016] José A. Carrillo, Young-Pil

Fully nonlinear development of the most unstable goertler vortex in a three dimensional boundary layer

NASA Technical Reports Server (NTRS)

Otto, S. R.; Bassom, Andrew P.

1992-01-01

The nonlinear development is studied of the most unstable Gortler mode within a general 3-D boundary layer upon a suitably concave surface. The structure of this mode was first identified by Denier, Hall and Seddougui (1991) who demonstrated that the growth rate of this instability is O(G sup 3/5) where G is the Gortler number (taken to be large here), which is effectively a measure of the curvature of the surface. Previous researchers have described the fate of the most unstable mode within a 2-D boundary layer. Denier and Hall (1992) discussed the fully nonlinear development of the vortex in this case and showed that the nonlinearity causes a breakdown of the flow structure. The effect of crossflow and unsteadiness upon an infinitesimal unstable mode was elucidated by Bassom and Hall (1991). They demonstrated that crossflow tends to stabilize the most unstable Gortler mode, and for certain crossflow/frequency combinations the Gortler mode may be made neutrally stable. These vortex configurations naturally lend themselves to a weakly nonlinear stability analysis; work which is described in a previous article by the present author. Here we extend the ideas of Denier and Hall (1992) to the three-dimensional boundary layer problem. It is found that the numerical solution of the fully nonlinear equations is best conducted using a method which is essentially an adaption of that utilized by Denier and Hall (1992). The influence of crossflow and unsteadiness upon the breakdown of the flow is described.

On the nonlinear development of the most unstable Goertler vortex mode

NASA Technical Reports Server (NTRS)

Denier, James P.; Hall, Philip

1991-01-01

The nonlinear development of the most unstable Gortler vortex mode in boundary layer flows over curved walls is investigated. The most unstable Gortler mode is confined to a viscous wall layer of thickness O(G -1/5) and has spanwise wavelength O(G 11/5); it is, of course, most relevant to flow situations where the Gortler number G is much greater than 1. The nonlinear equations covering the evolution of this mode over an O(G -3/5) streamwise lengthscale are derived and are found to be of a fully nonparallel nature. The solution of these equations is achieved by making use of the numerical scheme used by Hall (1988) for the numerical solution of the nonlinear Gortler equations valid for O(1) Gortler numbers. Thus, the spanwise dependence of the flow is described by a Fourier expansion, whereas the streamwise and normal variations of the flow are dealt with by employing a suitable finite difference discretization of the governing equations. Our calculations demonstrate that, given a suitable initial disturbance, after a brief interval of decay, the energy in all the higher harmonics grows until a singularity is encountered at some downstream position. The structure of the flowfield as this singularity is approached suggests that the singularity is responsible for the vortices, which are initially confined to the thin viscous wall layer, moving away from the wall and into the core of the boundary layer.

Efficient numerical method for solving Cauchy problem for the Gamma equation

NASA Astrophysics Data System (ADS)

Koleva, Miglena N.

2011-12-01

In this work we consider Cauchy problem for the so called Gamma equation, derived by transforming the fully nonlinear Black-Scholes equation for option price into a quasilinear parabolic equation for the second derivative (Greek) Γ = VSS of the option price V. We develop an efficient numerical method for solving the model problem concerning different volatility terms. Using suitable change of variables the problem is transformed on finite interval, keeping original behavior of the solution at the infinity. Then we construct Picard-Newton algorithm with adaptive mesh step in time, which can be applied also in the case of non-differentiable functions. Results of numerical simulations are given.

A simple, stable, and accurate linear tetrahedral finite element for transient, nearly, and fully incompressible solid dynamics: A dynamic variational multiscale approach [A simple, stable, and accurate tetrahedral finite element for transient, nearly incompressible, linear and nonlinear elasticity: A dynamic variational multiscale approach

DOE PAGES

Scovazzi, Guglielmo; Carnes, Brian; Zeng, Xianyi; ...

2015-11-12

Here, we propose a new approach for the stabilization of linear tetrahedral finite elements in the case of nearly incompressible transient solid dynamics computations. Our method is based on a mixed formulation, in which the momentum equation is complemented by a rate equation for the evolution of the pressure field, approximated with piece-wise linear, continuous finite element functions. The pressure equation is stabilized to prevent spurious pressure oscillations in computations. Incidentally, it is also shown that many stabilized methods previously developed for the static case do not generalize easily to transient dynamics. Extensive tests in the context of linear andmore » nonlinear elasticity are used to corroborate the claim that the proposed method is robust, stable, and accurate.« less

A simple, stable, and accurate linear tetrahedral finite element for transient, nearly, and fully incompressible solid dynamics: A dynamic variational multiscale approach [A simple, stable, and accurate tetrahedral finite element for transient, nearly incompressible, linear and nonlinear elasticity: A dynamic variational multiscale approach

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

Scovazzi, Guglielmo; Carnes, Brian; Zeng, Xianyi

Here, we propose a new approach for the stabilization of linear tetrahedral finite elements in the case of nearly incompressible transient solid dynamics computations. Our method is based on a mixed formulation, in which the momentum equation is complemented by a rate equation for the evolution of the pressure field, approximated with piece-wise linear, continuous finite element functions. The pressure equation is stabilized to prevent spurious pressure oscillations in computations. Incidentally, it is also shown that many stabilized methods previously developed for the static case do not generalize easily to transient dynamics. Extensive tests in the context of linear andmore » nonlinear elasticity are used to corroborate the claim that the proposed method is robust, stable, and accurate.« less

The Numerical Studies Program for the Atmospheric General Circulation Experiment (AGCE) for Spacelab Flights

NASA Technical Reports Server (NTRS)

Fowlis, W. W. (Editor); Davis, M. H. (Editor)

1981-01-01

The atmospheric general circulation experiment (AGCE) numerical design for Spacelab flights was studied. A spherical baroclinic flow experiment which models the large scale circulations of the Earth's atmosphere was proposed. Gravity is simulated by a radial dielectric body force. The major objective of the AGCE is to study nonlinear baroclinic wave flows in spherical geometry. Numerical models must be developed which accurately predict the basic axisymmetric states and the stability of nonlinear baroclinic wave flows. A three dimensional, fully nonlinear, numerical model and the AGCE based on the complete set of equations is required. Progress in the AGCE numerical design studies program is reported.

The initial instability and finite-amplitude stability of alternate bars in straight channels

USGS Publications Warehouse

Nelson, J.M.

1990-01-01

The initial instability and fully developed stability of alternate bars in straight channels are investigated using linearized and nonlinear analyses. The fundamental instability leading to these features is identified through a linear stability analysis of the equations governing the flow and sediment transport fields. This instability is explained in terms of topographically induced steering of the flow and the associated pattern of erosion and deposition on the bed. While the linear theory is useful for examining the instability mechanism, this approach is shown to yield relatively little information about well-developed alternate bars and, specifically, the linear analysis is shown to yield poor predictions of the fully developed bar wavelength. A fully nonlinear approach is presented that permits computation of the evolution of these bed features from an initial perturbation to their fully developed morphology. This analysis indicates that there is typically substantial elongation of the bar wavelength during the evolution process, a result that is consistent with observations of bar development in flumes and natural channels. The nonlinear approach demonstrates that the eventual stability of these features is a result of the interplay between topographic steering effects, secondary flow production as a result of streamline curvature, and gravitationally induced modifications of sediment fluxes over a sloping bed. ?? 1990.

Inside black holes with synchronized hair

NASA Astrophysics Data System (ADS)

Brihaye, Yves; Herdeiro, Carlos; Radu, Eugen

2016-09-01

Recently, various examples of asymptotically flat, rotating black holes (BHs) with synchronized hair have been explicitly constructed, including Kerr BHs with scalar or Proca hair, and Myers-Perry BHs with scalar hair and a mass gap, showing there is a general mechanism at work. All these solutions have been found numerically, integrating the fully non-linear field equations of motion from the event horizon outwards. Here, we address the spacetime geometry of these solutions inside the event horizon. Firstly, we provide arguments, within linear theory, that there is no regular inner horizon for these solutions. Then, we address this question fully non-linearly, using as a tractable model five dimensional, equal spinning, Myers-Perry hairy BHs. We find that, for non-extremal solutions: (1) the inside spacetime geometry in the vicinity of the event horizon is smooth and the equations of motion can be integrated inwards; (2) before an inner horizon is reached, the spacetime curvature grows (apparently) without bound. In all cases, our results suggest the absence of a smooth Cauchy horizon, beyond which the metric can be extended, for hairy BHs with synchronized hair.

Neoclassical Simulation of Tokamak Plasmas using Continuum Gyrokinetc Code TEMPEST

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

Xu, X Q

We present gyrokinetic neoclassical simulations of tokamak plasmas with self-consistent electric field for the first time using a fully nonlinear (full-f) continuum code TEMPEST in a circular geometry. A set of gyrokinetic equations are discretized on a five dimensional computational grid in phase space. The present implementation is a Method of Lines approach where the phase-space derivatives are discretized with finite differences and implicit backwards differencing formulas are used to advance the system in time. The fully nonlinear Boltzmann model is used for electrons. The neoclassical electric field is obtained by solving gyrokinetic Poisson equation with self-consistent poloidal variation. Withmore » our 4D ({psi}, {theta}, {epsilon}, {mu}) version of the TEMPEST code we compute radial particle and heat flux, the Geodesic-Acoustic Mode (GAM), and the development of neoclassical electric field, which we compare with neoclassical theory with a Lorentz collision model. The present work provides a numerical scheme and a new capability for self-consistently studying important aspects of neoclassical transport and rotations in toroidal magnetic fusion devices.« less

The influence of dynamic inflow and torsional flexibility on rotor damping in forward flight from symbolically generated equations

NASA Technical Reports Server (NTRS)

Reddy, T. S. R.; Warmbrodt, W.

1985-01-01

The combined effects of blade torsion and dynamic inflow on the aeroelastic stability of an elastic rotor blade in forward flight are studied. The governing sets of equations of motion (fully nonlinear, linearized, and multiblade equations) used in this study are derived symbolically using a program written in FORTRAN. Stability results are presented for different structural models with and without dynamic inflow. A combination of symbolic and numerical programs at the proper stage in the derivation process makes the obtainment of final stability results an efficient and straightforward procedure.

Modulational stability of periodic solutions of the Kuramoto-Sivaskinsky equation

NASA Technical Reports Server (NTRS)

Papageorgiou, Demetrios T.; Papanicolaou, George C.; Smyrlis, Yiorgos S.

1993-01-01

We study the long-wave, modulational, stability of steady periodic solutions of the Kuramoto-Sivashinsky equation. The analysis is fully nonlinear at first, and can in principle be carried out to all orders in the small parameter, which is the ratio of the spatial period to a characteristic length of the envelope perturbations. In the linearized regime, we recover a high-order version of the results of Frisch, She, and Thual, which shows that the periodic waves are much more stable than previously expected.

The ion-acoustic soliton: A gas-dynamic viewpoint

NASA Astrophysics Data System (ADS)

McKenzie, J. F.

2002-03-01

The properties of fully nonlinear ion-acoustic solitons are investigated by interpreting conservation of total momentum as the structure equation for the proton flow in the wave. In most studies momentum conservation is regarded as the first integral of the Poisson equation for the electric potential and is interpreted as being analogous to a particle moving in a pseudo-potential well. By adopting an essentially gas-dynamic viewpoint, which emphasizes momentum conservation and the properties of the Bernoulli-type energy equations, the crucial role played by the proton sonic point becomes apparent. The relationship (implied by energy conservation) between the electron and proton speeds in the transition yields a locus—the hodograph of the system-which shows that, in the first half of the soliton, the electrons initially lag behind the protons until the charge neutral point is reached, after which they run ahead of the protons. The system reaches an equilibrium point (the center of the soliton) before the proton flow goes sonic. It follows that the critical ion-acoustic Mach number, Mc, above which smooth, continuous solitons cannot be constructed, stems from the requirement that the two equilibrium points of the structure equation coalesce at the proton sonic point of the flow. In general the range of the ion-acoustic Mach numbers, Mep, in which solitons exist, is extended beyond the classical range 1

Theory of cavitons in complex plasmas.

PubMed

Shukla, P K; Eliasson, B; Sandberg, I

2003-08-15

Nonlinear coupling between Langmuir waves with finite amplitude dispersive dust acoustic perturbations is considered. It is shown that the interaction is governed by a pair of coupled nonlinear differential equations. Numerical results reveal the formation of Langmuir envelope solitons composed of the dust density depression created by the ponderomotive force of bell-shaped Langmuir wave envelops. The associated ambipolar potential is positive. The present nonlinear theory should be able to account for the trapping of large amplitude Langmuir waves in finite amplitude dust density holes. This scenario may appear in Saturn's dense rings, and the Cassini spacecraft should be able to observe fully nonlinear cavitons, as presented herein. Furthermore, we propose that new electron-beam plasma experiments should be conducted to verify our theoretical prediction.

A Fully Nonlinear, Dynamically Consistent Numerical Model for Solid-Body Ship Motion. I. Ship Motion with Fixed Heading

NASA Technical Reports Server (NTRS)

Lin, Ray-Quing; Kuang, Weijia

2011-01-01

In this paper, we describe the details of our numerical model for simulating ship solidbody motion in a given environment. In this model, the fully nonlinear dynamical equations governing the time-varying solid-body ship motion under the forces arising from ship wave interactions are solved with given initial conditions. The net force and moment (torque) on the ship body are directly calculated via integration of the hydrodynamic pressure over the wetted surface and the buoyancy effect from the underwater volume of the actual ship hull with a hybrid finite-difference/finite-element method. Neither empirical nor free parametrization is introduced in this model, i.e. no a priori experimental data are needed for modelling. This model is benchmarked with many experiments of various ship hulls for heave, roll and pitch motion. In addition to the benchmark cases, numerical experiments are also carried out for strongly nonlinear ship motion with a fixed heading. These new cases demonstrate clearly the importance of nonlinearities in ship motion modelling.

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

Nonlinear modeling of wave-topography interactions, shear instabilities and shear induced wave breaking using vortex method

NASA Astrophysics Data System (ADS)

Guha, Anirban

2017-11-01

Theoretical studies on linear shear instabilities as well as different kinds of wave interactions often use simple velocity and/or density profiles (e.g. constant, piecewise) for obtaining good qualitative and quantitative predictions of the initial disturbances. Moreover, such simple profiles provide a minimal model to obtain a mechanistic understanding of shear instabilities. Here we have extended this minimal paradigm into nonlinear domain using vortex method. Making use of unsteady Bernoulli's equation in presence of linear shear, and extending Birkhoff-Rott equation to multiple interfaces, we have numerically simulated the interaction between multiple fully nonlinear waves. This methodology is quite general, and has allowed us to simulate diverse problems that can be essentially reduced to the minimal system with interacting waves, e.g. spilling and plunging breakers, stratified shear instabilities (Holmboe, Taylor-Caulfield, stratified Rayleigh), jet flows, and even wave-topography interaction problem like Bragg resonance. We found that the minimal models capture key nonlinear features (e.g. wave breaking features like cusp formation and roll-ups) which are observed in experiments and/or extensive simulations with smooth, realistic profiles.

Generalized Courant-Snyder Theory and Kapchinskij-Vladimirskij Distribution For High-intensity Beams In A Coupled Transverse Focusing Lattice

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

Hong QIn, Ronald Davidson

2011-07-18

The Courant-Snyder (CS) theory and the Kapchinskij-Vladimirskij (KV) distribution for high-intensity beams in a uncoupled focusing lattice are generalized to the case of coupled transverse dynamics. The envelope function is generalized to an envelope matrix, and the envelope equation becomes a matrix envelope equation with matrix operations that are non-commutative. In an uncoupled lattice, the KV distribution function, first analyzed in 1959, is the only known exact solution of the nonlinear Vlasov-Maxwell equations for high-intensity beams including self-fields in a self-consistent manner. The KV solution is generalized to high-intensity beams in a coupled transverse lattice using the generalized CS invariant.more » This solution projects to a rotating, pulsating elliptical beam in transverse configuration space. The fully self-consistent solution reduces the nonlinear Vlasov-Maxwell equations to a nonlinear matrix ordinary differential equation for the envelope matrix, which determines the geometry of the pulsating and rotating beam ellipse. These results provide us with a new theoretical tool to investigate the dynamics of high-intensity beams in a coupled transverse lattice. A strongly coupled lattice, a so-called N-rolling lattice, is studied as an example. It is found that strong coupling does not deteriorate the beam quality. Instead, the coupling induces beam rotation, and reduces beam pulsation.« less

Generalized Courant-Snyder theory and Kapchinskij-Vladimirskij distribution for high-intensity beams in a coupled transverse focusing lattice

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

Qin Hong; Department of Modern Physics, University of Science and Technology of China, Hefei, Anhui 230026; Davidson, Ronald C.

2011-05-15

The Courant-Snyder (CS) theory and the Kapchinskij-Vladimirskij (KV) distribution for high-intensity beams in an uncoupled focusing lattice are generalized to the case of coupled transverse dynamics. The envelope function is generalized to an envelope matrix, and the envelope equation becomes a matrix envelope equation with matrix operations that are noncommutative. In an uncoupled lattice, the KV distribution function, first analyzed in 1959, is the only known exact solution of the nonlinear Vlasov-Maxwell equations for high-intensity beams including self-fields in a self-consistent manner. The KV solution is generalized to high-intensity beams in a coupled transverse lattice using the generalized CS invariant.more » This solution projects to a rotating, pulsating elliptical beam in transverse configuration space. The fully self-consistent solution reduces the nonlinear Vlasov-Maxwell equations to a nonlinear matrix ordinary differential equation for the envelope matrix, which determines the geometry of the pulsating and rotating beam ellipse. These results provide us with a new theoretical tool to investigate the dynamics of high-intensity beams in a coupled transverse lattice. A strongly coupled lattice, a so-called N-rolling lattice, is studied as an example. It is found that strong coupling does not deteriorate the beam quality. Instead, the coupling induces beam rotation and reduces beam pulsation.« less

Integrability and Linear Stability of Nonlinear Waves

NASA Astrophysics Data System (ADS)

Degasperis, Antonio; Lombardo, Sara; Sommacal, Matteo

2018-03-01

It is well known that the linear stability of solutions of 1+1 partial differential equations which are integrable can be very efficiently investigated by means of spectral methods. We present here a direct construction of the eigenmodes of the linearized equation which makes use only of the associated Lax pair with no reference to spectral data and boundary conditions. This local construction is given in the general N× N matrix scheme so as to be applicable to a large class of integrable equations, including the multicomponent nonlinear Schrödinger system and the multiwave resonant interaction system. The analytical and numerical computations involved in this general approach are detailed as an example for N=3 for the particular system of two coupled nonlinear Schrödinger equations in the defocusing, focusing and mixed regimes. The instabilities of the continuous wave solutions are fully discussed in the entire parameter space of their amplitudes and wave numbers. By defining and computing the spectrum in the complex plane of the spectral variable, the eigenfrequencies are explicitly expressed. According to their topological properties, the complete classification of these spectra in the parameter space is presented and graphically displayed. The continuous wave solutions are linearly unstable for a generic choice of the coupling constants.

Fourier imaging of non-linear structure formation

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

Brandbyge, Jacob; Hannestad, Steen, E-mail: jacobb@phys.au.dk, E-mail: sth@phys.au.dk

We perform a Fourier space decomposition of the dynamics of non-linear cosmological structure formation in ΛCDM models. From N -body simulations involving only cold dark matter we calculate 3-dimensional non-linear density, velocity divergence and vorticity Fourier realizations, and use these to calculate the fully non-linear mode coupling integrals in the corresponding fluid equations. Our approach allows for a reconstruction of the amount of mode coupling between any two wavenumbers as a function of redshift. With our Fourier decomposition method we identify the transfer of power from larger to smaller scales, the stable clustering regime, the scale where vorticity becomes important,more » and the suppression of the non-linear divergence power spectrum as compared to linear theory. Our results can be used to improve and calibrate semi-analytical structure formation models.« less

Fully implicit moving mesh adaptive algorithm

NASA Astrophysics Data System (ADS)

Serazio, C.; Chacon, L.; Lapenta, G.

2006-10-01

In many problems of interest, the numerical modeler is faced with the challenge of dealing with multiple time and length scales. The former is best dealt with with fully implicit methods, which are able to step over fast frequencies to resolve the dynamical time scale of interest. The latter requires grid adaptivity for efficiency. Moving-mesh grid adaptive methods are attractive because they can be designed to minimize the numerical error for a given resolution. However, the required grid governing equations are typically very nonlinear and stiff, and of considerably difficult numerical treatment. Not surprisingly, fully coupled, implicit approaches where the grid and the physics equations are solved simultaneously are rare in the literature, and circumscribed to 1D geometries. In this study, we present a fully implicit algorithm for moving mesh methods that is feasible for multidimensional geometries. Crucial elements are the development of an effective multilevel treatment of the grid equation, and a robust, rigorous error estimator. For the latter, we explore the effectiveness of a coarse grid correction error estimator, which faithfully reproduces spatial truncation errors for conservative equations. We will show that the moving mesh approach is competitive vs. uniform grids both in accuracy (due to adaptivity) and efficiency. Results for a variety of models 1D and 2D geometries will be presented. L. Chac'on, G. Lapenta, J. Comput. Phys., 212 (2), 703 (2006) G. Lapenta, L. Chac'on, J. Comput. Phys., accepted (2006)

Gaussian solitary waves and compactons in Fermi–Pasta–Ulam lattices with Hertzian potentials

PubMed Central

James, Guillaume; Pelinovsky, Dmitry

2014-01-01

We consider a class of fully nonlinear Fermi–Pasta–Ulam (FPU) lattices, consisting of a chain of particles coupled by fractional power nonlinearities of order α>1. This class of systems incorporates a classical Hertzian model describing acoustic wave propagation in chains of touching beads in the absence of precompression. We analyse the propagation of localized waves when α is close to unity. Solutions varying slowly in space and time are searched with an appropriate scaling, and two asymptotic models of the chain of particles are derived consistently. The first one is a logarithmic Korteweg–de Vries (KdV) equation and possesses linearly orbitally stable Gaussian solitary wave solutions. The second model consists of a generalized KdV equation with Hölder-continuous fractional power nonlinearity and admits compacton solutions, i.e. solitary waves with compact support. When , we numerically establish the asymptotically Gaussian shape of exact FPU solitary waves with near-sonic speed and analytically check the pointwise convergence of compactons towards the limiting Gaussian profile. PMID:24808748

Magneto-thermo-elastokinetics of Geometrically Nonlinear Laminated Composite Plates. Part 1: Foundation of the Theory

NASA Technical Reports Server (NTRS)

Hasanyan, Davresh; Librescu, Liviu; Qin, Zhanming; Ambur, Damodar R.

2006-01-01

A fully coupled magneto-thermo-elastokinetic model of laminated composite, finitely electroconductive plates incorporating geometrical nonlinearities and subjected to a combination of magnetic and thermal fields, as well as carrying an electrical current is developed, In this context. the first-order transversely shearable plate theory in conjunction with von-Karman geometrically nonlinear strain concept is adopted. Related to the distribution of electric and magnetic field disturbances within the plate, the assumptions proposed by Ambartsumyan and his collaborators are adopted. Based on the electromagnetic equations (i.e. the ones by Faraday, Ampere, Ohm, Maxwell and Lorentz), the modified Fourier's law of heat conduction and on the elastokinetic field equations, the 3-D coupled problem is reduced to an equivalent 2- D one. The theory developed herein provides a foundation for the investigation, both analytical and numerical, of the interacting effects among the magnetic, thermal and elastic fields in multi-layered thin plates made of anisotropic materials.

The strong nonlinear interaction of Tollmien-Schlichting waves and Taylor-Goertler vortices in curved channel flow

NASA Technical Reports Server (NTRS)

Bennett, J.; Hall, P.; Smith, F. T.

1988-01-01

Viscous fluid flows with curved streamlines can support both centrifugal and viscous traveling wave instabilities. Here the interaction of these instabilities in the context of the fully developed flow in a curved channel is discussed. The viscous (Tollmein-Schlichting) instability is described asymptotically at high Reynolds numbers and it is found that it can induce a Taylor-Goertler flow even at extremely small amplitudes. In this interaction, the Tollmein-Schlichting wave can drive a vortex state with wavelength either comparable with the channel width or the wavelength of lower branch viscous modes. The nonlinear equations which describe these interactions are solved for nonlinear equilibrium states.

Nonlinear rovibrational polarization response of water vapor to ultrashort long-wave infrared pulses

NASA Astrophysics Data System (ADS)

Schuh, K.; Rosenow, P.; Kolesik, M.; Wright, E. M.; Koch, S. W.; Moloney, J. V.

2017-10-01

We study the rovibrational polarization response of water vapor using a fully correlated optical Bloch equation approach employing data from the HITRAN database. For a 10 -μ m long-wave infrared pulse the resulting linear response is negative, with a negative nonlinear response at intermediate intensities and a positive value at higher intensities. For a model atmosphere comprised of the electronic response of argon combined with the rovibrational response of water vapor this leads to a weakened positive nonlinear response at intermediate intensities. Propagation simulations using a simplified noncorrelated approach show the resultant reduction in the peak filament intensity sustained during filamentation due to the presence of the water vapor.

A mass, momentum, and energy conserving, fully implicit, scalable algorithm for the multi-dimensional, multi-species Rosenbluth-Fokker-Planck equation

NASA Astrophysics Data System (ADS)

Taitano, W. T.; Chacón, L.; Simakov, A. N.; Molvig, K.

2015-09-01

In this study, we demonstrate a fully implicit algorithm for the multi-species, multidimensional Rosenbluth-Fokker-Planck equation which is exactly mass-, momentum-, and energy-conserving, and which preserves positivity. Unlike most earlier studies, we base our development on the Rosenbluth (rather than Landau) form of the Fokker-Planck collision operator, which reduces complexity while allowing for an optimal fully implicit treatment. Our discrete conservation strategy employs nonlinear constraints that force the continuum symmetries of the collision operator to be satisfied upon discretization. We converge the resulting nonlinear system iteratively using Jacobian-free Newton-Krylov methods, effectively preconditioned with multigrid methods for efficiency. Single- and multi-species numerical examples demonstrate the advertised accuracy properties of the scheme, and the superior algorithmic performance of our approach. In particular, the discretization approach is numerically shown to be second-order accurate in time and velocity space and to exhibit manifestly positive entropy production. That is, H-theorem behavior is indicated for all the examples we have tested. The solution approach is demonstrated to scale optimally with respect to grid refinement (with CPU time growing linearly with the number of mesh points), and timestep (showing very weak dependence of CPU time with time-step size). As a result, the proposed algorithm delivers several orders-of-magnitude speedup vs. explicit algorithms.

Computation of turbulent pipe and duct flow using third order upwind scheme

NASA Technical Reports Server (NTRS)

Kawamura, T.

1986-01-01

The fully developed turbulence in a circular pipe and in a square duct is simulated directly without using turbulence models in the Navier-Stokes equations. The utilized method employs a third-order upwind scheme for the approximation to the nonlinear term and the second-order Adams-Bashforth method for the time derivative in the Navier-Stokes equation. The computational results appear to capture the large-scale turbulent structures at least qualitatively. The significance of the artificial viscosity inherent in the present scheme is discussed.

Modeling charge collection efficiency degradation in partially depleted GaAs photodiodes using the 1- and 2-carrier Hecht equations

DOE PAGES

Auden, E. C.; Vizkelethy, G.; Serkland, D. K.; ...

2017-03-24

Here, the Hecht equation can be used to model the nonlinear degradation of charge collection efficiency (CCE) in response to radiation-induced displacement damage in both fully and partially depleted GaAs photodiodes. CCE degradation is measured for laser-generated photocurrent as a function of fluence and bias in Al 0.3Ga 0.7As/GaAs/Al 0.25Ga 0.75As p-i-n photodiodes which have been irradiated with 12 MeV C and 7.5 MeV Si ions. CCE is observed to degrade more rapidly with fluence in partially depleted photodiodes than in fully depleted photodiodes. When the intrinsic GaAs layer is fully depleted, the 2-carrier Hecht equation describes CCE degradation asmore » photogenerated electrons and holes recombine at defect sites created by radiation damage in the depletion region. If the GaAs layer is partially depleted, CCE degradation is more appropriately modeled as the sum of the 2-carrier Hecht equation applied to electrons and holes generated within the depletion region and the 1-carrier Hecht equation applied to minority carriers that diffuse from the field-free (non-depleted) region into the depletion region. Enhanced CCE degradation is attributed to holes that recombine within the field-free region of the partially depleted intrinsic GaAs layer before they can diffuse into the depletion region.« less

Modeling charge collection efficiency degradation in partially depleted GaAs photodiodes using the 1- and 2-carrier Hecht equations

NASA Astrophysics Data System (ADS)

Auden, E. C.; Vizkelethy, G.; Serkland, D. K.; Bossert, D. J.; Doyle, B. L.

2017-05-01

The Hecht equation can be used to model the nonlinear degradation of charge collection efficiency (CCE) in response to radiation-induced displacement damage in both fully and partially depleted GaAs photodiodes. CCE degradation is measured for laser-generated photocurrent as a function of fluence and bias in Al0.3Ga0.7As/GaAs/Al0.25Ga0.75As p-i-n photodiodes which have been irradiated with 12 MeV C and 7.5 MeV Si ions. CCE is observed to degrade more rapidly with fluence in partially depleted photodiodes than in fully depleted photodiodes. When the intrinsic GaAs layer is fully depleted, the 2-carrier Hecht equation describes CCE degradation as photogenerated electrons and holes recombine at defect sites created by radiation damage in the depletion region. If the GaAs layer is partially depleted, CCE degradation is more appropriately modeled as the sum of the 2-carrier Hecht equation applied to electrons and holes generated within the depletion region and the 1-carrier Hecht equation applied to minority carriers that diffuse from the field-free (non-depleted) region into the depletion region. Enhanced CCE degradation is attributed to holes that recombine within the field-free region of the partially depleted intrinsic GaAs layer before they can diffuse into the depletion region.

Modeling charge collection efficiency degradation in partially depleted GaAs photodiodes using the 1- and 2-carrier Hecht equations

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

Auden, E. C.; Vizkelethy, G.; Serkland, D. K.

Here, the Hecht equation can be used to model the nonlinear degradation of charge collection efficiency (CCE) in response to radiation-induced displacement damage in both fully and partially depleted GaAs photodiodes. CCE degradation is measured for laser-generated photocurrent as a function of fluence and bias in Al 0.3Ga 0.7As/GaAs/Al 0.25Ga 0.75As p-i-n photodiodes which have been irradiated with 12 MeV C and 7.5 MeV Si ions. CCE is observed to degrade more rapidly with fluence in partially depleted photodiodes than in fully depleted photodiodes. When the intrinsic GaAs layer is fully depleted, the 2-carrier Hecht equation describes CCE degradation asmore » photogenerated electrons and holes recombine at defect sites created by radiation damage in the depletion region. If the GaAs layer is partially depleted, CCE degradation is more appropriately modeled as the sum of the 2-carrier Hecht equation applied to electrons and holes generated within the depletion region and the 1-carrier Hecht equation applied to minority carriers that diffuse from the field-free (non-depleted) region into the depletion region. Enhanced CCE degradation is attributed to holes that recombine within the field-free region of the partially depleted intrinsic GaAs layer before they can diffuse into the depletion region.« less

MOOSE: A parallel computational framework for coupled systems of nonlinear equations.

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

Derek Gaston; Chris Newman; Glen Hansen

Systems of coupled, nonlinear partial differential equations (PDEs) often arise in simulation of nuclear processes. MOOSE: Multiphysics Object Oriented Simulation Environment, a parallel computational framework targeted at the solution of such systems, is presented. As opposed to traditional data-flow oriented computational frameworks, MOOSE is instead founded on the mathematical principle of Jacobian-free Newton-Krylov (JFNK) solution methods. Utilizing the mathematical structure present in JFNK, physics expressions are modularized into `Kernels,'' allowing for rapid production of new simulation tools. In addition, systems are solved implicitly and fully coupled, employing physics based preconditioning, which provides great flexibility even with large variance in timemore » scales. A summary of the mathematics, an overview of the structure of MOOSE, and several representative solutions from applications built on the framework are presented.« less

Subcritical Kelvin-Helmholtz instability in a Hele-Shaw cell.

PubMed

Meignin, L; Gondret, P; Ruyer-Quil, C; Rabaud, M

2003-06-13

We investigate experimentally the subcritical behavior of the Kelvin-Helmholtz instability for a gas-liquid shearing flow in a Hele-Shaw cell. The subcritical curve separating the solutions of a stable plane interface and a fully saturated nonlinear wave train is determined. Experimental results are fitted by a fifth order complex Ginzburg-Landau equation whose linear coefficients are compared to theoretical ones.

Nonlinear response and bistability of driven ion acoustic waves

NASA Astrophysics Data System (ADS)

Akbari-Moghanjoughi, M.

2017-08-01

The hydrodynamic model is used to obtain a generalized pseudoforce equation through which the nonlinear response of periodically driven ion acoustic waves is studied in an electron-ion plasma with isothermal and adiabatic ion fluids. The pseudotime series, corresponding to different driving frequencies, indicates that nonlinearity effects appear more strongly for smaller frequency values. The existence of extra harmonic resonances in the nonlinear amplitude spectrum is a clear indication of the interaction of an external force with harmonic components of the nonlinear ion acoustic waves. It is shown that many plasma parameters significantly and differently affect the nonlinear resonance spectrum of ion acoustic excitations. A heuristic but accurate model for the foldover effect is used which quite satisfactorily predicts the bistability of driven plasma oscillations. It is remarked that the characteristic resonance peak of isothermal ion plasma oscillations appears at lower frequencies but is stronger compared to that of adiabatic ions. Comparison of the exact numerical results for fully nonlinear and approximate (weakly nonlinear) models indicates that a weakly nonlinear model exaggerates the hysteresis and jump phenomenon for higher values of the external force amplitude.

Wave-induced response of a floating two-dimensional body with a moonpool

PubMed Central

Fredriksen, Arnt G.; Kristiansen, Trygve; Faltinsen, Odd M.

2015-01-01

Regular wave-induced behaviour of a floating stationary two-dimensional body with a moonpool is studied. The focus is on resonant piston-mode motion in the moonpool and rigid-body motions. Dedicated two-dimensional experiments have been performed. Two numerical hybrid methods, which have previously been applied to related problems, are further developed. Both numerical methods couple potential and viscous flow. The semi-nonlinear hybrid method uses linear free-surface and body-boundary conditions. The other one uses fully nonlinear free-surface and body-boundary conditions. The harmonic polynomial cell method solves the Laplace equation in the potential flow domain, while the finite volume method solves the Navier–Stokes equations in the viscous flow domain near the body. Results from the two codes are compared with the experimental data. The nonlinear hybrid method compares well with the data, while certain discrepancies are observed for the semi-nonlinear method. In particular, the roll motion is over-predicted by the semi-nonlinear hybrid method. Error sources in the semi-nonlinear hybrid method are discussed. The moonpool strongly affects heave motions in a frequency range around the piston-mode resonance frequency of the moonpool. No resonant water motions occur in the moonpool at the piston-mode resonance frequency. Instead large moonpool motions occur at a heave natural frequency associated with small damping near the piston-mode resonance frequency. PMID:25512594

Fully nonlinear theory of transcritical shallow-water flow past topography

NASA Astrophysics Data System (ADS)

El, Gennady; Grimshaw, Roger; Smyth, Noel

2010-05-01

In this talk recent results on the generation of undular bores in one-dimensional fully nonlinear shallow-water flows past localised topographies will be presented. The description is made in the framework of the forced Su-Gardner (a.k.a. 1D Green-Naghdi) system of equations, with a primary focus on the transcritical regime when the Froude number of the oncoming flow is close to unity. A combination of the local transcritical hydraulic solution over the localized topography, which produces upstream and downstream hydraulic jumps, and unsteady undular bore solutions describing the resolution of these hydraulic jumps, is used to describe various flow regimes depending on the combination of the topography height and the Froude number. We take advantage of the recently developed modulation theory of Su-Gardner undular bores to derive the main parameters of transcritical fully nonlinear shallow-water flow, such as the leading solitary wave amplitudes for the upstream and downstream undular bores, the speeds of the undular bores edges and the drag force. Our results confirm that most of the features of the previously developed description in the framework of the uni-directional forced KdV model hold up qualitatively for finite amplitude waves, while the quantitative description can be obtained in the framework of the bi-directional forced Su-Gardner system.

A scalable, fully implicit algorithm for the reduced two-field low-β extended MHD model

DOE PAGES

Chacon, Luis; Stanier, Adam John

2016-12-01

Here, we demonstrate a scalable fully implicit algorithm for the two-field low-β extended MHD model. This reduced model describes plasma behavior in the presence of strong guide fields, and is of significant practical impact both in nature and in laboratory plasmas. The model displays strong hyperbolic behavior, as manifested by the presence of fast dispersive waves, which make a fully implicit treatment very challenging. In this study, we employ a Jacobian-free Newton–Krylov nonlinear solver, for which we propose a physics-based preconditioner that renders the linearized set of equations suitable for inversion with multigrid methods. As a result, the algorithm ismore » shown to scale both algorithmically (i.e., the iteration count is insensitive to grid refinement and timestep size) and in parallel in a weak-scaling sense, with the wall-clock time scaling weakly with the number of cores for up to 4096 cores. For a 4096 × 4096 mesh, we demonstrate a wall-clock-time speedup of ~6700 with respect to explicit algorithms. The model is validated linearly (against linear theory predictions) and nonlinearly (against fully kinetic simulations), demonstrating excellent agreement.« less

Modeling Disturbance Dynamics in Transitional and Turbulent Boundary Layers

NASA Technical Reports Server (NTRS)

Grosch, C. E.; Gatski, T. B. (Technical Monitor)

2002-01-01

The dynamics of an ensemble of linear disturbances in boundary-layer flows at various Reynolds numbers is studied through an analysis of the transport equations for the mean disturbance kinetic energy and energy dissipation rate. Effects of adverse and favorable pressure-gradients on the disturbance dynamics are also included in the analysis. Unlike the fully turbulent regime where nonlinear phase scrambling of the fluctuations affects the flow field even in proximity to the wall, the early stage transition regime fluctuations studied here are influenced across the boundary layer by the solid boundary. In addition, the dominating dynamics in the disturbance kinetic energy equation is governed by the energy production, pressure-transport and viscous diffusion - also in contrast to the fully turbulent regime. For the disturbance dissipation rate, a dynamic balance exists between the destruction and diffusion of dissipation.

Fluid preconditioning for Newton–Krylov-based, fully implicit, electrostatic particle-in-cell simulations

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

Chen, G., E-mail: gchen@lanl.gov; Chacón, L.; Leibs, C.A.

2014-02-01

A recent proof-of-principle study proposes an energy- and charge-conserving, nonlinearly implicit electrostatic particle-in-cell (PIC) algorithm in one dimension [9]. The algorithm in the reference employs an unpreconditioned Jacobian-free Newton–Krylov method, which ensures nonlinear convergence at every timestep (resolving the dynamical timescale of interest). Kinetic enslavement, which is one key component of the algorithm, not only enables fully implicit PIC as a practical approach, but also allows preconditioning the kinetic solver with a fluid approximation. This study proposes such a preconditioner, in which the linearized moment equations are closed with moments computed from particles. Effective acceleration of the linear GMRES solvemore » is demonstrated, on both uniform and non-uniform meshes. The algorithm performance is largely insensitive to the electron–ion mass ratio. Numerical experiments are performed on a 1D multi-scale ion acoustic wave test problem.« less

Spatial Soliton Interactions for Photonic Switching. Part I

DTIC Science & Technology

2000-03-07

technique , a fully vectorial, first-order nonlinear wave equation that consistently includes terms two -orders beyond the slowly-varying amplitude , slowly...by using two tunable mode-locked Er-doped fiber lasers ," in Conference on Optical Fiber Communications, OSA Technical Digest Series, vol. 4, 1994...instead, based on optical logic gates. In addition, optical logic could be used for contention resolution, real-time encryption /decryption, and other

KEEN Wave Simulations: Comparing various PIC to various fixed grid Vlasov to Phase-Space Adaptive Sparse Tiling & Effective Lagrangian (PASTEL) Techniques

NASA Astrophysics Data System (ADS)

Afeyan, Bedros; Larson, David; Shadwick, Bradley; Sydora, Richard

2017-10-01

We compare various ways of solving the Vlasov-Poisson and Vlasov-Maxwell equations on rather demanding nonlinear kinetic phenomena associated with KEEN and KEEPN waves. KEEN stands for Kinetic, Electrostatic, Electron Nonlinear, and KEEPN, for electron-positron or pair plasmas analogs. Because these self-organized phase space structures are not steady-state, or single mode, or fluid or low order moment equation limited, typical techniques with low resolution or too much noise will distort the answer too much, too soon, and fail. This will be shown via Penrose criteria triggers for instability at the formation stage as well as particle orbit statistics in fully formed KEEN waves and KEEN-KEEN and KEEN-EPW interacting states. We will argue that PASTEL is a viable alternative to traditional methods with reasonable chances of success in higher dimensions. Work supported by a Grant from AFOSR PEEP.

A Penalty Method for the Numerical Solution of Hamilton-Jacobi-Bellman (HJB) Equations in Finance

NASA Astrophysics Data System (ADS)

Witte, J. H.; Reisinger, C.

2010-09-01

We present a simple and easy to implement method for the numerical solution of a rather general class of Hamilton-Jacobi-Bellman (HJB) equations. In many cases, the considered problems have only a viscosity solution, to which, fortunately, many intuitive (e.g. finite difference based) discretisations can be shown to converge. However, especially when using fully implicit time stepping schemes with their desireable stability properties, one is still faced with the considerable task of solving the resulting nonlinear discrete system. In this paper, we introduce a penalty method which approximates the nonlinear discrete system to an order of O(1/ρ), where ρ>0 is the penalty parameter, and we show that an iterative scheme can be used to solve the penalised discrete problem in finitely many steps. We include a number of examples from mathematical finance for which the described approach yields a rigorous numerical scheme and present numerical results.

SPHERES tethered formation flight testbed: advancements in enabling NASA's SPECS mission

NASA Astrophysics Data System (ADS)

Chung, Soon-Jo; Adams, Danielle; Saenz-Otero, Alvar; Kong, Edmund; Miller, David W.; Leisawitz, David; Lorenzini, Enrico; Sell, Steve

2006-06-01

This paper reports on efforts to control a tethered formation flight spacecraft array for NASA's SPECS mission using the SPHERES test-bed developed by the MIT Space Systems Laboratory. Specifically, advances in methodology and experimental results realized since the 2005 SPIE paper are emphasized. These include a new test-bed setup with a reaction wheel assembly, a novel relative attitude measurement system using force torque sensors, and modeling of non-ideal tethers to account for tether vibration modes. The nonlinear equations of motion of multi-vehicle tethered spacecraft with elastic flexible tethers are derived from Lagrange's equations. The controllability analysis indicates that both array resizing and spin-up are fully controllable by the reaction wheels and the tether motor, thereby saving thruster fuel consumption. Based upon this analysis, linear and nonlinear controllers have been successfully implemented on the tethered SPHERES testbed, and tested at the NASA MSFC's flat floor facility using two and three SPHERES configurations.

Direct numerical simulation of the laminar-turbulent transition at hypersonic flow speeds on a supercomputer

NASA Astrophysics Data System (ADS)

Egorov, I. V.; Novikov, A. V.; Fedorov, A. V.

2017-08-01

A method for direct numerical simulation of three-dimensional unsteady disturbances leading to a laminar-turbulent transition at hypersonic flow speeds is proposed. The simulation relies on solving the full three-dimensional unsteady Navier-Stokes equations. The computational technique is intended for multiprocessor supercomputers and is based on a fully implicit monotone approximation scheme and the Newton-Raphson method for solving systems of nonlinear difference equations. This approach is used to study the development of three-dimensional unstable disturbances in a flat-plate and compression-corner boundary layers in early laminar-turbulent transition stages at the free-stream Mach number M = 5.37. The three-dimensional disturbance field is visualized in order to reveal and discuss features of the instability development at the linear and nonlinear stages. The distribution of the skin friction coefficient is used to detect laminar and transient flow regimes and determine the onset of the laminar-turbulent transition.

A Fluid Dynamic Approach to the Dust-Acoustic Soliton

NASA Astrophysics Data System (ADS)

McKenzie, J. F.; Doyle, T. B.

2002-12-01

The properties of dust-acoustic solitons are derived from a fluid dynamic viewpoint in which conservation of total momentum, combined with the Bernoulli-like energy equations for each species, yields the structure equation for the heavy (or dust) speed in the stationary wave. This fully nonlinear approach reveals the crucial role played by the heavy sonic point in limiting the collective dust-acoustic Mach number, above which solitons cannot exist. An exact solution illustrates that the cold heavy species is compressed and this implies concomitant contraints on the potential and on the flow speed of the electrons and protons in the wave.

Linear and nonlinear ion-acoustic waves in nonrelativistic quantum plasmas with arbitrary degeneracy.

PubMed

Haas, Fernando; Mahmood, Shahzad

2015-11-01

Linear and nonlinear ion-acoustic waves are studied in a fluid model for nonrelativistic, unmagnetized quantum plasma with electrons with an arbitrary degeneracy degree. The equation of state for electrons follows from a local Fermi-Dirac distribution function and applies equally well both to fully degenerate and classical, nondegenerate limits. Ions are assumed to be cold. Quantum diffraction effects through the Bohm potential are also taken into account. A general coupling parameter valid for dilute and dense plasmas is proposed. The linear dispersion relation of the ion-acoustic waves is obtained and the ion-acoustic speed is discussed for the limiting cases of extremely dense or dilute systems. In the long-wavelength limit, the results agree with quantum kinetic theory. Using the reductive perturbation method, the appropriate Korteweg-de Vries equation for weakly nonlinear solutions is obtained and the corresponding soliton propagation is analyzed. It is found that soliton hump and dip structures are formed depending on the value of the quantum parameter for the degenerate electrons, which affect the phase velocities in the dispersive medium.

Linear and nonlinear ion-acoustic waves in nonrelativistic quantum plasmas with arbitrary degeneracy

NASA Astrophysics Data System (ADS)

Haas, Fernando; Mahmood, Shahzad

2015-11-01

Linear and nonlinear ion-acoustic waves are studied in a fluid model for nonrelativistic, unmagnetized quantum plasma with electrons with an arbitrary degeneracy degree. The equation of state for electrons follows from a local Fermi-Dirac distribution function and applies equally well both to fully degenerate and classical, nondegenerate limits. Ions are assumed to be cold. Quantum diffraction effects through the Bohm potential are also taken into account. A general coupling parameter valid for dilute and dense plasmas is proposed. The linear dispersion relation of the ion-acoustic waves is obtained and the ion-acoustic speed is discussed for the limiting cases of extremely dense or dilute systems. In the long-wavelength limit, the results agree with quantum kinetic theory. Using the reductive perturbation method, the appropriate Korteweg-de Vries equation for weakly nonlinear solutions is obtained and the corresponding soliton propagation is analyzed. It is found that soliton hump and dip structures are formed depending on the value of the quantum parameter for the degenerate electrons, which affect the phase velocities in the dispersive medium.

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

Chacon, Luis; Stanier, Adam John

Here, we demonstrate a scalable fully implicit algorithm for the two-field low-β extended MHD model. This reduced model describes plasma behavior in the presence of strong guide fields, and is of significant practical impact both in nature and in laboratory plasmas. The model displays strong hyperbolic behavior, as manifested by the presence of fast dispersive waves, which make a fully implicit treatment very challenging. In this study, we employ a Jacobian-free Newton–Krylov nonlinear solver, for which we propose a physics-based preconditioner that renders the linearized set of equations suitable for inversion with multigrid methods. As a result, the algorithm ismore » shown to scale both algorithmically (i.e., the iteration count is insensitive to grid refinement and timestep size) and in parallel in a weak-scaling sense, with the wall-clock time scaling weakly with the number of cores for up to 4096 cores. For a 4096 × 4096 mesh, we demonstrate a wall-clock-time speedup of ~6700 with respect to explicit algorithms. The model is validated linearly (against linear theory predictions) and nonlinearly (against fully kinetic simulations), demonstrating excellent agreement.« less

High-order solution methods for grey discrete ordinates thermal radiative transfer

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

Maginot, Peter G., E-mail: maginot1@llnl.gov; Ragusa, Jean C., E-mail: jean.ragusa@tamu.edu; Morel, Jim E., E-mail: morel@tamu.edu

This work presents a solution methodology for solving the grey radiative transfer equations that is both spatially and temporally more accurate than the canonical radiative transfer solution technique of linear discontinuous finite element discretization in space with implicit Euler integration in time. We solve the grey radiative transfer equations by fully converging the nonlinear temperature dependence of the material specific heat, material opacities, and Planck function. The grey radiative transfer equations are discretized in space using arbitrary-order self-lumping discontinuous finite elements and integrated in time with arbitrary-order diagonally implicit Runge–Kutta time integration techniques. Iterative convergence of the radiation equation ismore » accelerated using a modified interior penalty diffusion operator to precondition the full discrete ordinates transport operator.« less

High-order solution methods for grey discrete ordinates thermal radiative transfer

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

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

This paper presents a solution methodology for solving the grey radiative transfer equations that is both spatially and temporally more accurate than the canonical radiative transfer solution technique of linear discontinuous finite element discretization in space with implicit Euler integration in time. We solve the grey radiative transfer equations by fully converging the nonlinear temperature dependence of the material specific heat, material opacities, and Planck function. The grey radiative transfer equations are discretized in space using arbitrary-order self-lumping discontinuous finite elements and integrated in time with arbitrary-order diagonally implicit Runge–Kutta time integration techniques. Iterative convergence of the radiation equation ismore » accelerated using a modified interior penalty diffusion operator to precondition the full discrete ordinates transport operator.« less

High-order solution methods for grey discrete ordinates thermal radiative transfer

DOE PAGES

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

2016-09-29

This paper presents a solution methodology for solving the grey radiative transfer equations that is both spatially and temporally more accurate than the canonical radiative transfer solution technique of linear discontinuous finite element discretization in space with implicit Euler integration in time. We solve the grey radiative transfer equations by fully converging the nonlinear temperature dependence of the material specific heat, material opacities, and Planck function. The grey radiative transfer equations are discretized in space using arbitrary-order self-lumping discontinuous finite elements and integrated in time with arbitrary-order diagonally implicit Runge–Kutta time integration techniques. Iterative convergence of the radiation equation ismore » accelerated using a modified interior penalty diffusion operator to precondition the full discrete ordinates transport operator.« less

Implicit filtered P{sub N} for high-energy density thermal radiation transport using discontinuous Galerkin finite elements

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

Laboure, Vincent M., E-mail: vincent.laboure@tamu.edu; McClarren, Ryan G., E-mail: rgm@tamu.edu; Hauck, Cory D., E-mail: hauckc@ornl.gov

2016-09-15

In this work, we provide a fully-implicit implementation of the time-dependent, filtered spherical harmonics (FP{sub N}) equations for non-linear, thermal radiative transfer. We investigate local filtering strategies and analyze the effect of the filter on the conditioning of the system, showing in particular that the filter improves the convergence properties of the iterative solver. We also investigate numerically the rigorous error estimates derived in the linear setting, to determine whether they hold also for the non-linear case. Finally, we simulate a standard test problem on an unstructured mesh and make comparisons with implicit Monte Carlo (IMC) calculations.

Attractors of three-dimensional fast-rotating Navier-Stokes equations

NASA Astrophysics Data System (ADS)

Trahe, Markus

The three-dimensional (3-D) rotating Navier-Stokes equations describe the dynamics of rotating, incompressible, viscous fluids. In this work, they are considered with smooth, time-independent forces and the original statements implied by the classical "Taylor-Proudman Theorem" of geophysics are rigorously proved. It is shown that fully developed turbulence of 3-D fast-rotating fluids is essentially characterized by turbulence of two-dimensional (2-D) fluids in terms of numbers of degrees of freedom. In this context, the 3-D nonlinear "resonant limit equations", which arise in a non-linear averaging process as the rotation frequency O → infinity, are studied and optimal (2-D-type) upper bounds for fractal box and Hausdorff dimensions of the global attractor as well as upper bounds for box dimensions of exponential attractors are determined. Then, the convergence of exponential attractors for the full 3-D rotating Navier-Stokes equations to exponential attractors for the resonant limit equations as O → infinity in the sense of full Hausdorff-metric distances is established. This provides upper and lower semi-continuity of exponential attractors with respect to the rotation frequency and implies that the number of degrees of freedom (attractor dimension) of 3-D fast-rotating fluids is close to that of 2-D fluids. Finally, the algebraic-geometric structure of the Poincare curves, which control the resonances and small divisor estimates for partial differential equations, is further investigated; the 3-D nonlinear limit resonant operators are characterized by three-wave interactions governed by these curves. A new canonical transformation between those curves is constructed; with far-reaching consequences on the density of the latter.

Superposition of elliptic functions as solutions for a large number of nonlinear equations

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

Khare, Avinash; Saxena, Avadh

2014-03-15

For a large number of nonlinear equations, both discrete and continuum, we demonstrate a kind of linear superposition. We show that whenever a nonlinear equation admits solutions in terms of both Jacobi elliptic functions cn(x, m) and dn(x, m) with modulus m, then it also admits solutions in terms of their sum as well as difference. We have checked this in the case of several nonlinear equations such as the nonlinear Schrödinger equation, MKdV, a mixed KdV-MKdV system, a mixed quadratic-cubic nonlinear Schrödinger equation, the Ablowitz-Ladik equation, the saturable nonlinear Schrödinger equation, λϕ{sup 4}, the discrete MKdV as well asmore » for several coupled field equations. Further, for a large number of nonlinear equations, we show that whenever a nonlinear equation admits a periodic solution in terms of dn{sup 2}(x, m), it also admits solutions in terms of dn {sup 2}(x,m)±√(m) cn (x,m) dn (x,m), even though cn(x, m)dn(x, m) is not a solution of these nonlinear equations. Finally, we also obtain superposed solutions of various forms for several coupled nonlinear equations.« less

Algebraic approach to solve ttbar dilepton equations

NASA Astrophysics Data System (ADS)

Sonnenschein, Lars

2006-01-01

The set of non-linear equations describing the Standard Model kinematics of the top quark an- tiqark production system in the dilepton decay channel has at most a four-fold ambiguity due to two not fully reconstructed neutrinos. Its most precise and robust solution is of major importance for measurements of top quark properties like the top quark mass and t t spin correlations. Simple algebraic operations allow to transform the non-linear equations into a system of two polynomial equations with two unknowns. These two polynomials of multidegree eight can in turn be an- alytically reduced to one polynomial with one unknown by means of resultants. The obtained univariate polynomial is of degree sixteen and the coefficients are free of any singularity. The number of its real solutions is determined analytically by means of Sturm’s theorem, which is as well used to isolate each real solution into a unique pairwise disjoint interval. The solutions are polished by seeking the sign change of the polynomial in a given interval through binary brack- eting. Further a new Ansatz - exploiting an accidental cancelation in the process of transforming the equations - is presented. It permits to transform the initial system of equations into two poly- nomial equations with two unknowns. These two polynomials of multidegree two can be reduced to one univariate polynomial of degree four by means of resultants. The obtained quartic equation can be solved analytically. The analytical solution has singularities which can be circumvented by the algebraic approach described above.

Thin film flow along a periodically-stretched elastic beam

NASA Astrophysics Data System (ADS)

Boamah Mensah, Chris; Chini, Greg; Jensen, Oliver

2017-11-01

Motivated by an application to pulmonary alveolar micro-mechanics, a system of partial differential equations is derived that governs the motion of a thin liquid film lining both sides of an inertia-less elastic substrate. The evolution of the film mass distribution is described by invoking the usual lubrication approximation while the displacement of the substrate is determined by employing a kinematically nonlinear Euler-Bernoulli beam formulation. In the parameter regime of interest, the axial strain can be readily shown to be a linear function of arc-length specified completely by the motion of ends of the substrate. In contrast, the normal force balance on the beam yields an equation for the substrate curvature that is fully coupled to the time-dependent lubrication equation. Linear analyses of both a stationary and periodically-stretched flat substrate confirm the potential for buckling instabilities and reveal an upper bound on the dimensionless axial stiffness for which the coupled thin-film/inertial-less-beam model is well-posed. Numerical simulations of the coupled system are used to explore the nonlinear development of the buckling instabilities.

A monolithic mass tracking formulation for bubbles in incompressible flow

NASA Astrophysics Data System (ADS)

Aanjaneya, Mridul; Patkar, Saket; Fedkiw, Ronald

2013-08-01

We devise a novel method for treating bubbles in incompressible flow that relies on the conservative advection of bubble mass and an associated equation of state in order to determine pressure boundary conditions inside each bubble. We show that executing this algorithm in a traditional manner leads to stability issues similar to those seen for partitioned methods for solid-fluid coupling. Therefore, we reformulate the problem monolithically. This is accomplished by first proposing a new fully monolithic approach to coupling incompressible flow to fully nonlinear compressible flow including the effects of shocks and rarefactions, and then subsequently making a number of simplifying assumptions on the air flow removing not only the nonlinearities but also the spatial variations of both the density and the pressure. The resulting algorithm is quite robust, has been shown to converge to known solutions for test problems, and has been shown to be quite effective on more realistic problems including those with multiple bubbles, merging and pinching, etc. Notably, this approach departs from a standard two-phase incompressible flow model where the air flow preserves its volume despite potentially large forces and pressure differentials in the surrounding incompressible fluid that should change its volume. Our bubbles readily change volume according to an isothermal equation of state.

Perturbed soliton excitations of Rao-dust Alfvén waves in magnetized dusty plasmas

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

Kavitha, L., E-mail: louiskavitha@yahoo.co.in; The Abdus Salam International Centre for Theoretical Physics, Trieste; Lavanya, C.

We investigate the propagation dynamics of the perturbed soliton excitations in a three component fully ionized dusty magnetoplasma consisting of electrons, ions, and heavy charged dust particulates. We derive the governing equation of motion for the two dimensional Rao-dust magnetohydrodynamic (R-D-MHD) wave by employing the inertialess electron equation of motion, inertial ion equation of motion, the continuity equations in a plasma with immobile charged dust grains, together with the Maxwell's equations, by assuming quasi neutrality and neglecting the displacement current in Ampere's law. Furthermore, we assume the massive dust particles are practically immobile since we are interested in timescales muchmore » shorter than the dusty plasma period, thereby neglecting any damping of the modes due to the grain charge fluctuations. We invoke the reductive perturbation method to represent the governing dynamics by a perturbed cubic nonlinear Schrödinger (pCNLS) equation. We solve the pCNLS, along the lines of Kodama-Ablowitz multiple scale nonlinear perturbation technique and explored the R-D-MHD waves as solitary wave excitations in a magnetized dusty plasma. Since Alfvén waves play an important role in energy transport in driving field-aligned currents, particle acceleration and heating, solar flares, and the solar wind, this representation of R-D-MHD waves as soliton excitations may have extensive applications to study the lower part of the earth's ionosphere.« less

Wave-induced response of a floating two-dimensional body with a moonpool.

PubMed

Fredriksen, Arnt G; Kristiansen, Trygve; Faltinsen, Odd M

2015-01-28

Regular wave-induced behaviour of a floating stationary two-dimensional body with a moonpool is studied. The focus is on resonant piston-mode motion in the moonpool and rigid-body motions. Dedicated two-dimensional experiments have been performed. Two numerical hybrid methods, which have previously been applied to related problems, are further developed. Both numerical methods couple potential and viscous flow. The semi-nonlinear hybrid method uses linear free-surface and body-boundary conditions. The other one uses fully nonlinear free-surface and body-boundary conditions. The harmonic polynomial cell method solves the Laplace equation in the potential flow domain, while the finite volume method solves the Navier-Stokes equations in the viscous flow domain near the body. Results from the two codes are compared with the experimental data. The nonlinear hybrid method compares well with the data, while certain discrepancies are observed for the semi-nonlinear method. In particular, the roll motion is over-predicted by the semi-nonlinear hybrid method. Error sources in the semi-nonlinear hybrid method are discussed. The moonpool strongly affects heave motions in a frequency range around the piston-mode resonance frequency of the moonpool. No resonant water motions occur in the moonpool at the piston-mode resonance frequency. Instead large moonpool motions occur at a heave natural frequency associated with small damping near the piston-mode resonance frequency. © 2014 The Author(s) Published by the Royal Society. All rights reserved.

Solving Fluid Structure Interaction Problems with an Immersed Boundary Method

NASA Technical Reports Server (NTRS)

Barad, Michael F.; Brehm, Christoph; Kiris, Cetin C.

2016-01-01

An immersed boundary method for the compressible Navier-Stokes equations can be used for moving boundary problems as well as fully coupled fluid-structure interaction is presented. The underlying Cartesian immersed boundary method of the Launch Ascent and Vehicle Aerodynamics (LAVA) framework, based on the locally stabilized immersed boundary method previously presented by the authors, is extended to account for unsteady boundary motion and coupled to linear and geometrically nonlinear structural finite element solvers. The approach is validated for moving boundary problems with prescribed body motion and fully coupled fluid structure interaction problems. Keywords: Immersed Boundary Method, Higher-Order Finite Difference Method, Fluid Structure Interaction.

A fully implicit finite element method for bidomain models of cardiac electromechanics

PubMed Central

Dal, Hüsnü; Göktepe, Serdar; Kaliske, Michael; Kuhl, Ellen

2012-01-01

We propose a novel, monolithic, and unconditionally stable finite element algorithm for the bidomain-based approach to cardiac electromechanics. We introduce the transmembrane potential, the extracellular potential, and the displacement field as independent variables, and extend the common two-field bidomain formulation of electrophysiology to a three-field formulation of electromechanics. The intrinsic coupling arises from both excitation-induced contraction of cardiac cells and the deformation-induced generation of intra-cellular currents. The coupled reaction-diffusion equations of the electrical problem and the momentum balance of the mechanical problem are recast into their weak forms through a conventional isoparametric Galerkin approach. As a novel aspect, we propose a monolithic approach to solve the governing equations of excitation-contraction coupling in a fully coupled, implicit sense. We demonstrate the consistent linearization of the resulting set of non-linear residual equations. To assess the algorithmic performance, we illustrate characteristic features by means of representative three-dimensional initial-boundary value problems. The proposed algorithm may open new avenues to patient specific therapy design by circumventing stability and convergence issues inherent to conventional staggered solution schemes. PMID:23175588

An efficient solution technique for shockwave-boundary layer interactions with flow separation and slot suction effects

NASA Technical Reports Server (NTRS)

Edwards, Jack R.; Mcrae, D. Scott

1991-01-01

An efficient method for computing two-dimensional compressible Navier-Stokes flow fields is presented. The solution algorithm is a fully-implicit approximate factorization technique based on an unsymmetric line Gauss-Seidel splitting of the equation system Jacobian matrix. Convergence characteristics are improved by the addition of acceleration techniques based on Shamanskii's method for nonlinear equations and Broyden's quasi-Newton update. Characteristic-based differencing of the equations is provided by means of Van Leer's flux vector splitting. In this investigation, emphasis is placed on the fast and accurate computation of shock-wave-boundary layer interactions with and without slot suction effects. In the latter context, a set of numerical boundary conditions for simulating the transpiration flow in an open slot is devised. Both laminar and turbulent cases are considered, with turbulent closure provided by a modified Cebeci-Smith algebraic model. Comparisons with computational and experimental data sets are presented for a variety of interactions, and a fully-coupled simulation of a plenum chamber/inlet flowfield with shock interaction and suction is also shown and discussed.

A unifying framework for ghost-free Lorentz-invariant Lagrangian field theories

NASA Astrophysics Data System (ADS)

Li, Wenliang

2018-04-01

We propose a framework for Lorentz-invariant Lagrangian field theories where Ostrogradsky's scalar ghosts could be absent. A key ingredient is the generalized Kronecker delta. The general Lagrangians are reformulated in the language of differential forms. The absence of higher order equations of motion for the scalar modes stems from the basic fact that every exact form is closed. The well-established Lagrangian theories for spin-0, spin-1, p-form, spin-2 fields have natural formulations in this framework. We also propose novel building blocks for Lagrangian field theories. Some of them are novel nonlinear derivative terms for spin-2 fields. It is nontrivial that Ostrogradsky's scalar ghosts are absent in these fully nonlinear theories.

An efficient computational method for the approximate solution of nonlinear Lane-Emden type equations arising in astrophysics

NASA Astrophysics Data System (ADS)

Singh, Harendra

2018-04-01

The key purpose of this article is to introduce an efficient computational method for the approximate solution of the homogeneous as well as non-homogeneous nonlinear Lane-Emden type equations. Using proposed computational method given nonlinear equation is converted into a set of nonlinear algebraic equations whose solution gives the approximate solution to the Lane-Emden type equation. Various nonlinear cases of Lane-Emden type equations like standard Lane-Emden equation, the isothermal gas spheres equation and white-dwarf equation are discussed. Results are compared with some well-known numerical methods and it is observed that our results are more accurate.

Nonlinear unitary quantum collapse model with self-generated noise

NASA Astrophysics Data System (ADS)

Geszti, Tamás

2018-04-01

Collapse models including some external noise of unknown origin are routinely used to describe phenomena on the quantum-classical border; in particular, quantum measurement. Although containing nonlinear dynamics and thereby exposed to the possibility of superluminal signaling in individual events, such models are widely accepted on the basis of fully reproducing the non-signaling statistical predictions of quantum mechanics. Here we present a deterministic nonlinear model without any external noise, in which randomness—instead of being universally present—emerges in the measurement process, from deterministic irregular dynamics of the detectors. The treatment is based on a minimally nonlinear von Neumann equation for a Stern–Gerlach or Bell-type measuring setup, containing coordinate and momentum operators in a self-adjoint skew-symmetric, split scalar product structure over the configuration space. The microscopic states of the detectors act as a nonlocal set of hidden parameters, controlling individual outcomes. The model is shown to display pumping of weights between setup-defined basis states, with a single winner randomly selected and the rest collapsing to zero. Environmental decoherence has no role in the scenario. Through stochastic modelling, based on Pearle’s ‘gambler’s ruin’ scheme, outcome probabilities are shown to obey Born’s rule under a no-drift or ‘fair-game’ condition. This fully reproduces quantum statistical predictions, implying that the proposed non-linear deterministic model satisfies the non-signaling requirement. Our treatment is still vulnerable to hidden signaling in individual events, which remains to be handled by future research.

Modeling of Nonlinear Hydrodynamics of the Coastal Areas of the Black Sea by the Chain of the Proprietary and Open Source Models

NASA Astrophysics Data System (ADS)

Kantardgi, Igor; Zheleznyak, Mark; Demchenko, Raisa; Dykyi, Pavlo; Kivva, Sergei; Kolomiets, Pavlo; Sorokin, Maxim

2014-05-01

The nearshore hydrodynamic fields are produced by the nonlinear interactions of the shoaling waves of different time scales and currents. To simulate the wind wave and swells propagated to the coasts, wave generated near shore currents, nonlinear-dispersive wave transformation and wave diffraction in interaction with coastal and port structure, sediment transport and coastal erosion the chains of the models should be used. The objective of this presentation is to provide an overview of the results of the application of the model chains for the assessment of the wave impacts on new construction designed at the Black Sea coasts and the impacts of these constructions on the coastal erosion/ accretion processes to demonstrate needs for further development of the nonlinear models for the coastal engineering applications. The open source models Wave Watch III and SWAN has been used to simulate wave statistics of the dedicated areas of the Black Sea in high resolution to calculated the statistical parameters of the extreme wave approaching coastal zone construction in accordance with coastal engineering standards. As the main tool for the costal hydrodynamic simulations the modeling system COASTOX-MORPHO has been used, that includes the following models. HWAVE -code based on hyperbolic version of mild slope equations., HWAVE-S - spectral version of HWAVE., BOUSS-FNL - fully nonlinear system of Boussinesq equations for simulation wave nonlinear -dispersive wave transformation in coastal areas. COASTOX-CUR - the code provided the numerical solution of the Nonlinear Shallow Water Equations (NLSWE) by finite-volume methods on the unstructured grid describing the long wave transformation in the coastal zone with the efficient drying -wetting algorithms to simulate the inundation of the coastal areas including tsunami wave runup. Coastox -Cur equations with the radiation stress term calculated via near shore wave fields simulate the wave generated nearhore currents. COASTOX-SED - the module of the simulation of the sediment transport in which the suspended sediments are simulated on the basis of the solution of 2-D advection -diffusion equation and the bottom sediment transport calculations are provided the basis of a library of the most popular semi-empirical formulas. MORPH - the module of the simulation of the morphological transformation of coastal zone based on the mass balance equation, on the basis of the sediment fluxes, calculated in the SED module. MORPH management submodel is responsible for the execution of the model chain "waves- current- sediments - morphodynamics- waves". The open source model SWASH has been used to simulate nonlinear resonance phenomena in coastal waters. The model chain was applied to simulate the potential impact of the designed shore protection structures at the Sochi Olympic Park on coastal morphodynamics, the wave parameters and nonlinear oscillations in the new ports designed in Gelenddjik and Taman at North-East coast of the Black Sea. The modeling results are compared with the results of the physical modeling in the hydraulic flumes of Moscow University of Civil Engineering.

Modified Method of Simplest Equation Applied to the Nonlinear Schrödinger Equation

NASA Astrophysics Data System (ADS)

Vitanov, Nikolay K.; Dimitrova, Zlatinka I.

2018-03-01

We consider an extension of the methodology of the modified method of simplest equation to the case of use of two simplest equations. The extended methodology is applied for obtaining exact solutions of model nonlinear partial differential equations for deep water waves: the nonlinear Schrödinger equation. It is shown that the methodology works also for other equations of the nonlinear Schrödinger kind.

Proposed solution methodology for the dynamically coupled nonlinear geared rotor mechanics equations

NASA Technical Reports Server (NTRS)

Mitchell, L. D.; David, J. W.

1983-01-01

The equations which describe the three-dimensional motion of an unbalanced rigid disk in a shaft system are nonlinear and contain dynamic-coupling terms. Traditionally, investigators have used an order analysis to justify ignoring the nonlinear terms in the equations of motion, producing a set of linear equations. This paper will show that, when gears are included in such a rotor system, the nonlinear dynamic-coupling terms are potentially as large as the linear terms. Because of this, one must attempt to solve the nonlinear rotor mechanics equations. A solution methodology is investigated to obtain approximate steady-state solutions to these equations. As an example of the use of the technique, a simpler set of equations is solved and the results compared to numerical simulations. These equations represent the forced, steady-state response of a spring-supported pendulum. These equations were chosen because they contain the type of nonlinear terms found in the dynamically-coupled nonlinear rotor equations. The numerical simulations indicate this method is reasonably accurate even when the nonlinearities are large.

Consistent nonlinear deterministic and stochastic evolution equations for deep to shallow water wave shoaling

NASA Astrophysics Data System (ADS)

Vrecica, Teodor; Toledo, Yaron

2015-04-01

One-dimensional deterministic and stochastic evolution equations are derived for the dispersive nonlinear waves while taking dissipation of energy into account. The deterministic nonlinear evolution equations are formulated using operational calculus by following the approach of Bredmose et al. (2005). Their formulation is extended to include the linear and nonlinear effects of wave dissipation due to friction and breaking. The resulting equation set describes the linear evolution of the velocity potential for each wave harmonic coupled by quadratic nonlinear terms. These terms describe the nonlinear interactions between triads of waves, which represent the leading-order nonlinear effects in the near-shore region. The equations are translated to the amplitudes of the surface elevation by using the approach of Agnon and Sheremet (1997) with the correction of Eldeberky and Madsen (1999). The only current possibility for calculating the surface gravity wave field over large domains is by using stochastic wave evolution models. Hence, the above deterministic model is formulated as a stochastic one using the method of Agnon and Sheremet (1997) with two types of stochastic closure relations (Benney and Saffman's, 1966, and Hollway's, 1980). These formulations cannot be applied to the common wave forecasting models without further manipulation, as they include a non-local wave shoaling coefficients (i.e., ones that require integration along the wave rays). Therefore, a localization method was applied (see Stiassnie and Drimer, 2006, and Toledo and Agnon, 2012). This process essentially extracts the local terms that constitute the mean nonlinear energy transfer while discarding the remaining oscillatory terms, which transfer energy back and forth. One of the main findings of this work is the understanding that the approximated non-local coefficients behave in two essentially different manners. In intermediate water depths these coefficients indeed consist of rapidly oscillating terms, but as the water depth becomes shallow they change to an exponential growth (or decay) behavior. Hence, the formerly used localization technique cannot be justified for the shallow water region. A new formulation is devised for the localization in shallow water, it approximates the nonlinear non-local shoaling coefficient in shallow water and matches it to the one fitting to the intermediate water region. This allows the model behavior to be consistent from deep water to intermediate depths and up to the shallow water regime. Various simulations of the model were performed for the cases of intermediate, and shallow water, overall the model was found to give good results in both shallow and intermediate water depths. The essential difference between the shallow and intermediate nonlinear shoaling physics is explained via the dominating class III Bragg resonances phenomenon. By inspecting the resonance conditions and the nature of the dispersion relation, it is shown that unlike in the intermediate water regime, in shallow water depths the formation of resonant interactions is possible without taking into account bottom components. References Agnon, Y. & Sheremet, A. 1997 Stochastic nonlinear shoaling of directional spectra. J. Fluid Mech. 345, 79-99. Benney, D. J. & Saffman, P. G. 1966 Nonlinear interactions of random waves. Proc. R. Soc. Lond. A 289, 301-321. Bredmose, H., Agnon, Y., Madsen, P.A. & Schaffer, H.A. 2005 Wave transformation models with exact second-order transfer. European J. of Mech. - B/Fluids 24 (6), 659-682. Eldeberky, Y. & Madsen, P. A. 1999 Deterministic and stochastic evolution equations for fully dispersive and weakly nonlinear waves. Coastal Engineering 38, 1-24. Kaihatu, J. M. & Kirby, J. T. 1995 Nonlinear transformation of waves in infinite water depth. Phys. Fluids 8, 175-188. Holloway, G. 1980 Oceanic internal waves are not weak waves. J. Phys. Oceanogr. 10, 906-914. Stiassnie, M. & Drimer, N. 2006 Prediction of long forcing waves for harbor agitation studies. J. of waterways, port, coastal and ocean engineering 132(3), 166-171. Toledo, Y. & Agnon, Y. 2012 Stochastic evolution equations with localized nonlinear shoaling coefficients. European J. of Mech. - B/Fluids 34, 13-18.

Implicit solution of Navier-Stokes equations on staggered curvilinear grids using a Newton-Krylov method with a novel analytical Jacobian.

NASA Astrophysics Data System (ADS)

Borazjani, Iman; Asgharzadeh, Hafez

2015-11-01

Flow simulations involving complex geometries and moving boundaries suffer from time-step size restriction and low convergence rates with explicit and semi-implicit schemes. Implicit schemes can be used to overcome these restrictions. However, implementing implicit solver for nonlinear equations including Navier-Stokes is not straightforward. Newton-Krylov subspace methods (NKMs) are one of the most advanced iterative methods to solve non-linear equations such as implicit descritization of the Navier-Stokes equation. The efficiency of NKMs massively depends on the Jacobian formation method, e.g., automatic differentiation is very expensive, and matrix-free methods slow down as the mesh is refined. Analytical Jacobian is inexpensive method, but derivation of analytical Jacobian for Navier-Stokes equation on staggered grid is challenging. The NKM with a novel analytical Jacobian was developed and validated against Taylor-Green vortex and pulsatile flow in a 90 degree bend. The developed method successfully handled the complex geometries such as an intracranial aneurysm with multiple overset grids, and immersed boundaries. It is shown that the NKM with an analytical Jacobian is 3 to 25 times faster than the fixed-point implicit Runge-Kutta method, and more than 100 times faster than automatic differentiation depending on the grid (size) and the flow problem. The developed methods are fully parallelized with parallel efficiency of 80-90% on the problems tested.

Characteristics of steady vibration in a rotating hub-beam system

NASA Astrophysics Data System (ADS)

Zhao, Zhen; Liu, Caishan; Ma, Wei

2016-02-01

A rotating beam features a puzzling character in which its frequencies and modal shapes may vary with the hub's inertia and its rotating speed. To highlight the essential nature behind the vibration phenomena, we analyze the steady vibration of a rotating Euler-Bernoulli beam with a quasi-steady-state stretch. Newton's law is used to derive the equations governing the beam's elastic motion and the hub's rotation. A combination of these equations results in a nonlinear partial differential equation (PDE) that fully reflects the mutual interaction between the two kinds of motion. Via the Fourier series expansion within a finite interval of time, we reduce the PDE into an infinite system of a nonlinear ordinary differential equation (ODE) in spatial domain. We further nondimensionalize the ODE and discretize it via a difference method. The frequencies and modal shapes of a general rotating beam are then determined numerically. For a low-speed beam where the ignorance of geometric stiffening is feasible, the beam's vibration characteristics are solved analytically. We validate our numerical method and the analytical solutions by comparing with either the past experiments or the past numerical findings reported in existing literature. Finally, systematic simulations are performed to demonstrate how the beam's eigenfrequencies vary with the hub's inertia and rotating speed.

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).

Constructing general partial differential equations using polynomial and neural networks.

PubMed

Zjavka, Ladislav; Pedrycz, Witold

2016-01-01

Sum fraction terms can approximate multi-variable functions on the basis of discrete observations, replacing a partial differential equation definition with polynomial elementary data relation descriptions. Artificial neural networks commonly transform the weighted sum of inputs to describe overall similarity relationships of trained and new testing input patterns. Differential polynomial neural networks form a new class of neural networks, which construct and solve an unknown general partial differential equation of a function of interest with selected substitution relative terms using non-linear multi-variable composite polynomials. The layers of the network generate simple and composite relative substitution terms whose convergent series combinations can describe partial dependent derivative changes of the input variables. This regression is based on trained generalized partial derivative data relations, decomposed into a multi-layer polynomial network structure. The sigmoidal function, commonly used as a nonlinear activation of artificial neurons, may transform some polynomial items together with the parameters with the aim to improve the polynomial derivative term series ability to approximate complicated periodic functions, as simple low order polynomials are not able to fully make up for the complete cycles. The similarity analysis facilitates substitutions for differential equations or can form dimensional units from data samples to describe real-world problems. Copyright © 2015 Elsevier Ltd. All rights reserved.

Multi-component Wronskian solution to the Kadomtsev-Petviashvili equation

NASA Astrophysics Data System (ADS)

Xu, Tao; Sun, Fu-Wei; Zhang, Yi; Li, Juan

2014-01-01

It is known that the Kadomtsev-Petviashvili (KP) equation can be decomposed into the first two members of the coupled Ablowitz-Kaup-Newell-Segur (AKNS) hierarchy by the binary non-linearization of Lax pairs. In this paper, we construct the N-th iterated Darboux transformation (DT) for the second- and third-order m-coupled AKNS systems. By using together the N-th iterated DT and Cramer's rule, we find that the KPII equation has the unreduced multi-component Wronskian solution and the KPI equation admits a reduced multi-component Wronskian solution. In particular, based on the unreduced and reduced two-component Wronskians, we obtain two families of fully-resonant line-soliton solutions which contain arbitrary numbers of asymptotic solitons as y → ∓∞ to the KPII equation, and the ordinary N-soliton solution to the KPI equation. In addition, we find that the KPI line solitons propagating in parallel can exhibit the bound state at the moment of collision.

Green-Naghdi dynamics of surface wind waves in finite depth

NASA Astrophysics Data System (ADS)

Manna, M. A.; Latifi, A.; Kraenkel, R. A.

2018-04-01

The Miles’ quasi laminar theory of waves generation by wind in finite depth h is presented. In this context, the fully nonlinear Green-Naghdi model equation is derived for the first time. This model equation is obtained by the non perturbative Green-Naghdi approach, coupling a nonlinear evolution of water waves with the atmospheric dynamics which works as in the classic Miles’ theory. A depth-dependent and wind-dependent wave growth γ is drawn from the dispersion relation of the coupled Green-Naghdi model with the atmospheric dynamics. Different values of the dimensionless water depth parameter δ = gh/U 1, with g the gravity and U 1 a characteristic wind velocity, produce two families of growth rate γ in function of the dimensionless theoretical wave-age c 0: a family of γ with h constant and U 1 variable and another family of γ with U 1 constant and h variable. The allowed minimum and maximum values of γ in this model are exhibited.

Nonlinear water waves generated by impulsive motion of submerged obstacle

NASA Astrophysics Data System (ADS)

Makarenko, N.; Kostikov, V.

2012-04-01

The fully nonlinear problem on generation of unsteady water waves by impulsively moving obstacle is studied analytically. The method involves the reduction of basic Euler equations to the integral-differential system for the wave elevation together with normal and tangential fluid velocities at the free surface. Exact model equations are derived in explicit form when the isolated obstacle is presented by totally submerged circular- or elliptic cylinder. Small-time asymptotic solution is constructed for the cylinder which starts moving with constant acceleration from rest. It is demonstrated that the leading-order solution terms describe several wave regimes such as the formation of non-stationary splash jets by vertical rising or vertical submersion of the obstacle, as well as the generation of diverging waves by horizontal- and combined motion of the obstacle under free surface. This work was supported by RFBR (grant No 10-01-00447) and by Research Program of the Russian Government (grant No 11.G34.31.0035).

Analysis of Transition-Sensitized Turbulent Transport Equations

NASA Technical Reports Server (NTRS)

Rumsey, Christopher L.; Thacker, William D.; Gatski, Thomas B.; Grosch, Chester E,

2005-01-01

The dynamics of an ensemble of linear disturbances in boundary-layer flows at various Reynolds numbers is studied through an analysis of the transport equations for the mean disturbance kinetic energy and energy dissipation rate. Effects of adverse and favorable pressure-gradients on the disturbance dynamics are also included in the analysis Unlike the fully turbulent regime where nonlinear phase scrambling of the fluctuations affects the flow field even in proximity to the wall, the early stage transition regime fluctuations studied here are influenced cross the boundary layer by the solid boundary. The dominating dynamics in the disturbance kinetic energy and dissipation rate equations are described. These results are then used to formulate transition-sensitized turbulent transport equations, which are solved in a two-step process and applied to zero-pressure-gradient flow over a flat plate. Computed results are in good agreement with experimental data.

Free Convection Nanofluid Flow in the Stagnation-Point Region of a Three-Dimensional Body

PubMed Central

Farooq, Umer

2014-01-01

Analytical results are presented for a steady three-dimensional free convection flow in the stagnation point region over a general curved isothermal surface placed in a nanofluid. The momentum equations in x- and y-directions, energy balance equation, and nanoparticle concentration equation are reduced to a set of four fully coupled nonlinear differential equations under appropriate similarity transformations. The well known technique optimal homotopy analysis method (OHAM) is used to obtain the exact solution explicitly, whose convergence is then checked in detail. Besides, the effects of the physical parameters, such as the Lewis number, the Brownian motion parameter, the thermophoresis parameter, and the buoyancy ratio on the profiles of velocities, temperature, and concentration, are studied and discussed. Furthermore the local skin friction coefficients in x- and y-directions, the local Nusselt number, and the local Sherwood number are examined for various values of the physical parameters. PMID:25114954

A Survey of Ship Motion Reduction Devices

DTIC Science & Technology

1990-09-01

Bhattacharyya, R., Dynamics of Marine Vehicles, John Wiley & Sons, Inc. (1978) Chapter 14 - Motion Control. 24. Lewison , G. R. G. and J. M. Williams, "An...given in Miller et al." and Barr and Akudinov2 s . Lewison and Williams24 show the effectiveness of these tanks through case studies of commercial...fully understood. Rigorous derivations of the stabilized equations of motion tend to be non-linear with coefficients that are hard to define. Lewison 27

Linear stability analysis of the three-dimensional thermally-driven ocean circulation: application to interdecadal oscillations

NASA Astrophysics Data System (ADS)

Huck, Thierry; Vallis, Geoffrey K.

2001-08-01

What can we learn from performing a linear stability analysis of the large-scale ocean circulation? Can we predict from the basic state the occurrence of interdecadal oscillations, such as might be found in a forward integration of the full equations of motion? If so, do the structure and period of the linearly unstable modes resemble those found in a forward integration? We pursue here a preliminary study of these questions for a case in idealized geometry, in which the full nonlinear behavior can also be explored through forward integrations. Specifically, we perform a three-dimensional linear stability analysis of the thermally-driven circulation of the planetary geostrophic equations. We examine the resulting eigenvalues and eigenfunctions, comparing them with the structure of the interdecadal oscillations found in the fully nonlinear model in various parameter regimes. We obtain a steady state by running the time-dependent, nonlinear model to equilibrium using restoring boundary conditions on surface temperature. If the surface heat fluxes are then diagnosed, and these values applied as constant flux boundary conditions, the nonlinear model switches into a state of perpetual, finite amplitude, interdecadal oscillations. We construct a linearized version of the model by empirically evaluating the tangent linear matrix at the steady state, under both restoring and constant-flux boundary conditions. An eigen-analysis shows there are no unstable eigenmodes of the linearized model with restoring conditions. In contrast, under constant flux conditions, we find a single unstable eigenmode that shows a striking resemblance to the fully-developed oscillations in terms of three-dimensional structure, period and growth rate. The mode may be damped through either surface restoring boundary conditions or sufficiently large horizontal tracer diffusion. The success of this simple numerical method in idealized geometry suggests applications in the study of the stability of the ocean circulation in more realistic configurations, and the possibility of predicting potential oceanic modes, even weakly damped, that might be excited by stochastic atmospheric forcing or mesoscale ocean eddies.

Uncertainty Propagation for Turbulent, Compressible Flow in a Quasi-1D Nozzle Using Stochastic Methods

NASA Technical Reports Server (NTRS)

Zang, Thomas A.; Mathelin, Lionel; Hussaini, M. Yousuff; Bataille, Francoise

2003-01-01

This paper describes a fully spectral, Polynomial Chaos method for the propagation of uncertainty in numerical simulations of compressible, turbulent flow, as well as a novel stochastic collocation algorithm for the same application. The stochastic collocation method is key to the efficient use of stochastic methods on problems with complex nonlinearities, such as those associated with the turbulence model equations in compressible flow and for CFD schemes requiring solution of a Riemann problem. Both methods are applied to compressible flow in a quasi-one-dimensional nozzle. The stochastic collocation method is roughly an order of magnitude faster than the fully Galerkin Polynomial Chaos method on the inviscid problem.

Two nonlinear control schemes contrasted on a hydrodynamiclike model

NASA Technical Reports Server (NTRS)

Keefe, Laurence R.

1993-01-01

The principles of two flow control strategies, those of Huebler (Luescher and Huebler, 1989) and of Ott et al. (1990) are discussed, and the two schemes are compared for their ability to control shear flow, using fully developed and transitional solutions of the Ginzburg-Landau equation as models for such flows. It was found that the effectiveness of both methods in obtaining control of fully developed flows depended strongly on the 'distance' in state space between the uncontrolled flow and goal dynamics. There were conceptual difficulties in applying the Ott et al. method to transitional convectively unstable flows. On the other hand, the Huebler method worked well, within certain limitations, although at a large cost in energy terms.

Optimal antibunching in passive photonic devices based on coupled nonlinear resonators

NASA Astrophysics Data System (ADS)

Ferretti, S.; Savona, V.; Gerace, D.

2013-02-01

We propose the use of weakly nonlinear passive materials for prospective applications in integrated quantum photonics. It is shown that strong enhancement of native optical nonlinearities by electromagnetic field confinement in photonic crystal resonators can lead to single-photon generation only exploiting the quantum interference of two coupled modes and the effect of photon blockade under resonant coherent driving. For realistic system parameters in state of the art microcavities, the efficiency of such a single-photon source is theoretically characterized by means of the second-order correlation function at zero-time delay as the main figure of merit, where major sources of loss and decoherence are taken into account within a standard master equation treatment. These results could stimulate the realization of integrated quantum photonic devices based on non-resonant material media, fully integrable with current semiconductor technology and matching the relevant telecom band operational wavelengths, as an alternative to single-photon nonlinear devices based on cavity quantum electrodynamics with artificial atoms or single atomic-like emitters.

Verification of nonlinear particle simulation of radio frequency waves in fusion plasmas

NASA Astrophysics Data System (ADS)

Kuley, Animesh; Bao, Jian; Lin, Zhihong

2015-11-01

Nonlinear global particle simulation model has been developed in GTC to study the nonlinear interactions of radio frequency (RF) waves with plasmas in tokamak. In this model, ions are considered as fully kinetic particles using the Vlasov equation and electrons are treated as guiding centers using the drift kinetic. Boris push scheme for the ion motion has been implemented in the toroidal geometry using magnetic coordinates and successfully verified for the ion cyclotron, ion Bernstein and lower hybrid waves. The nonlinear GTC simulation of the lower hybrid wave shows that the amplitude of the electrostatic potential is oscillatory due to the trapping of resonant electrons by the electric field of the lower hybrid wave. The nonresonant parametric decay is observed an IBW sideband and an ion cyclotron quasimode (ICQM). The ICQM induces an ion perpendicular heating with a heating rate proportional to the pump wave intensity. This work is supported by PPPL subcontract number S013849-F and US Department of Energy (DOE) SciDAC GSEP Program.

Fast and high-order numerical algorithms for the solution of multidimensional nonlinear fractional Ginzburg-Landau equation

NASA Astrophysics Data System (ADS)

Mohebbi, Akbar

2018-02-01

In this paper we propose two fast and accurate numerical methods for the solution of multidimensional space fractional Ginzburg-Landau equation (FGLE). In the presented methods, to avoid solving a nonlinear system of algebraic equations and to increase the accuracy and efficiency of method, we split the complex problem into simpler sub-problems using the split-step idea. For a homogeneous FGLE, we propose a method which has fourth-order of accuracy in time component and spectral accuracy in space variable and for nonhomogeneous one, we introduce another scheme based on the Crank-Nicolson approach which has second-order of accuracy in time variable. Due to using the Fourier spectral method for fractional Laplacian operator, the resulting schemes are fully diagonal and easy to code. Numerical results are reported in terms of accuracy, computational order and CPU time to demonstrate the accuracy and efficiency of the proposed methods and to compare the results with the analytical solutions. The results show that the present methods are accurate and require low CPU time. It is illustrated that the numerical results are in good agreement with the theoretical ones.

A three-dimensional dynamical model for channeled lava flow with nonlinear rheology

NASA Astrophysics Data System (ADS)

Filippucci, Marilena; Tallarico, Andrea; Dragoni, Michele

2010-05-01

Recent laboratory studies on the rheology of lava samples from different volcanic areas have highlighted that the apparent viscosity depends on a power of the strain rate. Several authors agree in attributing this dependence to the crystal content of the sample and to temperature. Starting from these results, in this paper we studied the effect of a power law rheology on a gravity-driven lava flow. The equation of motion is nonlinear in the diffusion term, and an analytical solution does not seem to be possible. The finite-volume method has been applied to solve numerically the equation governing the fully developed laminar flow of a power law non-Newtonian fluid in an inclined rectangular channel. The convergence, the stability, and the order of approximation were tested for the Newtonian rheology case, comparing the numerical solution with the available analytical solution. Results indicate that the assumption on the rheology, whether linear or nonlinear, strongly affects the velocity and/or the thickness of the lava channel both for channels with fixed geometry and for channels with constant flow rate. Results on channels with fixed geometry are confirmed by some simulations for real lava channels. Finally, the study of the Reynolds number indicates that gravity-driven lava channel flows are always in laminar regime, except for strongly nonlinear pseudoplastic fluids with low fluid consistency and at high slopes.

Nonlinear Dynamics of Nanomechanical Resonators

NASA Astrophysics Data System (ADS)

Ramakrishnan, Subramanian; Gulak, Yuiry; Sundaram, Bala; Benaroya, Haym

2007-03-01

Nanoelectromechanical systems (NEMS) offer great promise for many applications including motion and mass sensing. Recent experimental results suggest the importance of nonlinear effects in NEMS, an issue which has not been addressed fully in theory. We report on a nonlinear extension of a recent analytical model by Armour et al [1] for the dynamics of a single-electron transistor (SET) coupled to a nanomechanical resonator. We consider the nonlinear resonator motion in both (a) the Duffing and (b) nonlinear pendulum regimes. The corresponding master equations are derived and solved numerically and we consider moment approximations as well. In the Duffing case with hardening stiffness, we observe that the resonator is damped by the SET at a significantly higher rate. In the cases of softening stiffness and the pendulum, there exist regimes where the SET adds energy to the resonator. To our knowledge, this is the first instance of a single model displaying both negative and positive resonator damping in different dynamical regimes. The implications of the results for SET sensitivity as well as for, as yet unexplained, experimental results will be discussed. 1. Armour et al. Phys.Rev.B (69) 125313 (2004).

Current structure of strongly nonlinear interfacial solitary waves

NASA Astrophysics Data System (ADS)

Semin, Sergey; Kurkina, Oxana; Kurkin, Andrey; Talipova, Tatiana; Pelinovsky, Efim; Churaev, Egor

2015-04-01

The characteristics of highly nonlinear solitary internal waves (solitons) in two-layer flow are computed within the fully nonlinear Navier-Stokes equations with use of numerical model of the Massachusetts Institute of Technology (MITgcm). The verification and adaptation of the model is based on the data from laboratory experiments [Carr & Davies, 2006]. The present paper also compares the results of our calculations with the computations performed in the framework of the fully nonlinear Bergen Ocean Model [Thiem et al, 2011]. The comparison of the computed soliton parameters with the predictions of the weakly nonlinear theory based on the Gardner equation is given. The occurrence of reverse flow in the bottom layer directly behind the soliton is confirmed in numerical simulations. The trajectories of Lagrangian particles in the internal soliton on the surface, on the interface and near the bottom are computed. The results demonstrated completely different trajectories at different depths of the model area. Thus, in the surface layer is observed the largest displacement of Lagrangian particles, which can be more than two and a half times larger than the characteristic width of the soliton. Located at the initial moment along the middle pycnocline fluid particles move along the elongated vertical loop at a distance of not more than one third of the width of the solitary wave. In the bottom layer of the fluid moves in the opposite direction of propagation of the internal wave, but under the influence of the reverse flow, when the bulk of the velocity field of the soliton ceases to influence the trajectory, it moves in the opposite direction. The magnitude of displacement of fluid particles in the bottom layer is not more than the half-width of the solitary wave. 1. Carr, M., and Davies, P.A. The motion of an internal solitary wave of depression over a fixed bottom boundary in a shallow, two-layer fluid. Phys. Fluids, 2006, vol. 18, No. 1, 1 - 10. 2. Thiem, O., Carr, M., Berntsen, J., and Davies, P.A. Numerical simulation of internal solitary wave-induced reverse flow and associated vortices in a shallow, two-layer fluid benthic boundary layer. Ocean Dynamics, 2011, vol. 61, No. 6, 857 - 872.

Nonlinear stability of non-stationary cross-flow vortices in compressible boundary layers

NASA Technical Reports Server (NTRS)

Gajjar, J. S. B.

1995-01-01

The nonlinear evolution of long wavelength non-stationary cross-flow vortices in a compressible boundary layer is investigated and the work extends that of Gajjar (1994) to flows involving multiple critical layers. The basic flow profile considered in this paper is that appropriate for a fully three-dimensional boundary layer with O(1) Mach number and with wall heating or cooling. The governing equations for the evolution of the cross-flow vortex are obtained and some special cases are discussed. One special case includes linear theory where exact analytic expressions for the growth rate of the vortices are obtained. Another special case is a generalization of the Bassom & Gajjar (1988) results for neutral waves to compressible flows. The viscous correction to the growth rate is derived and it is shown how the unsteady nonlinear critical layer structure merges with that for a Haberman type of viscous critical layer.

Fully decoupled monolithic projection method for natural convection problems

NASA Astrophysics Data System (ADS)

Pan, Xiaomin; Kim, Kyoungyoun; Lee, Changhoon; Choi, Jung-Il

2017-04-01

To solve time-dependent natural convection problems, we propose a fully decoupled monolithic projection method. The proposed method applies the Crank-Nicolson scheme in time and the second-order central finite difference in space. To obtain a non-iterative monolithic method from the fully discretized nonlinear system, we first adopt linearizations of the nonlinear convection terms and the general buoyancy term with incurring second-order errors in time. Approximate block lower-upper decompositions, along with an approximate factorization technique, are additionally employed to a global linearly coupled system, which leads to several decoupled subsystems, i.e., a fully decoupled monolithic procedure. We establish global error estimates to verify the second-order temporal accuracy of the proposed method for velocity, pressure, and temperature in terms of a discrete l2-norm. Moreover, according to the energy evolution, the proposed method is proved to be stable if the time step is less than or equal to a constant. In addition, we provide numerical simulations of two-dimensional Rayleigh-Bénard convection and periodic forced flow. The results demonstrate that the proposed method significantly mitigates the time step limitation, reduces the computational cost because only one Poisson equation is required to be solved, and preserves the second-order temporal accuracy for velocity, pressure, and temperature. Finally, the proposed method reasonably predicts a three-dimensional Rayleigh-Bénard convection for different Rayleigh numbers.

Numerical Simulations of Light Bullets, Using The Full Vector, Time Dependent, Nonlinear Maxwell Equations

NASA Technical Reports Server (NTRS)

Goorjian, Peter M.; Silberberg, Yaron; Kwak, Dochan (Technical Monitor)

1994-01-01

This paper will present results in computational nonlinear optics. An algorithm will be described that solves the full vector nonlinear Maxwell's equations exactly without the approximations that are currently made. Present methods solve a reduced scalar wave equation, namely the nonlinear Schrodinger equation, and neglect the optical carrier. Also, results will be shown of calculations of 2-D electromagnetic nonlinear waves computed by directly integrating in time the nonlinear vector Maxwell's equations. The results will include simulations of 'light bullet' like pulses. Here diffraction and dispersion will be counteracted by nonlinear effects. The time integration efficiently implements linear and nonlinear convolutions for the electric polarization, and can take into account such quantum effects as Kerr and Raman interactions. The present approach is robust and should permit modeling 2-D and 3-D optical soliton propagation, scattering, and switching directly from the full-vector Maxwell's equations.

Numerical Simulations of Light Bullets, Using The Full Vector, Time Dependent, Nonlinear Maxwell Equations

NASA Technical Reports Server (NTRS)

Goorjian, Peter M.; Silberberg, Yaron; Kwak, Dochan (Technical Monitor)

1995-01-01

This paper will present results in computational nonlinear optics. An algorithm will be described that solves the full vector nonlinear Maxwell's equations exactly without the approximations that we currently made. Present methods solve a reduced scalar wave equation, namely the nonlinear Schrodinger equation, and neglect the optical carrier. Also, results will be shown of calculations of 2-D electromagnetic nonlinear waves computed by directly integrating in time the nonlinear vector Maxwell's equations. The results will include simulations of 'light bullet' like pulses. Here diffraction and dispersion will be counteracted by nonlinear effects. The time integration efficiently implements linear and nonlinear convolutions for the electric polarization, and can take into account such quantum effects as Karr and Raman interactions. The present approach is robust and should permit modeling 2-D and 3-D optical soliton propagation, scattering, and switching directly from the full-vector Maxwell's equations.

High-accuracy power series solutions with arbitrarily large radius of convergence for the fractional nonlinear Schrödinger-type equations

NASA Astrophysics Data System (ADS)

Khawaja, U. Al; Al-Refai, M.; Shchedrin, Gavriil; Carr, Lincoln D.

2018-06-01

Fractional nonlinear differential equations present an interplay between two common and important effective descriptions used to simplify high dimensional or more complicated theories: nonlinearity and fractional derivatives. These effective descriptions thus appear commonly in physical and mathematical modeling. We present a new series method providing systematic controlled accuracy for solutions of fractional nonlinear differential equations, including the fractional nonlinear Schrödinger equation and the fractional nonlinear diffusion equation. The method relies on spatially iterative use of power series expansions. Our approach permits an arbitrarily large radius of convergence and thus solves the typical divergence problem endemic to power series approaches. In the specific case of the fractional nonlinear Schrödinger equation we find fractional generalizations of cnoidal waves of Jacobi elliptic functions as well as a fractional bright soliton. For the fractional nonlinear diffusion equation we find the combination of fractional and nonlinear effects results in a more strongly localized solution which nevertheless still exhibits power law tails, albeit at a much lower density.

Fully coupled simulation of multiple hydraulic fractures to propagate simultaneously from a perforated horizontal wellbore

NASA Astrophysics Data System (ADS)

Zeng, Qinglei; Liu, Zhanli; Wang, Tao; Gao, Yue; Zhuang, Zhuo

2018-02-01

In hydraulic fracturing process in shale rock, multiple fractures perpendicular to a horizontal wellbore are usually driven to propagate simultaneously by the pumping operation. In this paper, a numerical method is developed for the propagation of multiple hydraulic fractures (HFs) by fully coupling the deformation and fracturing of solid formation, fluid flow in fractures, fluid partitioning through a horizontal wellbore and perforation entry loss effect. The extended finite element method (XFEM) is adopted to model arbitrary growth of the fractures. Newton's iteration is proposed to solve these fully coupled nonlinear equations, which is more efficient comparing to the widely adopted fixed-point iteration in the literatures and avoids the need to impose fluid pressure boundary condition when solving flow equations. A secant iterative method based on the stress intensity factor (SIF) is proposed to capture different propagation velocities of multiple fractures. The numerical results are compared with theoretical solutions in literatures to verify the accuracy of the method. The simultaneous propagation of multiple HFs is simulated by the newly proposed algorithm. The coupled influences of propagation regime, stress interaction, wellbore pressure loss and perforation entry loss on simultaneous propagation of multiple HFs are investigated.

A numerical scheme for nonlinear Helmholtz equations with strong nonlinear optical effects.

PubMed

Xu, Zhengfu; Bao, Gang

2010-11-01

A numerical scheme is presented to solve the nonlinear Helmholtz (NLH) equation modeling second-harmonic generation (SHG) in photonic bandgap material doped with a nonlinear χ((2)) effect and the NLH equation modeling wave propagation in Kerr type gratings with a nonlinear χ((3)) effect in the one-dimensional case. Both of these nonlinear phenomena arise as a result of the combination of high electromagnetic mode density and nonlinear reaction from the medium. When the mode intensity of the incident wave is significantly strong, which makes the nonlinear effect non-negligible, numerical methods based on the linearization of the essentially nonlinear problem will become inadequate. In this work, a robust, stable numerical scheme is designed to simulate the NLH equations with strong nonlinearity.

Green functions and Langevin equations for nonlinear diffusion equations: A comment on ‘Markov processes, Hurst exponents, and nonlinear diffusion equations’ by Bassler et al.

NASA Astrophysics Data System (ADS)

Frank, T. D.

2008-02-01

We discuss two central claims made in the study by Bassler et al. [K.E. Bassler, G.H. Gunaratne, J.L. McCauley, Physica A 369 (2006) 343]. Bassler et al. claimed that Green functions and Langevin equations cannot be defined for nonlinear diffusion equations. In addition, they claimed that nonlinear diffusion equations are linear partial differential equations disguised as nonlinear ones. We review bottom-up and top-down approaches that have been used in the literature to derive Green functions for nonlinear diffusion equations and, in doing so, show that the first claim needs to be revised. We show that the second claim as well needs to be revised. To this end, we point out similarities and differences between non-autonomous linear Fokker-Planck equations and autonomous nonlinear Fokker-Planck equations. In this context, we raise the question whether Bassler et al.’s approach to financial markets is physically plausible because it necessitates the introduction of external traders and causes. Such external entities can easily be eliminated when taking self-organization principles and concepts of nonextensive thermostatistics into account and modeling financial processes by means of nonlinear Fokker-Planck equations.

A new solution procedure for a nonlinear infinite beam equation of motion

NASA Astrophysics Data System (ADS)

Jang, T. S.

2016-10-01

Our goal of this paper is of a purely theoretical question, however which would be fundamental in computational partial differential equations: Can a linear solution-structure for the equation of motion for an infinite nonlinear beam be directly manipulated for constructing its nonlinear solution? Here, the equation of motion is modeled as mathematically a fourth-order nonlinear partial differential equation. To answer the question, a pseudo-parameter is firstly introduced to modify the equation of motion. And then, an integral formalism for the modified equation is found here, being taken as a linear solution-structure. It enables us to formulate a nonlinear integral equation of second kind, equivalent to the original equation of motion. The fixed point approach, applied to the integral equation, results in proposing a new iterative solution procedure for constructing the nonlinear solution of the original beam equation of motion, which consists luckily of just the simple regular numerical integration for its iterative process; i.e., it appears to be fairly simple as well as straightforward to apply. A mathematical analysis is carried out on both natures of convergence and uniqueness of the iterative procedure by proving a contractive character of a nonlinear operator. It follows conclusively,therefore, that it would be one of the useful nonlinear strategies for integrating the equation of motion for a nonlinear infinite beam, whereby the preceding question may be answered. In addition, it may be worth noticing that the pseudo-parameter introduced here has double roles; firstly, it connects the original beam equation of motion with the integral equation, second, it is related with the convergence of the iterative method proposed here.

Obliquely propagating ion acoustic solitary structures in the presence of quantized magnetic field

NASA Astrophysics Data System (ADS)

Iqbal Shaukat, Muzzamal

2017-10-01

The effect of linear and nonlinear propagation of electrostatic waves have been studied in degenerate magnetoplasma taking into account the effect of electron trapping and finite temperature with quantizing magnetic field. The formation of solitary structures has been investigated by employing the small amplitude approximation both for fully and partially degenerate quantum plasma. It is observed that the inclusion of quantizing magnetic field significantly affects the propagation characteristics of the solitary wave. Importantly, the Zakharov-Kuznetsov equation under consideration has been found to allow the formation of compressive solitary structures only. The present investigation may be beneficial to understand the propagation of nonlinear electrostatic structures in dense astrophysical environments such as those found in white dwarfs.

Exact traveling wave solutions for system of nonlinear evolution equations.

PubMed

Khan, Kamruzzaman; Akbar, M Ali; Arnous, Ahmed H

2016-01-01

In this work, recently deduced generalized Kudryashov method is applied to the variant Boussinesq equations, and the (2 + 1)-dimensional breaking soliton equations. As a result a range of qualitative explicit exact traveling wave solutions are deduced for these equations, which motivates us to develop, in the near future, a new approach to obtain unsteady solutions of autonomous nonlinear evolution equations those arise in mathematical physics and engineering fields. It is uncomplicated to extend this method to higher-order nonlinear evolution equations in mathematical physics. And it should be possible to apply the same method to nonlinear evolution equations having more general forms of nonlinearities by utilizing the traveling wave hypothesis.

A scalable nonlinear fluid-structure interaction solver based on a Schwarz preconditioner with isogeometric unstructured coarse spaces in 3D

NASA Astrophysics Data System (ADS)

Kong, Fande; Cai, Xiao-Chuan

2017-07-01

Nonlinear fluid-structure interaction (FSI) problems on unstructured meshes in 3D appear in many applications in science and engineering, such as vibration analysis of aircrafts and patient-specific diagnosis of cardiovascular diseases. In this work, we develop a highly scalable, parallel algorithmic and software framework for FSI problems consisting of a nonlinear fluid system and a nonlinear solid system, that are coupled monolithically. The FSI system is discretized by a stabilized finite element method in space and a fully implicit backward difference scheme in time. To solve the large, sparse system of nonlinear algebraic equations at each time step, we propose an inexact Newton-Krylov method together with a multilevel, smoothed Schwarz preconditioner with isogeometric coarse meshes generated by a geometry preserving coarsening algorithm. Here "geometry" includes the boundary of the computational domain and the wet interface between the fluid and the solid. We show numerically that the proposed algorithm and implementation are highly scalable in terms of the number of linear and nonlinear iterations and the total compute time on a supercomputer with more than 10,000 processor cores for several problems with hundreds of millions of unknowns.

A scalable nonlinear fluid–structure interaction solver based on a Schwarz preconditioner with isogeometric unstructured coarse spaces in 3D

DOE PAGES

Kong, Fande; Cai, Xiao-Chuan

2017-03-24

Nonlinear fluid-structure interaction (FSI) problems on unstructured meshes in 3D appear many applications in science and engineering, such as vibration analysis of aircrafts and patient-specific diagnosis of cardiovascular diseases. In this work, we develop a highly scalable, parallel algorithmic and software framework for FSI problems consisting of a nonlinear fluid system and a nonlinear solid system, that are coupled monolithically. The FSI system is discretized by a stabilized finite element method in space and a fully implicit backward difference scheme in time. To solve the large, sparse system of nonlinear algebraic equations at each time step, we propose an inexactmore » Newton-Krylov method together with a multilevel, smoothed Schwarz preconditioner with isogeometric coarse meshes generated by a geometry preserving coarsening algorithm. Here ''geometry'' includes the boundary of the computational domain and the wet interface between the fluid and the solid. We show numerically that the proposed algorithm and implementation are highly scalable in terms of the number of linear and nonlinear iterations and the total compute time on a supercomputer with more than 10,000 processor cores for several problems with hundreds of millions of unknowns.« less

On symmetries, conservation laws and exact solutions of the nonlinear Schrödinger-Hirota equation

NASA Astrophysics Data System (ADS)

Akbulut, Arzu; Taşcan, Filiz

2018-04-01

In this paper, conservation laws and exact solution are found for nonlinear Schrödinger-Hirota equation. Conservation theorem is used for finding conservation laws. We get modified conservation laws for given equation. Modified simple equation method is used to obtain the exact solutions of the nonlinear Schrödinger-Hirota equation. It is shown that the suggested method provides a powerful mathematical instrument for solving nonlinear equations in mathematical physics and engineering.

An alternative explicit model expression equivalent to the integrated michaelis-menten equation and its application to nonlinear saturation pharmacokinetics.

PubMed

Goličnik, Marko

2011-06-01

Many pharmacodynamic processes can be described by the nonlinear saturation kinetics that are most frequently based on the hyperbolic Michaelis-Menten equation. Thus, various time-dependent solutions for drugs obeying such kinetics can be expressed in terms of the Lambert W(x)-omega function. However, unfortunately, computer programs that can perform the calculations for W(x) are not widely available. To avoid this problem, the replacement of the integrated Michaelis-Menten equation with an empiric integrated 1--exp alternative model equation was proposed recently by Keller et al. (Ther Drug Monit. 2009;31:783-785), although, as shown here, it was not necessary. Simulated concentrations of model drugs obeying Michaelis-Menten elimination kinetics were generated by two approaches: 1) calculation of time-course data based on an approximation equation W2*(x) performed using Microsoft Excel; and 2) calculation of reference time-course data based on an exact W(x) function built in to the Wolfram Mathematica. I show here that the W2*(x) function approximates the actual W(x) accurately. W2*(x) is expressed in terms of elementary mathematical functions and, consequently, it can be easily implemented using any of the widely available software. Hence, with the example of a hypothetical drug, I demonstrate here that an equation based on this approximation is far better, because it is nearly equivalent to the original solution, whereas the same characteristics cannot be fully confirmed for the 1--exp model equation. The W2*(x) equation proposed here might have an important role as a useful shortcut in optional software to estimate kinetic parameters from experimental data for drugs, and it might represent an easy and universal analytical tool for simulating and designing dosing regimens.

Fully Nonlinear Edge Gyrokinetic Simulations of Kinetic Geodesic-Acoustic Modes and Boundary Flows

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

Xu, X Q; Belli, E; Bodi, K

We present edge gyrokinetic neoclassical simulations of tokamak plasmas using the fully nonlinear (full-f) continuum code TEMPEST. A nonlinear Boltzmann model is used for the electrons. The electric field is obtained by solving the 2D gyrokinetic Poisson Equation. We demonstrate the following: (1) High harmonic resonances (n > 2) significantly enhance geodesic-acoustic mode (GAM) damping at high-q (tokamak safety factor), and are necessary to explain both the damping observed in our TEMPEST q-scans and experimental measurements of the scaling of the GAM amplitude with edge q{sub 95} in the absence of obvious evidence that there is a strong q dependencemore » of the turbulent drive and damping of the GAM. (2) The kinetic GAM exists in the edge for steep density and temperature gradients in the form of outgoing waves, its radial scale is set by the ion temperature profile, and ion temperature inhomogeneity is necessary for GAM radial propagation. (3) The development of the neoclassical electric field evolves through different phases of relaxation, including GAMs, their radial propagation, and their long-time collisional decay. (4) Natural consequences of orbits in the pedestal and scrape-off layer region in divertor geometry are substantial non-Maxwellian ion distributions and flow characteristics qualitatively like those observed in experiments.« less

New exact solutions of the Tzitzéica-type equations in non-linear optics using the expa function method

NASA Astrophysics Data System (ADS)

Hosseini, K.; Ayati, Z.; Ansari, R.

2018-04-01

One specific class of non-linear evolution equations, known as the Tzitzéica-type equations, has received great attention from a group of researchers involved in non-linear science. In this article, new exact solutions of the Tzitzéica-type equations arising in non-linear optics, including the Tzitzéica, Dodd-Bullough-Mikhailov and Tzitzéica-Dodd-Bullough equations, are obtained using the expa function method. The integration technique actually suggests a useful and reliable method to extract new exact solutions of a wide range of non-linear evolution equations.

The Nonlinear Steepest Descent Method to Long-Time Asymptotics of the Coupled Nonlinear Schrödinger Equation

NASA Astrophysics Data System (ADS)

Geng, Xianguo; Liu, Huan

2018-04-01

The Riemann-Hilbert problem for the coupled nonlinear Schrödinger equation is formulated on the basis of the corresponding 3× 3 matrix spectral problem. Using the nonlinear steepest descent method, we obtain leading-order asymptotics for the Cauchy problem of the coupled nonlinear Schrödinger equation.

Kinetic effects on Alfven wave nonlinearity. II - The modified nonlinear wave equation

NASA Technical Reports Server (NTRS)

Spangler, Steven R.

1990-01-01

A previously developed Vlasov theory is used here to study the role of resonant particle and other kinetic effects on Alfven wave nonlinearity. A hybrid fluid-Vlasov equation approach is used to obtain a modified version of the derivative nonlinear Schroedinger equation. The differences between a scalar model for the plasma pressure and a tensor model are discussed. The susceptibilty of the modified nonlinear wave equation to modulational instability is studied. The modulational instability normally associated with the derivative nonlinear Schroedinger equation will, under most circumstances, be restricted to left circularly polarized waves. The nonlocal term in the modified nonlinear wave equation engenders a new modulational instability that is independent of beta and the sense of circular polarization. This new instability may explain the occurrence of wave packet steepening for all values of the plasma beta in the vicinity of the earth's bow shock.

Numerical Simulations of Self-Focused Pulses Using the Nonlinear Maxwell Equations

NASA Technical Reports Server (NTRS)

Goorjian, Peter M.; Silberberg, Yaron; Kwak, Dochan (Technical Monitor)

1994-01-01

This paper will present results in computational nonlinear optics. An algorithm will be described that solves the full vector nonlinear Maxwell's equations exactly without the approximations that are currently made. Present methods solve a reduced scalar wave equation, namely the nonlinear Schrodinger equation, and neglect the optical carrier. Also, results will be shown of calculations of 2-D electromagnetic nonlinear waves computed by directly integrating in time the nonlinear vector Maxwell's equations. The results will include simulations of 'light bullet' like pulses. Here diffraction and dispersion will be counteracted by nonlinear effects. The time integration efficiently implements linear and nonlinear convolutions for the electric polarization, and can take into account such quantum effects as Kerr and Raman interactions. The present approach is robust and should permit modeling 2-D and 3-D optical soliton propagation, scattering, and switching directly from the full-vector Maxwell's equations. Abstract of a proposed paper for presentation at the meeting NONLINEAR OPTICS: Materials, Fundamentals, and Applications, Hyatt Regency Waikaloa, Waikaloa, Hawaii, July 24-29, 1994, Cosponsored by IEEE/Lasers and Electro-Optics Society and Optical Society of America

Breather Rogue Waves in Random Seas

NASA Astrophysics Data System (ADS)

Wang, J.; Ma, Q. W.; Yan, S.; Chabchoub, A.

2018-01-01

Rogue or freak waves are extreme wave events that have heights exceeding 8 times the standard deviation of surrounding waves and emerge, for instance, in the ocean as well as in other physical dispersive wave guides, such as in optical fibers. One effective and convenient way to model such an extreme dynamics in laboratory environments within a controlled framework as well as for short process time and length scales is provided through the breather formalism. Breathers are pulsating localized structures known to model extreme waves in several nonlinear dispersive media in which the initial underlying process is assumed to be narrow banded. On the other hand, several recent studies suggest that breathers can also persist in more complex environments, such as in random seas, beyond the attributed physical limitations. In this work, we study the robustness of the Peregrine breather (PB) embedded in Joint North Sea Wave Project (JONSWAP) configurations using fully nonlinear hydrodynamic numerical simulations in order to validate its practicalness for ocean engineering applications. We provide a specific range for both the spectral bandwidth of the dynamical process as well as the background wave steepness and, thus, quantify the applicability of the PB in modeling rogue waves in realistic oceanic conditions. Our results may motivate analogous studies in fields of physics such as optics and plasma to quantify the limitations of exact weakly nonlinear models, such as solitons and breathers, within the framework of the fully nonlinear governing equations of the corresponding medium.

Comparison of Quasi-Conservative Pressure-Based and Fully-Conservative Formulations for the Simulation of Transcritical Flows

NASA Astrophysics Data System (ADS)

Lacaze, Guilhem; Oefelein, Joseph

2016-11-01

High-pressure flows are known to be challenging to simulate due to thermodynamic non-linearities occurring in the vicinity of the pseudo-boiling line. This study investigates the origin of this issue by analyzing the behavior of thermodynamic processes at elevated pressure and low temperature. We show that under transcritical conditions, non-linearities significantly amplify numerical errors associated with construction of fluxes. These errors affect the local density and energy balances, which in turn creates pressure oscillations. For that reason, solvers based on a conservative system of equations that transport density and total energy are subject to unphysical pressure variations in gradient regions. These perturbations hinder numerical stability and degrade the accuracy of predictions. To circumvent this problem, the governing system can be reformulated to a pressure-based treatment of energy. We present comparisons between the pressure-based and fully conservative formulations using a progressive set of canonical cases, including a cryogenic turbulent mixing layer at rocket engine conditions. Department of Energy, Office of Science, Basic Energy Sciences Program.

Mathematical nonlinear optics

NASA Astrophysics Data System (ADS)

McLaughlin, David W.

1995-08-01

The principal investigator, together with a post-doctoral fellows Tetsuji Ueda and Xiao Wang, several graduate students, and colleagues, has applied the modern mathematical theory of nonlinear waves to problems in nonlinear optics and to equations directly relevant to nonlinear optics. Projects included the interaction of laser light with nematic liquid crystals and chaotic, homoclinic, small dispersive, and random behavior of solutions of the nonlinear Schroedinger equation. In project 1, the extremely strong nonlinear response of a continuous wave laser beam in a nematic liquid crystal medium has produced striking undulation and filamentation of the laser beam which has been observed experimentally and explained theoretically. In project 2, qualitative properties of the nonlinear Schroedinger equation (which is the fundamental equation for nonlinear optics) have been identified and studied. These properties include optical shocking behavior in the limit of very small dispersion, chaotic and homoclinic behavior in discretizations of the partial differential equation, and random behavior.

Nonlinear grid error effects on numerical solution of partial differential equations

NASA Technical Reports Server (NTRS)

Dey, S. K.

1980-01-01

Finite difference solutions of nonlinear partial differential equations require discretizations and consequently grid errors are generated. These errors strongly affect stability and convergence properties of difference models. Previously such errors were analyzed by linearizing the difference equations for solutions. Properties of mappings of decadence were used to analyze nonlinear instabilities. Such an analysis is directly affected by initial/boundary conditions. An algorithm was developed, applied to nonlinear Burgers equations, and verified computationally. A preliminary test shows that Navier-Stokes equations may be treated similarly.

Nonlinear Waves.

DTIC Science & Technology

1988-02-01

in Multi- dimensions II, P.M. Santini and A.S. Fokas, preprint INS#67, 1986. The Recursion Operator of the Kadomtsev - Petviashvili Equation and the...solitons, multidimensional inverse problems, Painleve equations , direct linearizations of certain nonlinear wave equations , DBAR problems, Riemann...the Navy is (a) the recent discovery that many of the equations describing ship hydrodynamics in channels of finite depth obey nonlinear equations

Cubication of Conservative Nonlinear Oscillators

ERIC Educational Resources Information Center

Belendez, Augusto; Alvarez, Mariela L.; Fernandez, Elena; Pascual, Immaculada

2009-01-01

A cubication procedure of the nonlinear differential equation for conservative nonlinear oscillators is analysed and discussed. This scheme is based on the Chebyshev series expansion of the restoring force, and this allows us to approximate the original nonlinear differential equation by a Duffing equation in which the coefficients for the linear…

FAST TRACK COMMUNICATION Quasi self-adjoint nonlinear wave equations

NASA Astrophysics Data System (ADS)

Ibragimov, N. H.; Torrisi, M.; Tracinà, R.

2010-11-01

In this paper we generalize the classification of self-adjoint second-order linear partial differential equation to a family of nonlinear wave equations with two independent variables. We find a class of quasi self-adjoint nonlinear equations which includes the self-adjoint linear equations as a particular case. The property of a differential equation to be quasi self-adjoint is important, e.g. for constructing conservation laws associated with symmetries of the differential equation.

S{sub 2}SA preconditioning for the S{sub n} equations with strictly non negative spatial discretization

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

Bruss, D. E.; Morel, J. E.; Ragusa, J. C.

2013-07-01

Preconditioners based upon sweeps and diffusion-synthetic acceleration have been constructed and applied to the zeroth and first spatial moments of the 1-D S{sub n} transport equation using a strictly non negative nonlinear spatial closure. Linear and nonlinear preconditioners have been analyzed. The effectiveness of various combinations of these preconditioners are compared. In one dimension, nonlinear sweep preconditioning is shown to be superior to linear sweep preconditioning, and DSA preconditioning using nonlinear sweeps in conjunction with a linear diffusion equation is found to be essentially equivalent to nonlinear sweeps in conjunction with a nonlinear diffusion equation. The ability to use amore » linear diffusion equation has important implications for preconditioning the S{sub n} equations with a strictly non negative spatial discretization in multiple dimensions. (authors)« less

Wiggly tails: A gravitational wave signature of massive fields around black holes

NASA Astrophysics Data System (ADS)

Degollado, Juan Carlos; Herdeiro, Carlos A. R.

2014-09-01

Massive fields can exist in long-lived configurations around black holes. We examine how the gravitational wave signal of a perturbed black hole is affected by such "dirtiness" within linear theory. As a concrete example, we consider the gravitational radiation emitted by the infall of a massive scalar field into a Schwarzschild black hole. Whereas part of the scalar field is absorbed/scattered by the black hole and triggers gravitational wave emission, another part lingers in long-lived quasibound states. Solving numerically the Teukolsky master equation for gravitational perturbations coupled to the massive Klein-Gordon equation, we find a characteristic gravitational wave signal, composed by a quasinormal ringing followed by a late time tail. In contrast to "clean" black holes, however, the late time tail contains small amplitude wiggles with the frequency of the dominating quasibound state. Additionally, an observer dependent beating pattern may also be seen. These features were already observed in fully nonlinear studies; our analysis shows they are present at linear level, and, since it reduces to a 1+1 dimensional numerical problem, allows for cleaner numerical data. Moreover, we discuss the power law of the tail and that it only becomes universal sufficiently far away from the dirty black hole. The wiggly tails, by constrast, are a generic feature that may be used as a smoking gun for the presence of massive fields around black holes, either as a linear cloud or as fully nonlinear hair.

A family of nonlinear Schrödinger equations admitting q-plane wave solutions

NASA Astrophysics Data System (ADS)

Nobre, F. D.; Plastino, A. R.

2017-08-01

Nonlinear Schrödinger equations with power-law nonlinearities have attracted considerable attention recently. Two previous proposals for these types of equations, corresponding respectively to the Gross-Pitaievsky equation and to the one associated with nonextensive statistical mechanics, are here unified into a single, parameterized family of nonlinear Schrödinger equations. Power-law nonlinear terms characterized by exponents depending on a real index q, typical of nonextensive statistical mechanics, are considered in such a way that the Gross-Pitaievsky equation is recovered in the limit q → 1. A classical field theory shows that, due to these nonlinearities, an extra field Φ (x → , t) (besides the usual one Ψ (x → , t)) must be introduced for consistency. The new field can be identified with Ψ* (x → , t) only when q → 1. For q ≠ 1 one has a pair of coupled nonlinear wave equations governing the joint evolution of the complex valued fields Ψ (x → , t) and Φ (x → , t). These equations reduce to the usual pair of complex-conjugate ones only in the q → 1 limit. Interestingly, the nonlinear equations obeyed by Ψ (x → , t) and Φ (x → , t) exhibit a common, soliton-like, traveling solution, which is expressible in terms of the q-exponential function that naturally emerges within nonextensive statistical mechanics.

Nonlinear Poisson Equation for Heterogeneous Media

PubMed Central

Hu, Langhua; Wei, Guo-Wei

2012-01-01

The Poisson equation is a widely accepted model for electrostatic analysis. However, the Poisson equation is derived based on electric polarizations in a linear, isotropic, and homogeneous dielectric medium. This article introduces a nonlinear Poisson equation to take into consideration of hyperpolarization effects due to intensive charges and possible nonlinear, anisotropic, and heterogeneous media. Variational principle is utilized to derive the nonlinear Poisson model from an electrostatic energy functional. To apply the proposed nonlinear Poisson equation for the solvation analysis, we also construct a nonpolar solvation energy functional based on the nonlinear Poisson equation by using the geometric measure theory. At a fixed temperature, the proposed nonlinear Poisson theory is extensively validated by the electrostatic analysis of the Kirkwood model and a set of 20 proteins, and the solvation analysis of a set of 17 small molecules whose experimental measurements are also available for a comparison. Moreover, the nonlinear Poisson equation is further applied to the solvation analysis of 21 compounds at different temperatures. Numerical results are compared to theoretical prediction, experimental measurements, and those obtained from other theoretical methods in the literature. A good agreement between our results and experimental data as well as theoretical results suggests that the proposed nonlinear Poisson model is a potentially useful model for electrostatic analysis involving hyperpolarization effects. PMID:22947937

Solving Nonlinear Coupled Differential Equations

NASA Technical Reports Server (NTRS)

Mitchell, L.; David, J.

1986-01-01

Harmonic balance method developed to obtain approximate steady-state solutions for nonlinear coupled ordinary differential equations. Method usable with transfer matrices commonly used to analyze shaft systems. Solution to nonlinear equation, with periodic forcing function represented as sum of series similar to Fourier series but with form of terms suggested by equation itself.

Nonlinear saturation of the slab ITG instability and zonal flow generation with fully kinetic ions

NASA Astrophysics Data System (ADS)

Miecnikowski, Matthew T.; Sturdevant, Benjamin J.; Chen, Yang; Parker, Scott E.

2018-05-01

Fully kinetic turbulence models are of interest for their potential to validate or replace gyrokinetic models in plasma regimes where the gyrokinetic expansion parameters are marginal. Here, we demonstrate fully kinetic ion capability by simulating the growth and nonlinear saturation of the ion-temperature-gradient instability in shearless slab geometry assuming adiabatic electrons and including zonal flow dynamics. The ion trajectories are integrated using the Lorentz force, and the cyclotron motion is fully resolved. Linear growth and nonlinear saturation characteristics show excellent agreement with analogous gyrokinetic simulations across a wide range of parameters. The fully kinetic simulation accurately reproduces the nonlinearly generated zonal flow. This work demonstrates nonlinear capability, resolution of weak gradient drive, and zonal flow physics, which are critical aspects of modeling plasma turbulence with full ion dynamics.

Theoretical and computational studies of the sheath of a planar wall

NASA Astrophysics Data System (ADS)

Giraudo, Martina; Camporeale, Enrico; Delzanno, Gian Luca; Lapenta, Giovanni

2012-03-01

We present an investigation of the stability and nonlinear evolution of the sheath of a planar wall. We focus on the electrostatic limit. The stability analysis is conducted with a fluid model where continuity and momentum equations for the electrons and ions are coupled through Poisson's equation. The effect of electron emission from the wall is studied parametrically. Our results show that a sheath instability associated with the emitted electrons can exist. Following Ref. [1], it is interpreted as a Rayleigh-Taylor instability driven by the favorable combination of the sheath electron density gradient and electric field. Fully kinetic Particle-In-Cell (PIC) simulations will also be presented to investigate whether this instability indeed exists and to study the nonlinear effect of electron emission on the sheath profiles. The simulations will be conducted with CPIC, a new electrostatic PIC code that couples the standard PIC algorithm with strategies for generation and adaptation of the computational grid. [4pt] [1] G.L. Delzanno, ``A paradigm for the stability of the plasma sheath against fluid perturbations,'' Phys. Plasmas 18, 103508 (2011).

Preliminary Analysis of a Submerged Wave Energy Device

NASA Astrophysics Data System (ADS)

Wagner, J. R.; Wagner, J. J.; Hayatdavoodi, M.; Ertekin, R. C.

2016-02-01

Preliminary analysis of a submerged wave energy harvesting device is presented. The device is composed of a thin, horizontally submerged plate that is restricted to heave oscillations under the influence of surface waves. The submerged plate is oscillating, and it can be attached to a fixed rotor, or a piston, to harvest the wave energy. A fully submerged wave energy converter is preferred over a surface energy convertor due to its durability and less visual and physical distractions it presents. In this study, the device is subject to nonlinear shallow-water waves. Wave loads on the submerged oscillating plate are obtained via the Level I Green-Naghdi equations. The unsteady motion of the plate is obtained by solving the nonlinear equations of motion. The results are obtained for a range of waves with varying heights and periods. The amplitude and period of plate oscillations are analyzed as functions of the wave parameters and plate width. Particular attention is given to the selection of the site of desired wave field. Initial estimation on the amount of energy extraction from the device, located near shore at a given site, is provided.

Applications of the ETEM for obtaining optical soliton solutions for the Lakshmanan-Porsezian-Daniel model

NASA Astrophysics Data System (ADS)

Manafian, Jalil; Foroutan, Mohammadreza; Guzali, Aref

2017-11-01

This paper examines the effectiveness of an integration scheme which is called the extended trial equation method (ETEM) for solving a well-known nonlinear equation of partial differential equations (PDEs). In this respect, the Lakshmanan-Porsezian-Daniel (LPD) equation with Kerr and power laws of nonlinearity which describes higher-order dispersion, full nonlinearity and spatiotemporal dispersion is considered, and as an achievement, a series of exact travelling-wave solutions for the aforementioned equation is formally extracted. Explicit new exact solutions are derived in different form such as dark solitons, bright solitons, solitary wave, periodic solitary wave, rational function, and elliptic function solutions of LPD equation. The movement of obtained solutions is shown graphically, which helps to understand the physical phenomena of this optical soliton equation. Many other such types of nonlinear equations arising in basic fabric of communications network technology and nonlinear optics can also be solved by this method.

A procedure to construct exact solutions of nonlinear fractional differential equations.

PubMed

Güner, Özkan; Cevikel, Adem C

2014-01-01

We use the fractional transformation to convert the nonlinear partial fractional differential equations with the nonlinear ordinary differential equations. The Exp-function method is extended to solve fractional partial differential equations in the sense of the modified Riemann-Liouville derivative. We apply the Exp-function method to the time fractional Sharma-Tasso-Olver equation, the space fractional Burgers equation, and the time fractional fmKdV equation. As a result, we obtain some new exact solutions.

Self-Consistent Frequency Sweeping of TAE mode

NASA Astrophysics Data System (ADS)

Wang, Ge

2012-03-01

We have extended our intuitive Toroidal Alfven Wave (TAE) model [1] for describing spontaneous frequency sweeping by a destabilizing component of energetic particles. Now a fully developed self-consistent description for frequency sweeping of an isolated TAE mode has been developed. As in [1], we use the Rosenbluth, Berk,Van Dam tip theory [2], valid for low beta, large aspect ratio, circular tokamaks, to describe the evolution of the TAE wave equation. The wave is coupled to the particle dynamics that uses the Berk, Breizman, Ye map model [3] to construct the particle/wave Lagrangian associated with a phase space dependent mode structure. Then together with the appropriate Vlasov equation for describing the particle dynamics, a set of equations determining the dynamics of the system has been formulated. Adiabatic solutions have been obtained and work is underway in simulating the exact nonlinear dynamics. A status report of our results will be given at the meeting. [4pt] [1] G. Wang and H. L. Berk, Communication in Nonlinear Science and Numerical Simulation 17, 2179 (2012) [0pt] [2] M. N. Rosenbluth,; H. L. Berk, J. Van Dam and D. M. Lingberg, Phys. Rev. Lett. 68, 596 (1992). [0pt] [3] Berk, H.L.; Breizman, B.N.; Ye, H. In: Physics of Fluids B 51993, 1506 (1993)

Time-dependent theoretical treatments of the dynamics of electrons and nuclei in molecular systems

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

Deumens, E.; Diz, A.; Longo, R.

1994-07-01

An overview is presented of methods for time-dependent treatments of molecules as systems of electrons and nuclei. The theoretical details of these methods are reviewed and contrasted in the light of a recently developed time-dependent method called electron-nuclear dynamics. Electron-nuclear dynamics (END) is a formulation of the complete dynamics of electrons and nuclei of a molecular system that eliminates the necessity of constructing potential-energy surfaces. Because of its general formulation, it encompasses many aspects found in other formulations and can serve as a didactic device for clarifying many of the principles and approximations relevant in time-dependent treatments of molecular systems.more » The END equations are derived from the time-dependent variational principle applied to a chosen family of efficiently parametrized approximate state vectors. A detailed analysis of the END equations is given for the case of a single-determinantal state for the electrons and a classical treatment of the nuclei. The approach leads to a simple formulation of the fully nonlinear time-dependent Hartree-Fock theory including nuclear dynamics. The nonlinear END equations with the [ital ab] [ital initio] Coulomb Hamiltonian have been implemented at this level of theory in a computer program, ENDyne, and have been shown feasible for the study of small molecular systems. Implementation of the Austin Model 1 semiempirical Hamiltonian is discussed as a route to large molecular systems. The linearized END equations at this level of theory are shown to lead to the random-phase approximation for the coupled system of electrons and nuclei. The qualitative features of the general nonlinear solution are analyzed using the results of the linearized equations as a first approximation. Some specific applications of END are presented, and the comparison with experiment and other theoretical approaches is discussed.« less

Nonlinear Waves In A Stenosed Elastic Tube Filled With Viscous Fluid: Forced Perturbed Korteweg-De Vries Equation

NASA Astrophysics Data System (ADS)

Gaik*, Tay Kim; Demiray, Hilmi; Tiong, Ong Chee

In the present work, treating the artery as a prestressed thin-walled and long circularly cylindrical elastic tube with a mild symmetrical stenosis and the blood as an incompressible Newtonian fluid, we have studied the pro pagation of weakly nonlinear waves in such a composite medium, in the long wave approximation, by use of the reductive perturbation method. By intro ducing a set of stretched coordinates suitable for the boundary value type of problems and expanding the field variables into asymptotic series of the small-ness parameter of nonlinearity and dispersion, we obtained a set of nonlinear differential equations governing the terms at various order. By solving these nonlinear differential equations, we obtained the forced perturbed Korteweg-de Vries equation with variable coefficient as the nonlinear evolution equation. By use of the coordinate transformation, it is shown that this type of nonlinear evolution equation admits a progressive wave solution with variable wave speed.

Dispersive optical soliton solutions for the hyperbolic and cubic-quintic nonlinear Schrödinger equations via the extended sinh-Gordon equation expansion method

NASA Astrophysics Data System (ADS)

Seadawy, Aly R.; Kumar, Dipankar; Chakrabarty, Anuz Kumar

2018-05-01

The (2+1)-dimensional hyperbolic and cubic-quintic nonlinear Schrödinger equations describe the propagation of ultra-short pulses in optical fibers of nonlinear media. By using an extended sinh-Gordon equation expansion method, some new complex hyperbolic and trigonometric functions prototype solutions for two nonlinear Schrödinger equations were derived. The acquired new complex hyperbolic and trigonometric solutions are expressed by dark, bright, combined dark-bright, singular and combined singular solitons. The obtained results are more compatible than those of other applied methods. The extended sinh-Gordon equation expansion method is a more powerful and robust mathematical tool for generating new optical solitary wave solutions for many other nonlinear evolution equations arising in the propagation of optical pulses.

Nonlinear Transient Growth and Boundary Layer Transition

NASA Technical Reports Server (NTRS)

Paredes, Pedro; Choudhari, Meelan M.; Li, Fei

2016-01-01

Parabolized stability equations (PSE) are used in a variational approach to study the optimal, non-modal disturbance growth in a Mach 3 at plate boundary layer and a Mach 6 circular cone boundary layer. As noted in previous works, the optimal initial disturbances correspond to steady counter-rotating streamwise vortices, which subsequently lead to the formation of streamwise-elongated structures, i.e., streaks, via a lift-up effect. The nonlinear evolution of the linearly optimal stationary perturbations is computed using the nonlinear plane-marching PSE for stationary perturbations. A fully implicit marching technique is used to facilitate the computation of nonlinear streaks with large amplitudes. To assess the effect of the finite-amplitude streaks on transition, the linear form of plane- marching PSE is used to investigate the instability of the boundary layer flow modified by spanwise periodic streaks. The onset of bypass transition is estimated by using an N- factor criterion based on the amplification of the streak instabilities. Results show that, for both flow configurations of interest, streaks of sufficiently large amplitude can lead to significantly earlier onset of transition than that in an unperturbed boundary layer without any streaks.

Drag reduction in channel flow using nonlinear control

NASA Technical Reports Server (NTRS)

Keefe, Laurence R.

1993-01-01

Two nonlinear control schemes have been applied to the problem of drag reduction in channel flow. Both schemes have been tested using numerical simulations at a mass flux Reynolds numbers of 4408, utilizing 2D nonlinear neutral modes for goal dynamics. The OGY-method, which requires feedback, reduces drag to 60-80 percent of the turbulent value at the same Reynolds number, and employs forcing only within a thin region near the wall. The H-method, or model-based control, fails to achieve any drag reduction when starting from a fully turbulent initial condition, but shows potential for suppressing or retarding laminar-to-turbulent transition by imposing instead a transition to a low drag, nonlinear traveling wave solution to the Navier-Stokes equation. The drag in this state corresponds to that achieved by the OGY-method. Model-based control requires no feedback, but in experiments to date has required the forcing be imposed within a thicker layer than the OGY-method. Control energy expenditures in both methods are small, representing less than 0.1 percent of the uncontrolled flow's energy.

A Leonard-Sanders-Budiansky-Koiter-Type Nonlinear Shell Theory with a Hierarchy of Transverse-Shearing Deformations

NASA Technical Reports Server (NTRS)

Nemeth, Michael P.

2013-01-01

A detailed exposition on a refined nonlinear shell theory suitable for nonlinear buckling analyses of laminated-composite shell structures is presented. This shell theory includes the classical nonlinear shell theory attributed to Leonard, Sanders, Koiter, and Budiansky as an explicit proper subset. This approach is used in order to leverage the exisiting experience base and to make the theory attractive to industry. In addition, the formalism of general tensors is avoided in order to expose the details needed to fully understand and use the theory. The shell theory is based on "small" strains and "moderate" rotations, and no shell-thinness approximations are used. As a result, the strain-displacement relations are exact within the presumptions of "small" strains and "moderate" rotations. The effects of transverse-shearing deformations are included in the theory by using analyst-defined functions to describe the through-the-thickness distributions of transverse-shearing strains. Constitutive equations for laminated-composite shells are derived without using any shell-thinness approximations, and simplified forms and special cases are presented.

Nonlinear Gyro-Landau-Fluid Equations

NASA Astrophysics Data System (ADS)

Raskolnikov, I.; Mattor, Nathan; Parker, Scott E.

1996-11-01

We present fluid equations which describe the effects of both linear and nonlinear Landau damping (wave-particle-wave effects). These are derived using a recently developed analytical method similar to renormalization group theory. (Scott E. Parker and Daniele Carati, Phys. Rev. Lett. 75), 441 (1995). In this technique, the phase space structure inherent in Landau damping is treated analytically by building a ``renormalized collisionality'' onto a bare collisionality (which may be taken as vanishingly small). Here we apply this technique to the nonlinear ion gyrokinetic equation in slab geometry, obtaining nonlinear fluid equations for density, parallel momentum and heat. Wave-particle resonances are described by two functions appearing in the heat equation: a renormalized ``collisionality'' and a renormalized nonlinear coupling coeffient. It will be shown that these new equations may correct a deficiency in existing gyrofluid equations, (G. W. Hammett and F. W. Perkins, Phys. Rev. Lett. 64,) 3019 (1990). which can severely underestimate the strength of nonlinear interaction in regimes where linear resonance is strong. (N. Mattor, Phys. Fluids B 4,) 3952 (1992).

Flow model for open-channel reach or network

USGS Publications Warehouse

Schaffranek, R.W.

1987-01-01

Formulation of a one-dimensional model for simulating unsteady flow in a single open-channel reach or in a network of interconnected channels is presented. The model is both general and flexible in that it can be used to simulate a wide range of flow conditions for various channel configurations. It is based on a four-point (box), implicit, finite-difference approximation of the governing nonlinear flow equations with user-definable weighting coefficients to permit varying the solution scheme from box-centered to fully forward. Unique transformation equations are formulated that permit correlation of the unknowns at the extremities of the channels, thereby reducing coefficient matrix and execution time requirements. Discharges and water-surface elevations computed at intermediate locations within a channel are determined following solution of the transformation equations. The matrix of transformation and boundary-condition equations is solved by Gauss elimination using maximum pivot strategy. Two diverse applications of the model are presented to illustrate its broad utility. (USGS)

Nonlinear and dissipative constitutive equations for coupled first-order acoustic field equations that are consistent with the generalized Westervelt equation

NASA Astrophysics Data System (ADS)

Verweij, Martin D.; Huijssen, Jacob

2006-05-01

In diagnostic medical ultrasound, it has become increasingly important to evaluate the nonlinear field of an acoustic beam that propagates in a weakly nonlinear, dissipative medium and that is steered off-axis up to very wide angles. In this case, computations cannot be based on the widely used KZK equation since it applies only to small angles. To benefit from successful computational schemes from elastodynamics and electromagnetics, we propose to use two first-order acoustic field equations, accompanied by two constitutive equations, as an alternative basis. This formulation quite naturally results in the contrast source formalism, makes a clear distinction between fundamental conservation laws and medium behavior, and allows for a straightforward inclusion of any medium inhomogenities. This paper is concerned with the derivation of relevant constitutive equations. We take a pragmatic approach and aim to find those constitutive equations that represent the same medium as implicitly described by the recognized, full wave, nonlinear equations such as the generalized Westervelt equation. We will show how this is achieved by considering the nonlinear case without attenuation, the linear case with attenuation, and the nonlinear case with attenuation. As a result we will obtain surprisingly simple constitutive equations for the full wave case.

A Procedure to Construct Exact Solutions of Nonlinear Fractional Differential Equations

PubMed Central

Güner, Özkan; Cevikel, Adem C.

2014-01-01

We use the fractional transformation to convert the nonlinear partial fractional differential equations with the nonlinear ordinary differential equations. The Exp-function method is extended to solve fractional partial differential equations in the sense of the modified Riemann-Liouville derivative. We apply the Exp-function method to the time fractional Sharma-Tasso-Olver equation, the space fractional Burgers equation, and the time fractional fmKdV equation. As a result, we obtain some new exact solutions. PMID:24737972

Bayesian Analysis of Structural Equation Models with Nonlinear Covariates and Latent Variables

ERIC Educational Resources Information Center

Song, Xin-Yuan; Lee, Sik-Yum

2006-01-01

In this article, we formulate a nonlinear structural equation model (SEM) that can accommodate covariates in the measurement equation and nonlinear terms of covariates and exogenous latent variables in the structural equation. The covariates can come from continuous or discrete distributions. A Bayesian approach is developed to analyze the…

Exact multisoliton solutions of general nonlinear Schrödinger equation with derivative.

PubMed

Li, Qi; Duan, Qiu-yuan; Zhang, Jian-bing

2014-01-01

Multisoliton solutions are derived for a general nonlinear Schrödinger equation with derivative by using Hirota's approach. The dynamics of one-soliton solution and two-soliton interactions are also illustrated. The considered equation can reduce to nonlinear Schrödinger equation with derivative as well as the solutions.

Convergence of Galerkin approximations for operator Riccati equations: A nonlinear evolution equation approach

NASA Technical Reports Server (NTRS)

Rosen, I. G.

1988-01-01

An approximation and convergence theory was developed for Galerkin approximations to infinite dimensional operator Riccati differential equations formulated in the space of Hilbert-Schmidt operators on a separable Hilbert space. The Riccati equation was treated as a nonlinear evolution equation with dynamics described by a nonlinear monotone perturbation of a strongly coercive linear operator. A generic approximation result was proven for quasi-autonomous nonlinear evolution system involving accretive operators which was then used to demonstrate the Hilbert-Schmidt norm convergence of Galerkin approximations to the solution of the Riccati equation. The application of the results was illustrated in the context of a linear quadratic optimal control problem for a one dimensional heat equation.

Current interactions from the one-form sector of nonlinear higher-spin equations

NASA Astrophysics Data System (ADS)

Gelfond, O. A.; Vasiliev, M. A.

2018-06-01

The form of higher-spin current interactions in the sector of one-forms is derived from the nonlinear higher-spin equations in AdS4. Quadratic corrections to higher-spin equations are shown to be independent of the phase of the parameter η = exp ⁡ iφ in the full nonlinear higher-spin equations. The current deformation resulting from the nonlinear higher-spin equations is represented in the canonical form with the minimal number of space-time derivatives. The non-zero spin-dependent coupling constants of the resulting currents are determined in terms of the higher-spin coupling constant η η bar . Our results confirm the conjecture that (anti-)self-dual nonlinear higher-spin equations result from the full system at (η = 0) η bar = 0.

Nonlinear Poisson equation for heterogeneous media.

PubMed

Hu, Langhua; Wei, Guo-Wei

2012-08-22

The Poisson equation is a widely accepted model for electrostatic analysis. However, the Poisson equation is derived based on electric polarizations in a linear, isotropic, and homogeneous dielectric medium. This article introduces a nonlinear Poisson equation to take into consideration of hyperpolarization effects due to intensive charges and possible nonlinear, anisotropic, and heterogeneous media. Variational principle is utilized to derive the nonlinear Poisson model from an electrostatic energy functional. To apply the proposed nonlinear Poisson equation for the solvation analysis, we also construct a nonpolar solvation energy functional based on the nonlinear Poisson equation by using the geometric measure theory. At a fixed temperature, the proposed nonlinear Poisson theory is extensively validated by the electrostatic analysis of the Kirkwood model and a set of 20 proteins, and the solvation analysis of a set of 17 small molecules whose experimental measurements are also available for a comparison. Moreover, the nonlinear Poisson equation is further applied to the solvation analysis of 21 compounds at different temperatures. Numerical results are compared to theoretical prediction, experimental measurements, and those obtained from other theoretical methods in the literature. A good agreement between our results and experimental data as well as theoretical results suggests that the proposed nonlinear Poisson model is a potentially useful model for electrostatic analysis involving hyperpolarization effects. Copyright © 2012 Biophysical Society. Published by Elsevier Inc. All rights reserved.

Relations between nonlinear Riccati equations and other equations in fundamental physics

NASA Astrophysics Data System (ADS)

Schuch, Dieter

2014-10-01

Many phenomena in the observable macroscopic world obey nonlinear evolution equations while the microscopic world is governed by quantum mechanics, a fundamental theory that is supposedly linear. In order to combine these two worlds in a common formalism, at least one of them must sacrifice one of its dogmas. Linearizing nonlinear dynamics would destroy the fundamental property of this theory, however, it can be shown that quantum mechanics can be reformulated in terms of nonlinear Riccati equations. In a first step, it will be shown that the information about the dynamics of quantum systems with analytical solutions can not only be obtainable from the time-dependent Schrödinger equation but equally-well from a complex Riccati equation. Comparison with supersymmetric quantum mechanics shows that even additional information can be obtained from the nonlinear formulation. Furthermore, the time-independent Schrödinger equation can also be rewritten as a complex Riccati equation for any potential. Extension of the Riccati formulation to include irreversible dissipative effects is straightforward. Via (real and complex) Riccati equations, other fields of physics can also be treated within the same formalism, e.g., statistical thermodynamics, nonlinear dynamical systems like those obeying a logistic equation as well as wave equations in classical optics, Bose- Einstein condensates and cosmological models. Finally, the link to abstract "quantizations" such as the Pythagorean triples and Riccati equations connected with trigonometric and hyperbolic functions will be shown.

Modified Taylor series method for solving nonlinear differential equations with mixed boundary conditions defined on finite intervals.

PubMed

Vazquez-Leal, Hector; Benhammouda, Brahim; Filobello-Nino, Uriel Antonio; Sarmiento-Reyes, Arturo; Jimenez-Fernandez, Victor Manuel; Marin-Hernandez, Antonio; Herrera-May, Agustin Leobardo; Diaz-Sanchez, Alejandro; Huerta-Chua, Jesus

2014-01-01

In this article, we propose the application of a modified Taylor series method (MTSM) for the approximation of nonlinear problems described on finite intervals. The issue of Taylor series method with mixed boundary conditions is circumvented using shooting constants and extra derivatives of the problem. In order to show the benefits of this proposal, three different kinds of problems are solved: three-point boundary valued problem (BVP) of third-order with a hyperbolic sine nonlinearity, two-point BVP for a second-order nonlinear differential equation with an exponential nonlinearity, and a two-point BVP for a third-order nonlinear differential equation with a radical nonlinearity. The result shows that the MTSM method is capable to generate easily computable and highly accurate approximations for nonlinear equations. 34L30.

An analytically soluble problem in fully nonlinear statistical gravitational lensing

NASA Technical Reports Server (NTRS)

Schneider, P.

1987-01-01

The amplification probability distribution p(I)dI for a point source behind a random star field which acts as the deflector exhibits a I exp-3 behavior for large amplification, as can be shown from the universality of the lens equation near critical lines. In this paper it is shown that the amplitude of the I exp-3 tail can be derived exactly for arbitrary mass distribution of the stars, surface mass density of stars and smoothly distributed matter, and large-scale shear. This is then compared with the corresponding linear result.

Solutions of the cylindrical nonlinear Maxwell equations.

PubMed

Xiong, Hao; Si, Liu-Gang; Ding, Chunling; Lü, Xin-You; Yang, Xiaoxue; Wu, Ying

2012-01-01

Cylindrical nonlinear optics is a burgeoning research area which describes cylindrical electromagnetic wave propagation in nonlinear media. Finding new exact solutions for different types of nonlinearity and inhomogeneity to describe cylindrical electromagnetic wave propagation is of great interest and meaningful for theory and application. This paper gives exact solutions for the cylindrical nonlinear Maxwell equations and presents an interesting connection between the exact solutions for different cylindrical nonlinear Maxwell equations. We also provide some examples and discussion to show the application of the results we obtained. Our results provide the basis for solving complex systems of nonlinearity and inhomogeneity with simple systems.

Fully vectorial laser resonator modeling of continuous-wave solid-state lasers including rate equations, thermal lensing and stress-induced birefringence.

PubMed

Asoubar, Daniel; Wyrowski, Frank

2015-07-27

The computer-aided design of high quality mono-mode, continuous-wave solid-state lasers requires fast, flexible and accurate simulation algorithms. Therefore in this work a model for the calculation of the transversal dominant mode structure is introduced. It is based on the generalization of the scalar Fox and Li algorithm to a fully-vectorial light representation. To provide a flexible modeling concept of different resonator geometries containing various optical elements, rigorous and approximative solutions of Maxwell's equations are combined in different subdomains of the resonator. This approach allows the simulation of plenty of different passive intracavity components as well as active media. For the numerically efficient simulation of nonlinear gain, thermal lensing and stress-induced birefringence effects in solid-state active crystals a semi-analytical vectorial beam propagation method is discussed in detail. As a numerical example the beam quality and output power of a flash-lamp-pumped Nd:YAG laser are improved. To that end we compensate the influence of stress-induced birefringence and thermal lensing by an aspherical mirror and a 90° quartz polarization rotator.

Chirped solitary pulses for a nonic nonlinear Schrödinger equation on a continuous-wave background

NASA Astrophysics Data System (ADS)

Triki, Houria; Porsezian, K.; Choudhuri, Amitava; Dinda, P. Tchofo

2016-06-01

A class of derivative nonlinear Schrödinger equation with cubic-quintic-septic-nonic nonlinear terms describing the propagation of ultrashort optical pulses through a nonlinear medium with higher-order Kerr responses is investigated. An intensity-dependent chirp ansatz is adopted for solving the two coupled amplitude-phase nonlinear equations of the propagating wave. We find that the dynamics of field amplitude in this system is governed by a first-order nonlinear ordinary differential equation with a tenth-degree nonlinear term. We demonstrate that this system allows the propagation of a very rich variety of solitary waves (kink, dark, bright, and gray solitary pulses) which do not coexist in the conventional nonlinear systems that have appeared so far in the literature. The stability of the solitary wave solution under some violation on the parametric conditions is investigated. Moreover, we show that, unlike conventional systems, the nonlinear Schrödinger equation considered here meets the special requirements for the propagation of a chirped solitary wave on a continuous-wave background, involving a balance among group velocity dispersion, self-steepening, and higher-order nonlinearities of different nature.

Bright and singular soliton solutions of the conformable time-fractional Klein-Gordon equations with different nonlinearities

NASA Astrophysics Data System (ADS)

Hosseini, Kamyar; Mayeli, Peyman; Ansari, Reza

2018-07-01

Finding the exact solutions of nonlinear fractional differential equations has gained considerable attention, during the past two decades. In this paper, the conformable time-fractional Klein-Gordon equations with quadratic and cubic nonlinearities are studied. Several exact soliton solutions, including the bright (non-topological) and singular soliton solutions are formally extracted by making use of the ansatz method. Results demonstrate that the method can efficiently handle the time-fractional Klein-Gordon equations with different nonlinearities.

On solutions of the fifth-order dispersive equations with porous medium type non-linearity

NASA Astrophysics Data System (ADS)

Kocak, Huseyin; Pinar, Zehra

2018-07-01

In this work, we focus on obtaining the exact solutions of the fifth-order semi-linear and non-linear dispersive partial differential equations, which have the second-order diffusion-like (porous-type) non-linearity. The proposed equations were not studied in the literature in the sense of the exact solutions. We reveal solutions of the proposed equations using the classical Riccati equations method. The obtained exact solutions, which can play a key role to simulate non-linear waves in the medium with dispersion and diffusion, are illustrated and discussed in details.

A more fundamental approach to the derivation of nonlinear acoustic wave equations with fractional loss operators (L).

PubMed

Prieur, Fabrice; Vilenskiy, Gregory; Holm, Sverre

2012-10-01

A corrected derivation of nonlinear wave propagation equations with fractional loss operators is presented. The fundamental approach is based on fractional formulations of the stress-strain and heat flux definitions but uses the energy equation and thermodynamic identities to link density and pressure instead of an erroneous fractional form of the entropy equation as done in Prieur and Holm ["Nonlinear acoustic wave equations with fractional loss operators," J. Acoust. Soc. Am. 130(3), 1125-1132 (2011)]. The loss operator of the obtained nonlinear wave equations differs from the previous derivations as well as the dispersion equation, but when approximating for low frequencies the expressions for the frequency dependent attenuation and velocity dispersion remain unchanged.

Fluid equations with nonlinear wave-particle resonances^

NASA Astrophysics Data System (ADS)

Mattor, Nathan

1997-11-01

We have derived fluid equations that include linear and nonlinear wave-particle resonance effects. This greatly extends previous ``Landau-fluid'' closures, which include linear Landau damping. (G.W. Hammett and F.W. Perkins, Phys. Rev. Lett. 64,) 3019 (1990).^, (Z. Chang and J. D. Callen, Phys. Fluids B 4,) 1167 (1992). The new fluid equations are derived with no approximation regarding nonlinear kinetic interaction, and so additionally include numerous nonlinear kinetic effects. The derivation starts with the electrostatic drift kinetic equation for simplicity, with a Maxwellian distribution function. Fluid closure is accomplished through a simple integration trick applied to the drift kinetic equation, using the property that the nth moment of Maxwellian distribution is related to the nth derivative. The result is a compact closure term appearing in the highest moment equation, a term which involves a plasma dispersion function of the electrostatic field and its derivatives. The new term reduces to the linear closures in appropriate limits, so both approaches retain linear Landau damping. But the nonlinearly closed equations have additional desirable properties. Unlike linear closures, the nonlinear closure retains the time-reversibility of the original kinetic equation. We have shown directly that the nonlinear closure retains at least two nonlinear resonance effects: wave-particle trapping and Compton scattering. Other nonlinear kinetic effects are currently under investigation. The new equations correct two previous discrepancies between kinetic and Landau-fluid predictions, including a propagator discrepancy (N. Mattor, Phys. Fluids B 4,) 3952 (1992). and a numerical discrepancy for the 3-mode shearless bounded slab ITG problem. (S. E. Parker et al.), Phys. Plasmas 1, 1461 (1994). ^* In collaboration with S. E. Parker, Department of Physics, University of Colorado, Boulder. ^ Work performed at LLNL under DoE contract No. W7405-ENG-48.

Overview of Edge Simulation Laboratory (ESL)

NASA Astrophysics Data System (ADS)

Cohen, R. H.; Dorr, M.; Hittinger, J.; Rognlien, T.; Umansky, M.; Xiong, A.; Xu, X.; Belli, E.; Candy, J.; Snyder, P.; Colella, P.; Martin, D.; Sternberg, T.; van Straalen, B.; Bodi, K.; Krasheninnikov, S.

2006-10-01

The ESL is a new collaboration to build a full-f electromagnetic gyrokinetic code for tokamak edge plasmas using continuum methods. Target applications are edge turbulence and transport (neoclassical and anomalous), and edge-localized modes. Initially the project has three major threads: (i) verification and validation of TEMPEST, the project's initial (electrostatic) edge code which can be run in 4D (neoclassical and transport-timescale applications) or 5D (turbulence); (ii) design of the next generation code, which will include more complete physics (electromagnetics, fluid equation option, improved collisions) and advanced numerics (fully conservative, high-order discretization, mapped multiblock grids, adaptivity), and (iii) rapid-prototype codes to explore the issues attached to solving fully nonlinear gyrokinetics with steep radial gradiens. We present a brief summary of the status of each of these activities.

From nonlinear Schrödinger hierarchy to some (2+1)-dimensional nonlinear pseudodifferential equations

NASA Astrophysics Data System (ADS)

Yang, Xiao; Du, Dianlou

2010-08-01

The Poisson structure on CN×RN is introduced to give the Hamiltonian system associated with a spectral problem which yields the nonlinear Schrödinger (NLS) hierarchy. The Hamiltonian system is proven to be Liouville integrable. Some (2+1)-dimensional equations including NLS equation, Kadomtesev-Petviashvili I (KPI) equation, coupled KPI equation, and modified Kadomtesev-Petviashvili (mKP) equation, are decomposed into Hamilton flows via the NLS hierarchy. The algebraic curve, Abel-Jacobi coordinates, and Riemann-Jacobi inversion are used to obtain the algebrogeometric solutions of these equations.

Lipschitz regularity results for nonlinear strictly elliptic equations and applications

NASA Astrophysics Data System (ADS)

Ley, Olivier; Nguyen, Vinh Duc

2017-10-01

Most of Lipschitz regularity results for nonlinear strictly elliptic equations are obtained for a suitable growth power of the nonlinearity with respect to the gradient variable (subquadratic for instance). For equations with superquadratic growth power in gradient, one usually uses weak Bernstein-type arguments which require regularity and/or convex-type assumptions on the gradient nonlinearity. In this article, we obtain new Lipschitz regularity results for a large class of nonlinear strictly elliptic equations with possibly arbitrary growth power of the Hamiltonian with respect to the gradient variable using some ideas coming from Ishii-Lions' method. We use these bounds to solve an ergodic problem and to study the regularity and the large time behavior of the solution of the evolution equation.

The fluid-dynamic paradigm of the dust-acoustic soliton

NASA Astrophysics Data System (ADS)

McKenzie, J. F.

2002-06-01

In most studies, the properties of dust-acoustic solitons are derived from the first integral of the Poisson equation, in which the shape of the pseudopotential determines both the conditions in which a soliton may exist and its amplitude. Here this first integral is interpreted as conservation of total momentum, which, along with the Bernoulli-like energy equations for each species, may be cast as the structure equation for the dust (or heavy-ion) speed in the wave. In this fluid-dynamic picture, the significance of the sonic points of each species becomes apparent. In the wave, the heavy-ion (or dust) flow speed is supersonic (relative to its sound speed), whereas the protons and electrons are subsonic (relative to their sound speeds), and the dust flow is driven towards its sonic point. It is this last feature that limits the strength (amplitude) of the wave, since the equilibrium point (the centre of the wave) must be reached before the dust speed becomes sonic. The wave is characterized by a compression in the heavies and a compression (rarefaction) in the electrons and a rarefaction (compression) in the protons if the heavies have positive (negative) charge, and the corresponding potential is a hump (dip). These features are elucidated by an exact analytical soliton, in a special case, which provides the fully nonlinear counterpoint to the weakly nonlinear sech2-type solitons associated with the Korteweg de Vries equation, and indicates the parameter regimes in which solitons may exist.

Explicit expressions for meromorphic solutions of autonomous nonlinear ordinary differential equations

NASA Astrophysics Data System (ADS)

Demina, Maria V.; Kudryashov, Nikolay A.

2011-03-01

Meromorphic solutions of autonomous nonlinear ordinary differential equations are studied. An algorithm for constructing meromorphic solutions in explicit form is presented. General expressions for meromorphic solutions (including rational, periodic, elliptic) are found for a wide class of autonomous nonlinear ordinary differential equations.

Exact traveling wave solutions of modified KdV-Zakharov-Kuznetsov equation and viscous Burgers equation.

PubMed

Islam, Md Hamidul; Khan, Kamruzzaman; Akbar, M Ali; Salam, Md Abdus

2014-01-01

Mathematical modeling of many physical systems leads to nonlinear evolution equations because most physical systems are inherently nonlinear in nature. The investigation of traveling wave solutions of nonlinear partial differential equations (NPDEs) plays a significant role in the study of nonlinear physical phenomena. In this article, we construct the traveling wave solutions of modified KDV-ZK equation and viscous Burgers equation by using an enhanced (G '/G) -expansion method. A number of traveling wave solutions in terms of unknown parameters are obtained. Derived traveling wave solutions exhibit solitary waves when special values are given to its unknown parameters. 35C07; 35C08; 35P99.

Oscillating solutions for nonlinear Helmholtz equations

NASA Astrophysics Data System (ADS)

Mandel, Rainer; Montefusco, Eugenio; Pellacci, Benedetta

2017-12-01

Existence results for radially symmetric oscillating solutions for a class of nonlinear autonomous Helmholtz equations are given and their exact asymptotic behaviour at infinity is established. Some generalizations to nonautonomous radial equations as well as existence results for nonradial solutions are found. Our theorems prove the existence of standing waves solutions of nonlinear Klein-Gordon or Schrödinger equations with large frequencies.

Hazardous Continuation Backward in Time in Nonlinear Parabolic Equations, and an Experiment in Deblurring Nonlinearly Blurred Imagery

PubMed Central

Carasso, Alfred S

2013-01-01

Identifying sources of ground water pollution, and deblurring nanoscale imagery as well as astronomical galaxy images, are two important applications involving numerical computation of parabolic equations backward in time. Surprisingly, very little is known about backward continuation in nonlinear parabolic equations. In this paper, an iterative procedure originating in spectroscopy in the 1930’s, is adapted into a useful tool for solving a wide class of 2D nonlinear backward parabolic equations. In addition, previously unsuspected difficulties are uncovered that may preclude useful backward continuation in parabolic equations deviating too strongly from the linear, autonomous, self adjoint, canonical model. This paper explores backward continuation in selected 2D nonlinear equations, by creating fictitious blurred images obtained by using several sharp images as initial data in these equations, and capturing the corresponding solutions at some positive time T. Successful backward continuation from t=T to t = 0, would recover the original sharp image. Visual recognition provides meaningful evaluation of the degree of success or failure in the reconstructed solutions. Instructive examples are developed, illustrating the unexpected influence of certain types of nonlinearities. Visually and statistically indistinguishable blurred images are presented, with vastly different deblurring results. These examples indicate that how an image is nonlinearly blurred is critical, in addition to the amount of blur. The equations studied represent nonlinear generalizations of Brownian motion, and the blurred images may be interpreted as visually expressing the results of novel stochastic processes. PMID:26401430

Hazardous Continuation Backward in Time in Nonlinear Parabolic Equations, and an Experiment in Deblurring Nonlinearly Blurred Imagery.

PubMed

Carasso, Alfred S

2013-01-01

Identifying sources of ground water pollution, and deblurring nanoscale imagery as well as astronomical galaxy images, are two important applications involving numerical computation of parabolic equations backward in time. Surprisingly, very little is known about backward continuation in nonlinear parabolic equations. In this paper, an iterative procedure originating in spectroscopy in the 1930's, is adapted into a useful tool for solving a wide class of 2D nonlinear backward parabolic equations. In addition, previously unsuspected difficulties are uncovered that may preclude useful backward continuation in parabolic equations deviating too strongly from the linear, autonomous, self adjoint, canonical model. This paper explores backward continuation in selected 2D nonlinear equations, by creating fictitious blurred images obtained by using several sharp images as initial data in these equations, and capturing the corresponding solutions at some positive time T. Successful backward continuation from t=T to t = 0, would recover the original sharp image. Visual recognition provides meaningful evaluation of the degree of success or failure in the reconstructed solutions. Instructive examples are developed, illustrating the unexpected influence of certain types of nonlinearities. Visually and statistically indistinguishable blurred images are presented, with vastly different deblurring results. These examples indicate that how an image is nonlinearly blurred is critical, in addition to the amount of blur. The equations studied represent nonlinear generalizations of Brownian motion, and the blurred images may be interpreted as visually expressing the results of novel stochastic processes.

Reactive Transport Modeling of Induced Calcite Precipitation Reaction Fronts in Porous Media Using A Parallel, Fully Coupled, Fully Implicit Approach

NASA Astrophysics Data System (ADS)

Guo, L.; Huang, H.; Gaston, D.; Redden, G. D.; Fox, D. T.; Fujita, Y.

2010-12-01

Inducing mineral precipitation in the subsurface is one potential strategy for immobilizing trace metal and radionuclide contaminants. Generating mineral precipitates in situ can be achieved by manipulating chemical conditions, typically through injection or in situ generation of reactants. How these reactants transport, mix and react within the medium controls the spatial distribution and composition of the resulting mineral phases. Multiple processes, including fluid flow, dispersive/diffusive transport of reactants, biogeochemical reactions and changes in porosity-permeability, are tightly coupled over a number of scales. Numerical modeling can be used to investigate the nonlinear coupling effects of these processes which are quite challenging to explore experimentally. Many subsurface reactive transport simulators employ a de-coupled or operator-splitting approach where transport equations and batch chemistry reactions are solved sequentially. However, such an approach has limited applicability for biogeochemical systems with fast kinetics and strong coupling between chemical reactions and medium properties. A massively parallel, fully coupled, fully implicit Reactive Transport simulator (referred to as “RAT”) based on a parallel multi-physics object-oriented simulation framework (MOOSE) has been developed at the Idaho National Laboratory. Within this simulator, systems of transport and reaction equations can be solved simultaneously in a fully coupled, fully implicit manner using the Jacobian Free Newton-Krylov (JFNK) method with additional advanced computing capabilities such as (1) physics-based preconditioning for solution convergence acceleration, (2) massively parallel computing and scalability, and (3) adaptive mesh refinements for 2D and 3D structured and unstructured mesh. The simulator was first tested against analytical solutions, then applied to simulating induced calcium carbonate mineral precipitation in 1D columns and 2D flow cells as analogs to homogeneous and heterogeneous porous media, respectively. In 1D columns, calcium carbonate mineral precipitation was driven by urea hydrolysis catalyzed by urease enzyme, and in 2D flow cells, calcium carbonate mineral forming reactants were injected sequentially, forming migrating reaction fronts that are typically highly nonuniform. The RAT simulation results for the spatial and temporal distributions of precipitates, reaction rates and major species in the system, and also for changes in porosity and permeability, were compared to both laboratory experimental data and computational results obtained using other reactive transport simulators. The comparisons demonstrate the ability of RAT to simulate complex nonlinear systems and the advantages of fully coupled approaches, over de-coupled methods, for accurate simulation of complex, dynamic processes such as engineered mineral precipitation in subsurface environments.

Nonlinear Problems in Fluid Dynamics and Inverse Scattering

DTIC Science & Technology

1993-05-31

nonlinear Kadomtsev - Petviashvili (KP) equations , have solutions which will become infinite in finite time. This phenomenon is sometimes referred to as...40 (November 1992). 4 7. Wave Collapse and Instability of Solitary Waves of a Generalized Nonlinear Kaoiomtsev- Petviashvili Equation , X.P. Wang, M.J...words) The inverse scattering of a class of differential-difference equations and multidimensional operators has been constructed. Solutions of nonlinear

Fractional analysis for nonlinear electrical transmission line and nonlinear Schroedinger equations with incomplete sub-equation

NASA Astrophysics Data System (ADS)

Fendzi-Donfack, Emmanuel; Nguenang, Jean Pierre; Nana, Laurent

2018-02-01

We use the fractional complex transform with the modified Riemann-Liouville derivative operator to establish the exact and generalized solutions of two fractional partial differential equations. We determine the solutions of fractional nonlinear electrical transmission lines (NETL) and the perturbed nonlinear Schroedinger (NLS) equation with the Kerr law nonlinearity term. The solutions are obtained for the parameters in the range (0<α≤1) of the derivative operator and we found the traditional solutions for the limiting case of α =1. We show that according to the modified Riemann-Liouville derivative, the solutions found can describe physical systems with memory effect, transient effects in electrical systems and nonlinear transmission lines, and other systems such as optical fiber.

A Bivariate Chebyshev Spectral Collocation Quasilinearization Method for Nonlinear Evolution Parabolic Equations

PubMed Central

Motsa, S. S.; Magagula, V. M.; Sibanda, P.

2014-01-01

This paper presents a new method for solving higher order nonlinear evolution partial differential equations (NPDEs). The method combines quasilinearisation, the Chebyshev spectral collocation method, and bivariate Lagrange interpolation. In this paper, we use the method to solve several nonlinear evolution equations, such as the modified KdV-Burgers equation, highly nonlinear modified KdV equation, Fisher's equation, Burgers-Fisher equation, Burgers-Huxley equation, and the Fitzhugh-Nagumo equation. The results are compared with known exact analytical solutions from literature to confirm accuracy, convergence, and effectiveness of the method. There is congruence between the numerical results and the exact solutions to a high order of accuracy. Tables were generated to present the order of accuracy of the method; convergence graphs to verify convergence of the method and error graphs are presented to show the excellent agreement between the results from this study and the known results from literature. PMID:25254252

A bivariate Chebyshev spectral collocation quasilinearization method for nonlinear evolution parabolic equations.

PubMed

Motsa, S S; Magagula, V M; Sibanda, P

2014-01-01

This paper presents a new method for solving higher order nonlinear evolution partial differential equations (NPDEs). The method combines quasilinearisation, the Chebyshev spectral collocation method, and bivariate Lagrange interpolation. In this paper, we use the method to solve several nonlinear evolution equations, such as the modified KdV-Burgers equation, highly nonlinear modified KdV equation, Fisher's equation, Burgers-Fisher equation, Burgers-Huxley equation, and the Fitzhugh-Nagumo equation. The results are compared with known exact analytical solutions from literature to confirm accuracy, convergence, and effectiveness of the method. There is congruence between the numerical results and the exact solutions to a high order of accuracy. Tables were generated to present the order of accuracy of the method; convergence graphs to verify convergence of the method and error graphs are presented to show the excellent agreement between the results from this study and the known results from literature.

Solitary wave for a nonintegrable discrete nonlinear Schrödinger equation in nonlinear optical waveguide arrays

NASA Astrophysics Data System (ADS)

Ma, Li-Yuan; Ji, Jia-Liang; Xu, Zong-Wei; Zhu, Zuo-Nong

2018-03-01

We study a nonintegrable discrete nonlinear Schrödinger (dNLS) equation with the term of nonlinear nearest-neighbor interaction occurred in nonlinear optical waveguide arrays. By using discrete Fourier transformation, we obtain numerical approximations of stationary and travelling solitary wave solutions of the nonintegrable dNLS equation. The analysis of stability of stationary solitary waves is performed. It is shown that the nonlinear nearest-neighbor interaction term has great influence on the form of solitary wave. The shape of solitary wave is important in the electric field propagating. If we neglect the nonlinear nearest-neighbor interaction term, much important information in the electric field propagating may be missed. Our numerical simulation also demonstrates the difference of chaos phenomenon between the nonintegrable dNLS equation with nonlinear nearest-neighbor interaction and another nonintegrable dNLS equation without the term. Project supported by the National Natural Science Foundation of China (Grant Nos. 11671255 and 11701510), the Ministry of Economy and Competitiveness of Spain (Grant No. MTM2016-80276-P (AEI/FEDER, EU)), and the China Postdoctoral Science Foundation (Grant No. 2017M621964).

Thermodynamics of viscoelastic rate-type fluids with stress diffusion

NASA Astrophysics Data System (ADS)

Málek, Josef; Průša, Vít; Skřivan, Tomáš; Süli, Endre

2018-02-01

We propose thermodynamically consistent models for viscoelastic fluids with a stress diffusion term. In particular, we derive variants of compressible/incompressible Maxwell/Oldroyd-B models with a stress diffusion term in the evolution equation for the extra stress tensor. It is shown that the stress diffusion term can be interpreted either as a consequence of a nonlocal energy storage mechanism or as a consequence of a nonlocal entropy production mechanism, while different interpretations of the stress diffusion mechanism lead to different evolution equations for the temperature. The benefits of the knowledge of the thermodynamical background of the derived models are documented in the study of nonlinear stability of equilibrium rest states. The derived models open up the possibility to study fully coupled thermomechanical problems involving viscoelastic rate-type fluids with stress diffusion.

Entanglement dynamics following a sudden quench: An exact solution

NASA Astrophysics Data System (ADS)

Ghosh, Supriyo; Gupta, Kumar S.; Srivastava, Shashi C. L.

2017-12-01

We present an exact and fully analytical treatment of the entanglement dynamics for an isolated system of N coupled oscillators following a sudden quench of the system parameters. The system is analyzed using the solutions of the time-dependent Schrodinger's equation, which are obtained by solving the corresponding nonlinear Ermakov equations. The entanglement entropies exhibit a multi-oscillatory behaviour, where the number of dynamically generated time scales increases with N. The harmonic chains exhibit entanglement revival and for larger values of N (> 10), we find near-critical logarithmic scaling for the entanglement entropy, which is modulated by a time-dependent factor. The N = 2 case is equivalent to the two-site Bose-Hubbard model in the tunneling regime, which is amenable to empirical realization in cold-atom systems.

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.

Optical analogue of relativistic Dirac solitons in binary waveguide arrays

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

Tran, Truong X., E-mail: truong.tran@mpl.mpg.de; Max Planck Institute for the Science of Light, Günther-Scharowsky str. 1, 91058 Erlangen; Longhi, Stefano

2014-01-15

We study analytically and numerically an optical analogue of Dirac solitons in binary waveguide arrays in the presence of Kerr nonlinearity. Pseudo-relativistic soliton solutions of the coupled-mode equations describing dynamics in the array are analytically derived. We demonstrate that with the found soliton solutions, the coupled mode equations can be converted into the nonlinear relativistic 1D Dirac equation. This paves the way for using binary waveguide arrays as a classical simulator of quantum nonlinear effects arising from the Dirac equation, something that is thought to be impossible to achieve in conventional (i.e. linear) quantum field theory. -- Highlights: •An opticalmore » analogue of Dirac solitons in nonlinear binary waveguide arrays is suggested. •Analytical solutions to pseudo-relativistic solitons are presented. •A correspondence of optical coupled-mode equations with the nonlinear relativistic Dirac equation is established.« less

Dispersive optical soliton solutions for higher order nonlinear Sasa-Satsuma equation in mono mode fibers via new auxiliary equation method

NASA Astrophysics Data System (ADS)

Khater, Mostafa M. A.; Seadawy, Aly R.; Lu, Dianchen

2018-01-01

In this research, we apply new technique for higher order nonlinear Schrödinger equation which is representing the propagation of short light pulses in the monomode optical fibers and the evolution of slowly varying packets of quasi-monochromatic waves in weakly nonlinear media that have dispersion. Nonlinear Schrödinger equation is one of the basic model in fiber optics. We apply new auxiliary equation method for nonlinear Sasa-Satsuma equation to obtain a new optical forms of solitary traveling wave solutions. Exact and solitary traveling wave solutions are obtained in different kinds like trigonometric, hyperbolic, exponential, rational functions, …, etc. These forms of solutions that we represent in this research prove the superiority of our new technique on almost thirteen powerful methods. The main merits of this method over the other methods are that it gives more general solutions with some free parameters.

Modified harmonic balance method for the solution of nonlinear jerk equations

NASA Astrophysics Data System (ADS)

Rahman, M. Saifur; Hasan, A. S. M. Z.

2018-03-01

In this paper, a second approximate solution of nonlinear jerk equations (third order differential equation) can be obtained by using modified harmonic balance method. The method is simpler and easier to carry out the solution of nonlinear differential equations due to less number of nonlinear equations are required to solve than the classical harmonic balance method. The results obtained from this method are compared with those obtained from the other existing analytical methods that are available in the literature and the numerical method. The solution shows a good agreement with the numerical solution as well as the analytical methods of the available literature.

N-body simulations for f(R) gravity using a self-adaptive particle-mesh code

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

Zhao Gongbo; Koyama, Kazuya; Li Baojiu

2011-02-15

We perform high-resolution N-body simulations for f(R) gravity based on a self-adaptive particle-mesh code MLAPM. The chameleon mechanism that recovers general relativity on small scales is fully taken into account by self-consistently solving the nonlinear equation for the scalar field. We independently confirm the previous simulation results, including the matter power spectrum, halo mass function, and density profiles, obtained by Oyaizu et al.[Phys. Rev. D 78, 123524 (2008)] and Schmidt et al.[Phys. Rev. D 79, 083518 (2009)], and extend the resolution up to k{approx}20 h/Mpc for the measurement of the matter power spectrum. Based on our simulation results, we discussmore » how the chameleon mechanism affects the clustering of dark matter and halos on full nonlinear scales.« less

A note on the generation of phase plane plots on a digital computer. [for solution of nonlinear differential equations

NASA Technical Reports Server (NTRS)

Simon, M. K.

1980-01-01

A technique is presented for generating phase plane plots on a digital computer which circumvents the difficulties associated with more traditional methods of numerical solving nonlinear differential equations. In particular, the nonlinear differential equation of operation is formulated.

Capillary waves in the subcritical nonlinear Schroedinger equation

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

Kozyreff, G.

2010-01-15

We expand recent results on the nonlinear Schroedinger equation with cubic-quintic nonlinearity to show that some solutions are described by the Bernoulli equation in the presence of surface tension. As a consequence, capillary waves are predicted and found numerically at the interface between regions of large and low amplitude.

Finding all solutions of nonlinear equations using the dual simplex method

NASA Astrophysics Data System (ADS)

Yamamura, Kiyotaka; Fujioka, Tsuyoshi

2003-03-01

Recently, an efficient algorithm has been proposed for finding all solutions of systems of nonlinear equations using linear programming. This algorithm is based on a simple test (termed the LP test) for nonexistence of a solution to a system of nonlinear equations using the dual simplex method. In this letter, an improved version of the LP test algorithm is proposed. By numerical examples, it is shown that the proposed algorithm could find all solutions of a system of 300 nonlinear equations in practical computation time.

Universal and integrable nonlinear evolution systems of equations in 2+1 dimensions

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

Maccari, A.

1997-08-01

Integrable systems of nonlinear partial differential equations (PDEs) are obtained from integrable equations in 2+1 dimensions, by means of a reduction method of broad applicability based on Fourier expansion and spatio{endash}temporal rescalings, which is asymptotically exact in the limit of weak nonlinearity. The integrability by the spectral transform is explicitly demonstrated, because the corresponding Lax pairs have been derived, applying the same reduction method to the Lax pair of the initial equation. These systems of nonlinear PDEs are likely to be of applicative relevance and have a {open_quotes}universal{close_quotes} character, inasmuch as they may be derived from a very large classmore » of nonlinear evolution equations with a linear dispersive part. {copyright} {ital 1997 American Institute of Physics.}« less

Hybrid Upwinding for Two-Phase Flow in Heterogeneous Porous Media with Buoyancy and Capillarity

NASA Astrophysics Data System (ADS)

Hamon, F. P.; Mallison, B.; Tchelepi, H.

2016-12-01

In subsurface flow simulation, efficient discretization schemes for the partial differential equations governing multiphase flow and transport are critical. For highly heterogeneous porous media, the temporal discretization of choice is often the unconditionally stable fully implicit (backward-Euler) method. In this scheme, the simultaneous update of all the degrees of freedom requires solving large algebraic nonlinear systems at each time step using Newton's method. This is computationally expensive, especially in the presence of strong capillary effects driven by abrupt changes in porosity and permeability between different rock types. Therefore, discretization schemes that reduce the simulation cost by improving the nonlinear convergence rate are highly desirable. To speed up nonlinear convergence, we present an efficient fully implicit finite-volume scheme for immiscible two-phase flow in the presence of strong capillary forces. In this scheme, the discrete viscous, buoyancy, and capillary spatial terms are evaluated separately based on physical considerations. We build on previous work on Implicit Hybrid Upwinding (IHU) by using the upstream saturations with respect to the total velocity to compute the relative permeabilities in the viscous term, and by determining the directionality of the buoyancy term based on the phase density differences. The capillary numerical flux is decomposed into a rock- and geometry-dependent transmissibility factor, a nonlinear capillary diffusion coefficient, and an approximation of the saturation gradient. Combining the viscous, buoyancy, and capillary terms, we obtain a numerical flux that is consistent, bounded, differentiable, and monotone for homogeneous one-dimensional flow. The proposed scheme also accounts for spatially discontinuous capillary pressure functions. Specifically, at the interface between two rock types, the numerical scheme accurately honors the entry pressure condition by solving a local nonlinear problem to compute the numerical flux. Heterogeneous numerical tests demonstrate that this extended IHU scheme is non-oscillatory and convergent upon refinement. They also illustrate the superior accuracy and nonlinear convergence rate of the IHU scheme compared with the standard phase-based upstream weighting approach.

Reflection and Transmission of a Focused Finite Amplitude Sound Beam Incident on a Curved Interface

NASA Astrophysics Data System (ADS)

Makin, Inder Raj Singh

Reflection and transmission of a finite amplitude focused sound beam at a weakly curved interface separating two fluid-like media are investigated. The KZK parabolic wave equation, which accounts for thermoviscous absorption, diffraction, and nonlinearity, is used to describe the high intensity focused beam. The first part of the work deals with the quasilinear analysis of a weakly nonlinear beam after its reflection and transmission from a curved interface. A Green's function approach is used to define the field integrals describing the primary and the nonlinearly generated second harmonic beam. Closed-form solutions are obtained for the primary and second harmonic beams when a Gaussian amplitude distribution at the source is assumed. The second part of the research uses a numerical frequency domain solution of the KZK equation for a fully nonlinear analysis of the reflected and transmitted fields. Both piston and Gaussian sources are considered. Harmonic components generated in the medium due to propagation of the focused beam are evaluated, and formation of shocks in the reflected and transmitted beams is investigated. A finite amplitude focused beam is observed to be modified due to reflection and transmission from a curved interface in a manner distinct from that in the case of a small signal beam. Propagation curves, beam patterns, phase plots and time waveforms for various parameters defining the source and media pairs are presented, highlighting the effect of the interface curvature on the reflected and transmitted beams. Relevance of the current work to biomedical applications of ultrasound is discussed.

A Bifurcation Problem for a Nonlinear Partial Differential Equation of Parabolic Type,

DTIC Science & Technology

NONLINEAR DIFFERENTIAL EQUATIONS, INTEGRATION), (*PARTIAL DIFFERENTIAL EQUATIONS, BOUNDARY VALUE PROBLEMS), BANACH SPACE , MAPPING (TRANSFORMATIONS), SET THEORY, TOPOLOGY, ITERATIONS, STABILITY, THEOREMS

The nonlinear modified equation approach to analyzing finite difference schemes

NASA Technical Reports Server (NTRS)

Klopfer, G. H.; Mcrae, D. S.

1981-01-01

The nonlinear modified equation approach is taken in this paper to analyze the generalized Lax-Wendroff explicit scheme approximation to the unsteady one- and two-dimensional equations of gas dynamics. Three important applications of the method are demonstrated. The nonlinear modified equation analysis is used to (1) generate higher order accurate schemes, (2) obtain more accurate estimates of the discretization error for nonlinear systems of partial differential equations, and (3) generate an adaptive mesh procedure for the unsteady gas dynamic equations. Results are obtained for all three areas. For the adaptive mesh procedure, mesh point requirements for equal resolution of discontinuities were reduced by a factor of five for a 1-D shock tube problem solved by the explicit MacCormack scheme.

The route to chaos for the Kuramoto-Sivashinsky equation

NASA Technical Reports Server (NTRS)

Papageorgiou, Demetrios T.; Smyrlis, Yiorgos

1990-01-01

The results of extensive numerical experiments of the spatially periodic initial value problem for the Kuramoto-Sivashinsky equation. This paper is concerned with the asymptotic nonlinear dynamics at the dissipation parameter decreases and spatio-temporal chaos sets in. To this end the initial condition is taken to be the same for all numerical experiments (a single sine wave is used) and the large time evolution of the system is followed numerically. Numerous computations were performed to establish the existence of windows, in parameter space, in which the solution has the following characteristics as the viscosity is decreased: a steady fully modal attractor to a steady bimodal attractor to another steady fully modal attractor to a steady trimodal attractor to a periodic attractor, to another steady fully modal attractor, to another periodic attractor, to a steady tetramodal attractor, to another periodic attractor having a full sequence of period-doublings (in parameter space) to chaos. Numerous solutions are presented which provide conclusive evidence of the period-doubling cascades which precede chaos for this infinite-dimensional dynamical system. These results permit a computation of the length of subwindows which in turn provide an estimate for their successive ratios as the cascade develops. A calculation based on the numerical results is also presented to show that the period doubling sequences found here for the Kuramoto-Sivashinsky equation, are in complete agreement with Feigenbaum's universal constant of 4,669201609... . Some preliminary work shows several other windows following the first chaotic one including periodic, chaotic, and a steady octamodal window; however, the windows shrink significantly in size to enable concrete quantitative conclusions to be made.

Evolution of basic equations for nearshore wave field

PubMed Central

ISOBE, Masahiko

2013-01-01

In this paper, a systematic, overall view of theories for periodic waves of permanent form, such as Stokes and cnoidal waves, is described first with their validity ranges. To deal with random waves, a method for estimating directional spectra is given. Then, various wave equations are introduced according to the assumptions included in their derivations. The mild-slope equation is derived for combined refraction and diffraction of linear periodic waves. Various parabolic approximations and time-dependent forms are proposed to include randomness and nonlinearity of waves as well as to simplify numerical calculation. Boussinesq equations are the equations developed for calculating nonlinear wave transformations in shallow water. Nonlinear mild-slope equations are derived as a set of wave equations to predict transformation of nonlinear random waves in the nearshore region. Finally, wave equations are classified systematically for a clear theoretical understanding and appropriate selection for specific applications. PMID:23318680

A parallel algorithm for nonlinear convection-diffusion equations

NASA Technical Reports Server (NTRS)

Scroggs, Jeffrey S.

1990-01-01

A parallel algorithm for the efficient solution of nonlinear time-dependent convection-diffusion equations with small parameter on the diffusion term is presented. The method is based on a physically motivated domain decomposition that is dictated by singular perturbation analysis. The analysis is used to determine regions where certain reduced equations may be solved in place of the full equation. The method is suitable for the solution of problems arising in the simulation of fluid dynamics. Experimental results for a nonlinear equation in two-dimensions are presented.

Nonlinear modes of the tensor Dirac equation and CPT violation

NASA Technical Reports Server (NTRS)

Reifler, Frank J.; Morris, Randall D.

1993-01-01

Recently, it has been shown that Dirac's bispinor equation can be expressed, in an equivalent tensor form, as a constrained Yang-Mills equation in the limit of an infinitely large coupling constant. It was also shown that the free tensor Dirac equation is a completely integrable Hamiltonian system with Lie algebra type Poisson brackets, from which Fermi quantization can be derived directly without using bispinors. The Yang-Mills equation for a finite coupling constant is investigated. It is shown that the nonlinear Yang-Mills equation has exact plane wave solutions in one-to-one correspondence with the plane wave solutions of Dirac's bispinor equation. The theory of nonlinear dispersive waves is applied to establish the existence of wave packets. The CPT violation of these nonlinear wave packets, which could lead to new observable effects consistent with current experimental bounds, is investigated.

Maximum Likelihood Analysis of Nonlinear Structural Equation Models with Dichotomous Variables

ERIC Educational Resources Information Center

Song, Xin-Yuan; Lee, Sik-Yum

2005-01-01

In this article, a maximum likelihood approach is developed to analyze structural equation models with dichotomous variables that are common in behavioral, psychological and social research. To assess nonlinear causal effects among the latent variables, the structural equation in the model is defined by a nonlinear function. The basic idea of the…

Robust iterative method for nonlinear Helmholtz equation

NASA Astrophysics Data System (ADS)

Yuan, Lijun; Lu, Ya Yan

2017-08-01

A new iterative method is developed for solving the two-dimensional nonlinear Helmholtz equation which governs polarized light in media with the optical Kerr nonlinearity. In the strongly nonlinear regime, the nonlinear Helmholtz equation could have multiple solutions related to phenomena such as optical bistability and symmetry breaking. The new method exhibits a much more robust convergence behavior than existing iterative methods, such as frozen-nonlinearity iteration, Newton's method and damped Newton's method, and it can be used to find solutions when good initial guesses are unavailable. Numerical results are presented for the scattering of light by a nonlinear circular cylinder based on the exact nonlocal boundary condition and a pseudospectral method in the polar coordinate system.

Orbital stability and energy estimate of ground states of saturable nonlinear Schrödinger equations with intensity functions in R2

NASA Astrophysics Data System (ADS)

Lin, Tai-Chia; Wang, Xiaoming; Wang, Zhi-Qiang

2017-10-01

Conventionally, the existence and orbital stability of ground states of nonlinear Schrödinger (NLS) equations with power-law nonlinearity (subcritical case) can be proved by an argument using strict subadditivity of the ground state energy and the concentration compactness method of Cazenave and Lions [4]. However, for saturable nonlinearity, such an argument is not applicable because strict subadditivity of the ground state energy fails in this case. Here we use a convexity argument to prove the existence and orbital stability of ground states of NLS equations with saturable nonlinearity and intensity functions in R2. Besides, we derive the energy estimate of ground states of saturable NLS equations with intensity functions using the eigenvalue estimate of saturable NLS equations without intensity function.

Nonlinear spherical perturbations in quintessence models of dark energy

NASA Astrophysics Data System (ADS)

Pratap Rajvanshi, Manvendra; Bagla, J. S.

2018-06-01

Observations have confirmed the accelerated expansion of the universe. The accelerated expansion can be modelled by invoking a cosmological constant or a dynamical model of dark energy. A key difference between these models is that the equation of state parameter w for dark energy differs from ‑1 in dynamical dark energy (DDE) models. Further, the equation of state parameter is not constant for a general DDE model. Such differences can be probed using the variation of scale factor with time by measuring distances. Another significant difference between the cosmological constant and DDE models is that the latter must cluster. Linear perturbation analysis indicates that perturbations in quintessence models of dark energy do not grow to have a significant amplitude at small length scales. In this paper we study the response of quintessence dark energy to non-linear perturbations in dark matter. We use a fully relativistic model for spherically symmetric perturbations. In this study we focus on thawing models. We find that in response to non-linear perturbations in dark matter, dark energy perturbations grow at a faster rate than expected in linear perturbation theory. We find that dark energy perturbation remains localised and does not diffuse out to larger scales. The dominant drivers of the evolution of dark energy perturbations are the local Hubble flow and a supression of gradients of the scalar field. We also find that the equation of state parameter w changes in response to perturbations in dark matter such that it also becomes a function of position. The variation of w in space is correlated with density contrast for matter. Variation of w and perturbations in dark energy are more pronounced in response to large scale perturbations in matter while the dependence on the amplitude of matter perturbations is much weaker.

Incorporating soil variability in continental soil water modelling: a trade-off between data availability and model complexity

NASA Astrophysics Data System (ADS)

Peeters, L.; Crosbie, R. S.; Doble, R.; van Dijk, A. I. J. M.

2012-04-01

Developing a continental land surface model implies finding a balance between the complexity in representing the system processes and the availability of reliable data to drive, parameterise and calibrate the model. While a high level of process understanding at plot or catchment scales may warrant a complex model, such data is not available at the continental scale. This data sparsity is especially an issue for the Australian Water Resources Assessment system, AWRA-L, a land-surface model designed to estimate the components of the water balance for the Australian continent. This study focuses on the conceptualization and parametrization of the soil drainage process in AWRA-L. Traditionally soil drainage is simulated with Richards' equation, which is highly non-linear. As general analytic solutions are not available, this equation is usually solved numerically. In AWRA-L however, we introduce a simpler function based on simulation experiments that solve Richards' equation. In the simplified function soil drainage rate, the ratio of drainage (D) over storage (S), decreases exponentially with relative water content. This function is controlled by three parameters, the soil water storage at field capacity (SFC), the drainage fraction at field capacity (KFC) and a drainage function exponent (β). [ ] D- -S- S = KF C exp - β (1 - SFC ) To obtain spatially variable estimates of these three parameters, the Atlas of Australian Soils is used, which lists soil hydraulic properties for each soil profile type. For each soil profile type in the Atlas, 10 days of draining an initially fully saturated, freely draining soil is simulated using HYDRUS-1D. With field capacity defined as the volume of water in the soil after 1 day, the remaining parameters can be obtained by fitting the AWRA-L soil drainage function to the HYDRUS-1D results. This model conceptualisation fully exploits the data available in the Atlas of Australian Soils, without the need to solve the non-linear Richards' equation for each time-step. The spatial distribution of long term recharge and baseflow obtained with a 30 year simulation of historic data using this parameterisation, corresponds well with the spatial patterns of groundwater recharge inferred from field measurements.

Travelling-wave solutions of a weakly nonlinear two-dimensional higher-order Kadomtsev-Petviashvili dynamical equation for dispersive shallow-water waves

NASA Astrophysics Data System (ADS)

Seadawy, Aly R.

2017-01-01

The propagation of three-dimensional nonlinear irrotational flow of an inviscid and incompressible fluid of the long waves in dispersive shallow-water approximation is analyzed. The problem formulation of the long waves in dispersive shallow-water approximation lead to fifth-order Kadomtsev-Petviashvili (KP) dynamical equation by applying the reductive perturbation theory. By using an extended auxiliary equation method, the solitary travelling-wave solutions of the two-dimensional nonlinear fifth-order KP dynamical equation are derived. An analytical as well as a numerical solution of the two-dimensional nonlinear KP equation are obtained and analyzed with the effects of external pressure flow.

A comparison of linear and non-linear data assimilation methods using the NEMO ocean model

NASA Astrophysics Data System (ADS)

Kirchgessner, Paul; Tödter, Julian; Nerger, Lars

2015-04-01

The assimilation behavior of the widely used LETKF is compared with the Equivalent Weight Particle Filter (EWPF) in a data assimilation application with an idealized configuration of the NEMO ocean model. The experiments show how the different filter methods behave when they are applied to a realistic ocean test case. The LETKF is an ensemble-based Kalman filter, which assumes Gaussian error distributions and hence implicitly requires model linearity. In contrast, the EWPF is a fully nonlinear data assimilation method that does not rely on a particular error distribution. The EWPF has been demonstrated to work well in highly nonlinear situations, like in a model solving a barotropic vorticity equation, but it is still unknown how the assimilation performance compares to ensemble Kalman filters in realistic situations. For the experiments, twin assimilation experiments with a square basin configuration of the NEMO model are performed. The configuration simulates a double gyre, which exhibits significant nonlinearity. The LETKF and EWPF are both implemented in PDAF (Parallel Data Assimilation Framework, http://pdaf.awi.de), which ensures identical experimental conditions for both filters. To account for the nonlinearity, the assimilation skill of the two methods is assessed by using different statistical metrics, like CRPS and Histograms.

Dynamic Nonlinear Elastic Stability of Helicopter Rotor Blades in Hover and in Forward Flight

NASA Technical Reports Server (NTRS)

Friedmann, P.; Tong, P.

1972-01-01

Equations for large coupled flap-lag motion of hingeless elastic helicopter blades are consistently derived. Only torsionally-rigid blades excited by quasi-steady aerodynamic loads are considered. The nonlinear equations of motion in the time and space variables are reduced to a system of coupled nonlinear ordinary differential equations with periodic coefficients, using Galerkin's method for the space variables. The nonlinearities present in the equations are those arising from the inclusion of moderately large deflections in the inertia and aerodynamic loading terms. The resulting system of nonlinear equations has been solved, using an asymptotic expansion procedure in multiple time scales. The stability boundaries, amplitudes of nonlinear response, and conditions for existence of limit cycles are obtained analytically. Thus, the different roles played by the forcing function, parametric excitation, and nonlinear coupling in affecting the solution can be easily identified, and the basic physical mechanism of coupled flap-lag response becomes clear. The effect of forward flight is obtained with the requirement of trimmed flight at fixed values of the thrust coefficient.

Low-Dispersion Scheme for Nonlinear Acoustic Waves in Nonuniform Flow

NASA Technical Reports Server (NTRS)

Baysal, Oktay; Kaushik, Dinesh K.; Idres, Moumen

1997-01-01

The linear dispersion-relation-preserving scheme and its boundary conditions have been extended to the nonlinear Euler equations. This allowed computing, a nonuniform flowfield and a nonlinear acoustic wave propagation in such a medium, by the same scheme. By casting all the equations, boundary conditions, and the solution scheme in generalized curvilinear coordinates, the solutions were made possible for non-Cartesian domains and, for the better deployment of the grid points, nonuniform grid step sizes could be used. It has been tested for a number of simple initial-value and periodic-source problems. A simple demonstration of the difference between a linear and nonlinear propagation was conducted. The wall boundary condition, derived from the momentum equations and implemented through a pressure at a ghost point, and the radiation boundary condition, derived from the asymptotic solution to the Euler equations, have proven to be effective for the nonlinear equations and nonuniform flows. The nonreflective characteristic boundary conditions also have shown success but limited to the nonlinear waves in no mean flow, and failed for nonlinear waves in nonuniform flow.

On the stability of lumps and wave collapse in water waves.

PubMed

Akylas, T R; Cho, Yeunwoo

2008-08-13

In the classical water-wave problem, fully localized nonlinear waves of permanent form, commonly referred to as lumps, are possible only if both gravity and surface tension are present. While much attention has been paid to shallow-water lumps, which are generalizations of Korteweg-de Vries solitary waves, the present study is concerned with a distinct class of gravity-capillary lumps recently found on water of finite or infinite depth. In the near linear limit, these lumps resemble locally confined wave packets with envelope and wave crests moving at the same speed, and they can be approximated in terms of a particular steady solution (ground state) of an elliptic equation system of the Benney-Roskes-Davey-Stewartson (BRDS) type, which governs the coupled evolution of the envelope along with the induced mean flow. According to the BRDS equations, however, initial conditions above a certain threshold develop a singularity in finite time, known as wave collapse, due to nonlinear focusing; the ground state, in fact, being exactly at the threshold for collapse suggests that the newly discovered lumps are unstable. In an effort to understand the role of this singularity in the dynamics of lumps, here we consider the fifth-order Kadomtsev-Petviashvili equation, a model for weakly nonlinear gravity-capillary waves on water of finite depth when the Bond number is close to one-third, which also admits lumps of the wave packet type. It is found that an exchange of stability occurs at a certain finite wave steepness, lumps being unstable below but stable above this critical value. As a result, a small-amplitude lump, which is linearly unstable and according to the BRDS equations would be prone to wave collapse, depending on the perturbation, either decays into dispersive waves or evolves into an oscillatory state near a finite-amplitude stable lump.

Nonlinear ordinary difference equations

NASA Technical Reports Server (NTRS)

Caughey, T. K.

1979-01-01

Future space vehicles will be relatively large and flexible, and active control will be necessary to maintain geometrical configuration. While the stresses and strains in these space vehicles are not expected to be excessively large, their cumulative effects will cause significant geometrical nonlinearities to appear in the equations of motion, in addition to the nonlinearities caused by material properties. Since the only effective tool for the analysis of such large complex structures is the digital computer, it will be necessary to gain a better understanding of the nonlinear ordinary difference equations which result from the time discretization of the semidiscrete equations of motion for such structures.

Implementation of an integrated op-amp based chaotic neuron model and observation of its chaotic dynamics

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

Jung, Jinwoo; Lee, Jewon; Song, Hanjung

2011-03-15

This paper presents a fully integrated circuit implementation of an operational amplifier (op-amp) based chaotic neuron model with a bipolar output function, experimental measurements, and analyses of its chaotic behavior. The proposed chaotic neuron model integrated circuit consists of several op-amps, sample and hold circuits, a nonlinear function block for chaotic signal generation, a clock generator, a nonlinear output function, etc. Based on the HSPICE (circuit program) simulation results, approximated empirical equations for analyses were formulated. Then, the chaotic dynamical responses such as bifurcation diagrams, time series, and Lyapunov exponent were calculated using these empirical equations. In addition, we performedmore » simulations about two chaotic neuron systems with four synapses to confirm neural network connections and got normal behavior of the chaotic neuron such as internal state bifurcation diagram according to the synaptic weight variation. The proposed circuit was fabricated using a 0.8-{mu}m single poly complementary metal-oxide semiconductor technology. Measurements of the fabricated single chaotic neuron with {+-}2.5 V power supplies and a 10 kHz sampling clock frequency were carried out and compared with the simulated results.« less

New Quantum Diffusion Monte Carlo Method for strong field time dependent problems

NASA Astrophysics Data System (ADS)

Kalinski, Matt

2017-04-01

We have recently formulated the Quantum Diffusion Quantum Monte Carlo (QDMC) method for the solution of the time-dependent Schrödinger equation when it is equivalent to the reaction-diffusion system coupled by the highly nonlinear potentials of the type of Shay. Here we formulate a new Time Dependent QDMC method free of the nonlinearities described by the constant stochastic process of the coupled diffusion with transmutation. As before two kinds of diffusing particles (color walkers) are considered but which can further also transmute one into the other. Each of the species undergoes the hypothetical Einstein random walk progression with transmutation. The progressed particles transmute into the particles of the other kind before contributing to or annihilating the other particles density. This fully emulates the Time Dependent Schrödinger equation for any number of quantum particles. The negative sign of the real and the imaginary parts of the wave function is handled by the ``spinor'' densities carrying the sign as the degree of freedom. We apply the method for the exact time-dependent observation of our discovered two-electron Langmuir configurations in the magnetic and circularly polarized fields.

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

NASA Astrophysics Data System (ADS)

Calini, A.; Schober, C. M.

2013-09-01

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