Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Semi-implicit”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

The application of Seasat-1 radar altimetry to continental shelf circulation modeling

Seasat-1 radar altimetry data have been used to verify a continental shelf circulation model subjected to tidal forcing on the open boundary. The model makes use of the semi-implicit mode of time integration, removing the surface gravity wave time step restriction. Both the altimetry and the model predictions are compared with empirically determined tidal fluctuations and generally good agreement is obtained.

Cornillon, P.↗

The solution of variable-step implicit difference equations for dynamic systems analysis

A semi-implicit iterative technique is given for the solution of implicit matrix difference equations which arise in variable-step direct time integration of dynamic systems equations. For all practical purposes, the technique avoids repeated refactoring of the stepsize-dependent solution matrix and achieves the solution convergence for engineering accuracy requirements within a permissible number of iterations. Extension of the present method to nonlinear problems treated by the pseudo-force procedure is straightforward. For nonlinear solution procedures based on the tangent modulus formulation, the technique requires modification to make use of the tangent stiffness matrix. It is noted that the question of whether to include the nonlinear terms within the iteration loop has not yet been resolved.

Park, K. C.↗

Laminar flow transition: A large-eddy simulation approach

A vectorized, semi-implicit code was developed for the solution of the time-dependent, three dimensional equations of motion in plane Poiseuille flow by the large-eddy simulation technique. The code is tested by comparing results with those obtained from the solutions of the Orr-Sommerfeld equation. Comparisons indicate that finite-differences employed along the cross-stream direction act as an implicit filter. This removes the necessity of explicit filtering along this direction (where a nonhomogeneous mesh is used) for the simulation of laminar flow transition into turbulence in which small scale turbulence will be accounted for by a subgrid scale turbulence model.

Biringen, S.↗

Iterative spectral methods and spectral solutions to compressible flows

A spectral multigrid scheme is described which can solve pseudospectral discretizations of self-adjoint elliptic problems in O(N log N) operations. An iterative technique for efficiently implementing semi-implicit time-stepping for pseudospectral discretizations of Navier-Stokes equations is discussed. This approach can handle variable coefficient terms in an effective manner. Pseudospectral solutions of compressible flow problems are presented. These include one dimensional problems and two dimensional Euler solutions. Results are given both for shock-capturing approaches and for shock-fitting ones.

Hussaini, M. Y.↗

A spectral collocation method for the Navier-Stokes equations

A Fourier-Chebyshev spectral method for the incompressible Navier-Stokes equations is described. It is applicable to a variety of problems including some with fluid properties which vary strongly both in the normal direction and in time. In this fully spectral algorithm, a preconditioned iterative technique is used for solving the implicit equations arising from semi-implicit treatment of pressure, mean advection and vertical diffusion terms. The algorithm is tested by applying it to hydrodynamic stability problems in channel flow and in external boundary layers with both constant and variable viscosity.

Malik, M. R.↗

A spectral collocation method for the Navier-Stokes equations

A Fourier-Chebyshev spectral method for the incompressible Navier-Stokes equations is described. It is applicable to a variety of problems including some with fluid properties which vary strongly both in the normal direction and in time. In this fully spectral algorithm, a preconditioned iterative technique is used for solving the implicit equations arising from semi-implicit treatment of pressure, mean advection and vertical diffusion terms. The algorithm is tested by applying it to hydrodynamic stability problems in channel flow and in external boundary layers with both constant and variable viscosity.

Malik, M. R.↗

The efficient calculation of chemically reacting flow

A semi-implicit finite volume formulation is used to study flows with chemical reactions. In this formulation the source terms resulting from the chemical reactions are treated implicitly and the resulting system of partial differential equations is solved using two time-stepping schemes. The first is based on the Runge-Kutta method while the second is based on an Adams predictor-corrector method. Results show that improvements in computational efficiency depend to a large extent on the manner in which the source term is treated. Further, analysis and computation indicate that the Runge-Kutta method is more efficient than the Adams methods. Finally, an adaptive time stepping scheme is developed to study problems involving shock ignition. Calculations for a hydrogen-air system agree well with other methods.

Eklund, D. R.↗

A relaxation technique for the parabolized Navier-Stokes (PNS) equations

