The study of the properties of fluids and gases in motion

learn more… | top users | synonyms (1)

1
vote
0answers
8 views

OpenFOAM precipitation/crystallisation reaction solver

This question is probably related to chemical engineers that are around this forum. I am looking into writing a solver for a precipitation reaction (struvite to be precise): A+B+C+D <=> E+F A, B, ...
0
votes
0answers
27 views

How to get started with MATLAB? [on hold]

How to get started with MATLAB specially for Numerical Analysis and Computational Fluid Dynamics?
1
vote
2answers
38 views

Estimation of time taken to reach steady state in an MD simulation of Poiseuille flow

I am trying to do a Molecular Dynamics Simulation of a complex fluid, confined between solid surfaces. I would like to find the flow rate as a function of fluid film thickness, $h$, for a plane ...
2
votes
0answers
41 views

Implementation of no-slip boundary conditions in lattice Boltzmann method fluid simulation

My faculty advisor recommended that I take a look at the lattice Boltzmann method as an introduction to scientific computing and potentially an undergraduate honors thesis topic. I cooked up a some ...
3
votes
2answers
48 views

Adding artificial fluctuations to trigger hydrodynamic instabilities.

I'm doing a 2D combustion hydrodynamic simulation and there's a hydrodynamic instability that should be triggered because of the particular physical properties of my system. The key to the instability ...
4
votes
2answers
129 views

How to discretize Burger's equation?

I am trying to solve the very simple one dimensional burgers equation which is: $$\frac{\delta U}{\delta t} + \frac{\delta F}{\delta x} = 0$$ where the flux F of some variable U is defined as$$ F= ...
0
votes
0answers
25 views

Tube with cylinder inside in OpenFoam

There is a solid irregularly shaped cylinder inside a tube. Fluid flows through the irregularly shaped tube and the cylinder. In other words, the radii of tube and cylinder are not constant. One end ...
0
votes
2answers
73 views

3D Divergence-free Stokes equations

I would like to have a divergence-free formulation of the Navier-Stokes equations for creeping flow (a.k.a. Stokes equation, N.-S. without the inertial terms) in 3D, for purposes of tracking particle ...
3
votes
0answers
64 views

Is it generally unstable to use in multidimensional simulations finite difference schemes with higher orders than 2?

I'm part of a team trying to generalize a 1D Advection-Diffussion-Reaction code we inherited by extending it to 2D by using dimensional splitting, i.e. solving advection and diffusion for x and y ...
5
votes
2answers
109 views

9-point stencil finite difference Laplacian with variable diffusion coefficients

So I'm trying to implement a 9-point stencil discretization to the 2D difussion equation. The stencil is here. However, most of the literature deals with a Laplacian that has a constant diffusion ...
4
votes
1answer
66 views

transverse component for multidimensional advection in method of lines

So I inherited from some people a code that solves the advection-diffusion-reaction equation for a particular system. The original code was first implemented in 1D which worked fine in cartesian ...
2
votes
1answer
114 views

Darcy flow finite elements

The Darcy equations for porous media flow are given by: $\frac{\mu}{\kappa}\mathbf{u} - \nabla p = \mathbf{0}$ $\nabla\cdot\mathbf{u} = 0$ where $\kappa$ is the permeability and can in general be ...
1
vote
0answers
30 views

ANSYS Fluent: Defining External Force on a Cell Zone [closed]

I am using ANSYS Fluent 15 to simulate a case of fluid flow. In my case, I have three cell zones in my simulation, let's call them top, middle and bottom. My requirement is to put a vertical constant ...
2
votes
1answer
69 views

CFD turbulence modelling mean pressures vs peak pressures

I have been experimenting with using Autodesk CFD to investigate facade/ cladding pressures on a (rectangular) building, comparing results with cladding/ facades pressures pressures from design codes ...
2
votes
1answer
79 views

What are acceptable boundary conditions for porous media flow?

I am attempting to simulate fluid flow through a porous foam. I would like to have no-slip boundary conditions on part of the boundary and free flow conditions on the inlet and outlet. Right now I am ...
0
votes
0answers
37 views

Computing Kolmogorov/Energy spectrum

I read a comment that you calculated energy spectrum for 3D DNS data obtained for isotropic turbulence which is equally spaced in all three directions. I have the same case, I am trying to do this ...
3
votes
3answers
188 views

Where do I begin learning CFD?

