Carlos Andres Rivera Villarreyes Global Product Specialist - FEFLOW

Dear Joachim,

Please notice that FePEST 6.1 was a quasi beta version compared to FePEST 6.2, which is fully developed. I strongly recommend you to use 6.2 version.
What exactly would you like to do? Could you be more especific in what kind of regularization setings you need to chance?

Cheers,

Carlos

Hi Debora and Adam,

You can take all the advantage of the graphical interface of the IFM plug-in. Below a small example to create two input boxes. You can reach this by right-click on the name of attached plug-in and then choose [b]Edit [/b] option.

[color=blue][i]static IfmPropExtInt aux= { 2, 50};

IfmProperty props[] = {
{ "Param1", IfmPROP_DOUBLE, &par1, &DefaultParam1, "Write value of Parameter 1", &aux },
{ "Param2", IfmPROP_DOUBLE, &par2, &DefaultParam2, "Write value of Parameter 2", &aux },
}

IfmEditProperties (pDoc, "Properties of IFM Plugin", "My parameters: ", props, 2); // last number is amount of parameters.
[/i][/color]

The best place of the code is on [b]OnEditDocument[/b]. Don't forget to declare the default values (&DefaultParam1) and the variables for the parameter itself (&par1).

More sophisticated manners are using fully Qt stuff (i.e. design your own *.ui file).

Cheers,

Carlos

Dear Debora,

I have in mind different possibilities to reach your goal:

1) Define several parameter definitions with different parameter bounds. For example, you can assign bounds with a decrease in conductivity over depth.

2) If the ratio is Kx/Kz is kept over the layer, you can provide several ratios which change over depth. Here you will need to use tied parameter when you create the definitions (e.g. Kz will be tied to Kx). The parameter is computed from the initial values assigned in the FEM file. FePEST will keep this ratio over the entire optimization.

3) Use Prior Information. The syntax of Prior Information lines in FePEST follows exactly the PEST description (see PEST Manual, section 4.2.10). An example below:

PIBL PIFAC * PARNAME + ... = PIVAL WEIGHT OBGNME

where PIBL is a character variable of the Prior Information Definition, PIFAC is a factor which affects the parameter, PARNME is the parameter name, PIVAL is a real variable to define the right-side of the PI equation, WEIGHT is the weight assigned to the article of prior information in the parameter estimation process (e.g. inverse proportional to standard deviation), and OBGNME is the observation group to which the prior observation belongs.

FePEST separates the different variables in different entry boxes for better visualization, you need just to define the equation part. Let's say you have conductivity values 100, 50, 25 and 5 decreasing over depth. In FePEST you can use PI equation to suggest that conductivity needs to decrease in depth during the optimization as follows:

1.0 * log(Kxx_1) = 2.00
1.0 * log(Kxx_2) = 1.70
1.0 * log(Kxx_3) = 1.40
1.0 * log(Kxx_4) = 0.70

Notice that factor 1.0 is required. If the parameter is log-transformed, you need to indicate it in the equation.

Hi,

Just few comments based on the previous discussion:

1) Are you running a Richards-based model with a steady-state solution? If so, I would highly recommend to use the transient mode. It is very hard to solve the non-linearity of the Richards' Equation under a steady-state solution.

2) It seems to me that model is solving fine. The fact that you see a "second" water table at the top, it reflects the typical wetting front propagation in an unsaturated medium. Ideally over the time, if recharge decreases, you may see the movement of this "second" phreatic level down to the regional water level. If the recharge is very large, eventually you will see that model becomes completely saturated after several iterations, therefore it will reach a topography-like phreatic surface.

3) Material properties: Are you sure that such van Genuchten values correspond to the assigned conductivity? As indicated in item 2, it seems to me that water is moving very slow. A higher conductivity necessarily does mean fast water movement, if van Genuchten properties are not the adequate. So, please make double-check.

4) You may try to decrease the non-linearity by means of using the Modified Van Genuchten relationship and assuming a linear form (i.e. what you named before as "sigma").

Two more alternatives (in case analysis should be done varying material properties):

1) Use FePEST in the command line mode to setup a SENSAN (Sensitivity Analysis) run using PEST. Here you can define the parameter to be modified and distribute the output results in different files. It may require some time until to get into the workflow, but then it is worthy since you can have additional feedbackfrom PEST, which may help you to understand the sensitivity analysis.

2) Download the old FEFLOW-specific groundwater utilities for PEST. There is a nice routine named FMAINFEM to assign automatically material values to a FEM file directly in the windows command line. This was used in the past to run PEST with FEFLOW, now with the FePEST interface is not required anymore. Link is http://www.feflow.com/miscellaneos.html

Sure. You need to define adequate parameters for the unsaturated flow modelling carried out by the Richards equation. This means you need an empirical model to relate saturation and negative pressure, saturation limits, etc. By the way, analogous to the residual water depth in the Phreatic approach, you will also need to setup an adequate residual saturation limit proper for the material properties.

Hi Michael,

"Pending" means that FEFLOW cannot retrieve primary variables values to compute the budget. You just need to start the simulation.
If you save the DAC file using the Record Properties option at the simulation toolbar, you can have all the entire model results to compute any kind of budget.

Cheers,

Carlos

There are two alternatives:

1) Create Observation Points at specific locations and connect them to the reference values. Then go to View - Chart - Hydraulic Head History and you will have the simulated and "reference" heads. Scatter plot is also available. You can find a lot of information how to add Observation Points in the Help System.

2) Export the FEFLOW simulated hydraulic heads. Simply a right-click on Hydraulic Head and the Export Data. You can export values from either all nodes or nodal selection.

Hi,

The slices below should be indicated as "Dependent". The "free" first slice will move according to the fluctuations of the phreatic level. Here, pressure will be always zero and FEFLOW will adjust the slice elevation to fulfill the statement. The next slices are adapted to changes of free slice.
I would recommend you to take a look on White Paper I (http://www.feflow.info/manuals.html) for further details.

Cheers,

Carlos

Dear Perl,

Personally I believe there is not an exact rule for deciding about interpolation techniques. The selection depends on many variables such as data type, data transformation, distribution of your points, geometry of domain, etc. These last two are the most critical. I think the best way is a trial-and-error approach here. You can make double-check plotting both interpolated fields and measurements at the same view.

Regards,

Carlos