A rapidly converging relaxation technique for the parabolized Navier-Stokes equations has been devised. The scheme is applicable in both supersonic and subsonic flows, but it is discussed here in the context of supersonic flows. The upstream propagating acoustic influence in the subsonic part of the flow is introduced semi-implicitly through the streamwise momentum equation applied on the body, and through a forward-differencing on the streamwise pressure gradient term in the interior. This procedure yields a new boundary condition on the energy in the total energy equation. The pressure-velocity system in the subsonic layer is coupled, but the positive time-like marching characteristic of the governing equations is still maintained. The relaxation technique is demontrated to work for a three-dimensional flow over a cone-flare in supersonic flight.

Kaul, Upender K.↗

Navier-Stokes cascade analysis with a stiff Kappa-Epsilon turbulence solver

The two dimensional, compressible, thin layer Navier-Stokes equations with the Baldwin-Lomax turbulence model and the kinetic energy-energy dissipation (k-epsilon) model are solved numerically to simulate the flow through a cascade. The governing equations are solved for the entire flow domain, without the boundary layer assumptions. The stiffness of the k-epsilon equations is discussed. A semi-implicit, Runge-Kutta, time-marching scheme is developed to solve the k-epsilon equations. The impact of the k-epsilon solver on the explicit Runge-Kutta Navier-Stokes solver is discussed. Numerical solutions are presented for two dimensional turbulent flow over a flat plate and a double circular arc cascade and compared with experimental data.

Liu, Jong-Shang↗

Studies of pressure-velocity coupling schemes for analysis of incompressible and compressible flows

Two pressure-velocity coupling schemes, both of which solve the fully implicit discretized equations governing the flow of fluids were examined, and the capability of performing large Reynolds number, low Mach number compressible flow calculations were assessed. The semi-implicit iterative SIMPLE algorithm is extended to handle transient compressible flow calculations. This extension takes into account a strong coupling between the pressure and temperature through a correction procedure, based on the equations of state. Results obtained from the extended SIMPLE algorithm are then compared to similar results obtained from the non-iterative PISO algorithm. Both time-dependent and steady state calculations were performed using an axisymmetric 2:1 pipe expansion geometry and laminar flow conditions corresponding to Reynolds number of 1000 and Mach number of 2.0. For calculations simulating a time-dependent compression/expansion process, both schemes exhibit transient features in excellent agreement with each other, and moreover, the PISO method shows a significant computational time reduction of 60 percent over the SIMPLE scheme, regardless of the time step size or grid size employed. The effects of numerical diffusion are shown to be significant in these calculations. For steady state compressible flows, however, the SIMPLE algorithm displays increasing computational efficiency over the PISO method as the time step sizes employed to reach steady state conditions are decreased.

Burgreen, Gregory Wayne↗

Navier-Stokes cascade analysis with a stiff k-epsilon turbulence solver

The two dimensional, compressible, thin layer Navier-Stokes equations with the Baldwin-Lomax turbulence model and the kinetic energy-energy dissipation (k-epsilon) model are solved numerically to simulate the flow through a cascade. The governing equations are solved for the entire flow domain, without the boundary layer assumptions. The stiffness of the k-epsilon equations is discussed. A semi-implicit, Runge-Kutta, time-marching scheme is developed to solve the k-epsilon equations. The impact of the k-epsilon solver on the explicit Runge-Kutta Navier-Stokes solver is discussed. Numerical solutions are presented for two dimensional turbulent flow over a flat plate and a double circular arc cascade and compared with experimental data.

Liu, Jong-Shang↗

Semi-Lagrangian integration of a gridpoint shallow water model on the sphere

A stable, semi-Lagrangian, semi-implicit, two-time-level, gridpoint integration scheme for the shallow water equations on the sphere is presented. A rotated spherical coordinate system is used to integrate the equations of motion at each gridpoint poleward of a certain latitude, thus overcoming problems associated with the polar singularity. The results of medium term integrations of large scale test patterns using a long time step are presented.

Mcdonald, A.↗

Numerical experiments on transition control in wall-bounded shear flows