I am working on some CFD codes - generating adjoint of a CFD model, for example. My background however is in computer science. I have absolutely no clue what the model represents or what it is doing. ...
0
votes
1answer
51 views

Visualise data under Tecplot

I have a set of data that I need to visualise under Tecplot, the file format is for example: ...
0
votes
0answers
40 views

modify COMSOL codes

I am a physics student and I am planning to simulate microfluidics and capillary waves on a fluid film and I decided to use COMSOL. I am not so much aware of COMSOL abilities. I wonder if I can ...
1
vote
0answers
50 views

Conservation at grid interface in adaptive mesh refinement

I am using adaptive mesh refinement to solve one dimensional inviscid Burgers equation. However I am facing some difficulty to handle grid interfaces which are not uniform (coarse-fine grid ...
1
vote
1answer
94 views

openfoam - Programming customized PDEs

I am looking for a method to automatically solve custom PDEs on a custom control volume. Specifically I would like to solve equations similar (but not exactly alike) to: $$\frac{\partial y}{\partial ...
1
vote
1answer
101 views

GPU-enabled Lattice Boltzman solvers?

Is anybody aware of any GPU-enabled Lattice Boltzman solvers (preferably on C++/OpenCL and open-source) that would be recommended? I have found Advanced Simulation Library, but it seems to be very ...
0
votes
0answers
19 views

Plotting two different results in Tecplot together

I want to merge two results together and display them in Tecplot. For example: I have one DAT file which contains temperature of 64 points in each x and y direction and the other contains temperature ...
1
vote
2answers
131 views

How to deal with nonlinear term in Navier Stokes equations (finite element code)

I am trying to solve the Navier Stokes equations using the finite element method. I plan on using the pressure correction method to deal with the pressure and an implicit time stepping scheme for ...
0
votes
1answer
101 views

Adaptive mesh refinement basic conceptual problem

I am a beginner in adaptive mesh refinement (AMR). After I am done with the first two papers by Dr. Marsha Berger, I was trying to write my own code for a problem which has a parabolic partial ...
0
votes
1answer
52 views

Representation of polynomial order in CFD codes

I am currently working on a CFD code over a cubic grid. Now, the number of elements used in the simulation is decomposed among the number of processors. Each of those processors (a section of the ...
0
votes
4answers
229 views

Why are upwind schemes stable in convection flow calculation?

It is well known that upwind schemes are stable when calculating convection flows with $|\text{Pe}|>2$, $\text{Pe}$ is the Peclet number. Why is that, and why is central difference unstable? Is ...
1
vote
1answer
110 views

A Question About the Rhie-Chow Interpolation Used for Solving the Incompressible Navier-Stokes Equations on Unstructured Grids

When using the SIMPLE method on a mesh with a collocated variable arrangement, the following interpolation is used for the advecting velocities: \begin{equation} u_f = \overline{u}_f - ...
2
votes
2answers
112 views

Direct Numerical Simulation of reacting flows

I need to write a DNS code for simulating reacting flows in turbulent conditions. The code has to be highly scalable, because the computational cost of simulations is expected very large. My idea is ...
3
votes
1answer
100 views

2D Poisson Solver for Taylor Green Vortex Problem

I am trying to write a 2D Navier Stokes solver using an RK3 for time advancement in python. For debugging, I have converted the RK3 to an Euler step for simplicity. Checking my divergence for my ...
0
votes
0answers
47 views

Best way to average surface data or function data in 2d grid

I have pressure data on a 2D triangular surface. I can compute average(data). What is the best way to average the surface data?
5
votes
2answers
191 views

applications of computational geometry in fields such as CFD?

Out of curiosity, I was recently trying to search what skills are required to be successful as developer in scientific computing field (e.g. CFD or similar). And to do so, I was going to through ...
1
vote
1answer
77 views

alternatives to moving mesh technique

If one has to simulate oscillating plate/solid wall ( sinusoidal function of time) in a domain (a simple piston movement in 'y' direction), the obvious way would be scenario 1. For brevity, assume ...
1
vote
3answers
314 views

2D Stokes equation Code

Does anyone know where could I find a code (in Matlab or Mathematica, for example) for he Stokes equation in 2D? It has been solved numerically by so many people and referenced in so many paper that I ...
2
votes
1answer
144 views

Add User-defined/custom differential equations in OpenFoam (CFD)

I am new to OpenFoam. And I am trying to add a set (user defined) of differential equations to OpenFoam. I want to solve this user defined set of equations at each time point in addition to standard ...
0
votes
0answers
48 views

Can COMSOL solve a fluids problem given complicated N-S forcing terms as the input?

I am very much a novice at COMSOL but after pouring over its documentation, I have yet to build a working fluid model with this method. Here's the background: I am attempting to model a turbulent jet ...
0
votes
0answers
47 views

Help about Fluid-Fluid coupling techniques

I need a little help and advice with a project I want to do: the idea is to "couple" (I don't know whether I can call it like this) a conservative Navier-Stokes Solver (Fractional-step, 2nd order FDM) ...
0
votes
2answers
46 views

The “standard map” and Jupiter's red spot [closed]

Has there been an attempt to "explain" Jupiter's Red Spot via what is known as the standard map in chaos theory? I cannot be the first to notice the visual similarities...       ...
0
votes
0answers
23 views

Computational time in fluent using PBM+REACTION

I am having a problem of huge computation time. My system is a reaction, nucleation and growth. Time step i have taken is 10^-6 sec. it takes 10 minutes to reach 10^-2 sec. It will nearly take 70 ...
1
vote
2answers
112 views

How to interpret/convert density in 2D simulations?

Density has units of $kg/m^3$. In simulations often scaling is done to deal with non-dimensional values, so the physical density is converted to the non-dimensional density so that it can be used in ...
1
vote
1answer
46 views

Value of density when there are no or very few neighbours in SPH simulation

Sorry for the noob question. I am trying to implement SPH using the directions shown in this paper. The density needs to be updated using the formula The smoothing kernel is If there are no ...
2
votes
2answers
264 views

Manufactured solution for pressure based 3d incompressible Navier-Stokes solver with wall boundaries

I already successfully verified my solver (SIMPLE-type FVM-method) with the following manufactured solution (3d Taylor-Green vortex) on the solution domain $[-1,1]^3$ with Dirichlet boundary ...
0
votes
1answer
186 views

How to solve Energy Balance equation by numerical method

Good Day I am new to heat transfer technique please give me some suggestion on solving energy balance equation $$a \frac{\partial T_p}{\partial t}=\frac{\partial}{\partial x}\left(b\frac{\partial ...
1
vote
2answers
123 views

How to prove the strong formulation for Discontinuous Galerkin is conservative?

I've been trying to figure out how to prove that the following equation is conservative $$\int_{D^k} \partial_t(u^k) v_j^k + \partial_x(f) v_j^k dx = [(f -f^*) v^k_j]_{x^k}^{x^{k+1}}$$ where $v_j^k$ ...
2
votes
0answers
90 views

Arbitrary Choosing of the Solution Domain - Navier Stokes and Manufactured Solutions

I want to verify a finite-volume solver (SIMPLE-Algorithm) for the incompressible Navier-Stokes equations by using a manufactured solution. I use Dirichlet boundary conditions for the velocity at all ...
4
votes
3answers
145 views

Variable viscosity Stokes equation

One very efficient way to solve Stokes equation with periodic boundary conditions \begin{equation} -\eta \nabla^{2} \bf{v} + \nabla p = f \\ \nabla \cdot \bf{v} = 0 \end{equation} is using the ...
3
votes
2answers
377 views

Solution blows up when using Runge-Kutta to solve simultaneous ODEs for liquid film equations

This is a continuation of a problem I asked over at physics exchange and math exchange. Basically I have two ODEs that I am solving in order to calculate the radial and tangential velocity of liquid ...
2
votes
1answer
223 views

What's the difference between grid-based and mesh-based methods for PDEs?

I am into computational fluid dynamics and so far I've found that the most common approaches to solve for the governing equations are Eulerian and Lagrangian. The former samples the domain at fixed ...
3
votes
1answer
64 views

Zero-k mode in Pseudo-spectral solution of Stokes Flow

I'm trying to solve a Stokes flow problem with a pseudo-spectral method in periodic boundary conditions. The equations of interest are $-\nabla^2 \bf{v} + \nabla p = \bf{f} \\ \nabla \cdot \bf{v} = ...
5
votes
1answer
155 views

ENO/WENO component-wise vs characteristic-wise

Can someone give some references to understand what's the differences between a component-wise and a characteristic-wise ENO scheme? If I'm right, the characteristic variables come from the ...