Results are presented from a numerical simulation of transition control in plane channel and boundary layer flows. The analysis is based on a pseudo-spectral/finite difference semi-implicit solution procedure employed to numerically integrate the time-dependent, three-dimensional, incompressible Navier-Stokes equations in a doubly periodic domain. In the channel flow, it was found that the active periodic suction/blowing method was effective in controlling strongly three-dimensional disturbances. In the boundary layer, the preliminary analysis indicated that in the early stages, passive control by suction is as effective as active control to suppress instabilities. The current work is focused on a detailed comparison of active and passive control by suction/blowing in the boundary layer.

Biringen, S.↗

Oscillatory flow with heat transfer in a square cavity

A computational study is presented for the flow inside an oscillatory cavity. The numerical scheme employs a semi-implicit, time-splitting method to integrate the two-dimensional full Navier-Stokes equations satisfying continuity to machine accuracy. The efficient use of direct solvers for the uncoupled momentum and pressure equations is demonstrated. The oscillatory cavity flow is studied considering the effects of heat transfer, Reynolds number, and oscillatory Stokes number.

Biringen, S.↗

Numerical simulation of 3-D Benard convection with gravitational modulation

In this numerical study, randomly and sinusoidally modulated gravitational fields imposed on three-dimensional Rayleigh-Benard convection are investigated in an effort to understand the effects of vibration (G-Jitter) on fluid systems. The time-dependent, Navier-Stokes equations and the energy equation with Boussinesq approximations are solved by a semi-implicit, pseudospectral procedure. An analysis of energy balances indicates that with increasing modulation amplitude, transition from synchronous to relaxation oscillation goes through the subharmonic response. Random modulations are found to be less stabilizing than sinusoidal and are shown to impose three-dimensionality on the flow for some parameter ranges both at terrestrial and zero base gravity conditions.

Biringen, S.↗

Integration of the shallow water equations on the sphere using a vector semi-Lagrangian scheme with a multigrid solver

A vector semi-Lagrangian semi-implicit two-time-level finite-difference integration scheme for the shallow water equations on the sphere is presented. A C-grid is used for the spatial differencing. The trajectory-centered discretization of the momentum equation in vector form eliminates pole problems and, at comparable cost, gives greater accuracy than a previous semi-Lagrangian finite-difference scheme which used a rotated spherical coordinate system. In terms of the insensitivity of the results to increasing timestep, the new scheme is as successful as recent spectral semi-Lagrangian schemes. In addition, the use of a multigrid method for solving the elliptic equation for the geopotential allows efficient integration with an operation count which, at high resolution, is of lower order than in the case of the spectral models. The properties of the new scheme should allow finite-difference models to compete with spectral models more effectively than has previously been possible.

Bates, J. R.↗

A Fourier-Chebyshev pseudospectral method for solving steady 3-D Navier-Stokes and heat equations in cylindrical cavities

A Fourier-Chebyshev pseudospectral method for solving steady 3D Navier-Stokes equations in cylindrical cavities is presented and discussed. The general method is pseudo-unsteady and uses a semi-implicit finite difference scheme for the time discretization. The generalized ADI (Alternating Direction Implicit) procedure is then applied to reduce the problem to a successive solution of one-dimensional problems. The spatial approximation uses a Fourier-Galerkin approximation in the periodic direction and a Chebyshev-collocation approximation in the other directions. Difficulties related to the pressure are surmounted by using the artifical compressibility method. A suitable variable change was chosen to avoid the problem of singularity at the axis generated by cylindrical coordinates. The method is first tested on an advection-diffusion equation and then on the Navier-Stokes equations. Finally, the method is illustrated by a convection problem of a differentially heated fluid.

Pulicani, J. P.↗

A Fourier collocation time domain method for numerically solving Maxwell's equations

A new method for solving Maxwell's equations in the time domain for arbitrary values of permittivity, conductivity, and permeability is presented. Spatial derivatives are found by a Fourier transform method and time integration is performed using a second order, semi-implicit procedure. Electric and magnetic fields are collocated on the same grid points, rather than on interleaved points, as in the Finite Difference Time Domain (FDTD) method. Numerical results are presented for the propagation of a 2-D Transverse Electromagnetic (TEM) mode out of a parallel plate waveguide and into a dielectric and conducting medium.

Shebalin, John V.↗